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ABSTRACT 


A  genetic  study  was  initiated  in  the  tropical  legume 
Pesmodlum  to  establish  whether  there  existed  genetic 
variation  in  the  symbiotic  relationship  between  I>e£raodi.um 
and  Rh lzoblum>  Particular  emphasis  was  placed  upon  effective 
vs •  ineffective  nodulation  and  quantitative  variation  for 
the  estimates  of  nitrogen  fixed* 

Analysis  of  F2  and  backcross  populations  from  two 
interspecific  crosses  suggested  that  a  strain-specific 
ineffective  nodulation  interaction  was  conditioned  by  a 
dominant  allele  at  a  single  locus* 

Significant  genetic  variation  was  found  among  genotypes 
for  all  nitrogen  fixation  estimates*  The  major  portion  of 
the  variation  among  genotypes  was  non-additive  variation  and 
was  attributable  to  a  significant  Varieties  vs*  Crosses 
genetic  component*  The  mean  of  the  F 1  crosses  was  found  to 
be  80%  larger  than  the  mean  of  the  inbred  parents  Indicating 
that  most  F*  crosses  exceeded  their  larger  parent  for  the 
quantitative  estimation  of  nitrogen  fixation*  Since 
inbreeding  depression  for  forage  yield  is  generally 


considered 

to 

be 

absent 

in  the  presence  of 

adequate 

nitrogen, 

the 

Fi 

hyb  ri d 

vigor  for  the  nitrogen 

f ixat 1  on 

estimates  was  due  solely  to  these  differences  when  the 
plants  were  grown  in  a  nitrogen  limiting  environment*  The 
significance  of  these  results  was  discussed* 
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Large  cotyledons  and  juvenile  leaves  appeared  to  be 
controlled  by  dominant  genes  at  different  loci*  Plants  with 
large  cotyledons  or  juvenile  leaves  fixed  significantly  more 
nitrogen  than  plants  with  small  ones*  Adjusting  the 
quantitative  estimates  for  differences  due  to  cotyledon  and 
juvenile  leaf  size  produced  a  totally  additive  model  for  the 
inheritance  of  the  estimates  of  nitrogen  fixation*  The  non¬ 
additive  variance  for  the  unadjusted  estimates  could  be 
divided  into  two  components:  (a)  additive  variation  due  to 
plant  vigor,  and  (b)  additive  variation  due  to  ’efficiency' 
of  nitrogen  fixation*  A  breeding  program  utilizing  these 
additive  components  was  discussed* 

Genotype  x  environment  interactions  were  of  significant 
importance  suggesting  that  study  concerning  quantitative 
variation  for  amount  of  nitrogen  fixed  should  be  conducted 
in  different  environments  with  various  strains* 
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INTRODUCTION 


The  development  of  tropical  agriculture  has  recently 
commanded  much  interest  in  an  effort  to  solve  pending 
problems  of  world  food  shortages.  The  improvement  of 
tropical  legumes  can  contribute  to  this  development  through 
the  improved  symbiotic  fixation  of  atmospheric  nitrogen.  It 
has  been  conceded  that  the  amounts  of  nitrogen  fixed  by 
tropical  legumes  were  generally  less  than  for  temperate 
legumes  (Vincent,  1974).  Since  temperate  and  tropical  legume 
species  differed  in  the  amount  of  nitrogen  fixed,  it  would 
be  pertinent  to  Inquire  as  to  whether  genotypes  within  a 
species  also  differed.  If  genotypes  differed  significantly, 
can  this  variation  within  the  host  for  symbiotic  nitrogen 
fixation  be  exploited  in  a  breeding  program  for  the 
Improvement  of  tropical  legume  species?  It  would  be  of 
particular  interest  to  initiate  a  genetical  study  of  the 
inheritance  of  symbiotic  nitrogen  fixation  in  the  tropical 
legume  Desmodium  because  of  its  great  variation  in  the 
morphological  differentiation  of  the  numerous  species 
(Bryan,  1966)  and  its  potential  as  a  forage  legume  (Bryan, 
1969  ). 

Norris  (1956)  postulated  that  the  legume-ff fol^obLuB 
symbiosis  originated  in  tropical  rain  forests.  As  the  legume 
adapted  to  the  subtropical  and  temperate  regions,  the 
symbiotic  interaction  appeared  to  become  more  specific. 
Since  the  Rhlzobl  urn  which  nodulates  £e£|&Q.gjLii2i  is  a  member  of 
the  "Cow-pea"  group,  a  promiscuous  interaction  group  of 
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£hJL2L2.b ium »  the  DesmodLum  -  Rhlzoblum  symbiotic  interaction 
may  represent  a  more  primitive  form  of  sytrblosis*  However, 
Bryan  ( 1966)  reported  that  this  Interaction  appeared  to  be 
specific  as  measured,  by  numerous  ineffective  nodulations* 

A*  Generalized  description  of  Desmodimn. 

Desmodiun)  is  a  self-fertile  tropical  legume  with  a 
somatic  chromosome  number  of  2n=22  (  Rotar  and  Urata,  1967  )* 
Rotar  and  Urata  (1967)  reported  collections  of  many 
different  species  from  tropical  countries  including  two 
species  from  Canada*  McWhirter  (  1969)  suggested  that  the 
center  of  origin  of  Desmodlum  was  in  South  America*  The 
taxonomic  identity  of  this  genus  is  uncertain  (McWhlrter, 
1969)  and  voucher  specimens  have  been  deposited  by  some 
Desmodlum  researchers  (Rotar  and  Urata,  1967;  McWhlrter, 
pers*  com*)  at  the  Harvard  University  Herbarium  for 
identifications  by  Dr*  Bernice  Shubert,  Arnold  Arboretum, 
Harvard  University* 

Four  Desmodlum  species  have  attracted  attention  as 
forage  legumes*  These  species  are  £•  sandwl cense  E*  Mey,  D* 
cftflum  (Gmel),  £>•  a  pa  ri  nes  (Link)  D*C**  and  £•  ilHC_LD.E-i.iiE. 
(Jacq)  D*  C • •  D*  qparines  and  D*  uncinatum  have  been  released 
as  the  commercial  cultivars  Green  leaf  and  Silverleaf 
desmodlum,  respectively,  in  Australia  ( Imrl e ,  1973)  as  they 

are  much  more  productive  and  robust  ( stem  length  is  5  to  20 


*Also  called  £•  in tor turn  (Mill)  Urb* 


•  •-  y*.  .  t 


_ 

-  -  - 

.i  -w.  -It  »•«  •'  «  "  v"' 

. 

•  £>W< 


- 

■ 


------  v ,  ~  • 


,«•■•  .w!Wte» 


3 


meters  )  than  either  D*  can  urn  or  D*  sandwi cense*  pfipupi  has 
attracted  attention  in  Hawaii  in  pasture  swards  (Younge,  g.t 
&I«  t  1964)*  X>*  sandw Icen se  is  the  least  robust  of  these  four 
species  with  the  stem  length  rarely  exceeding  two  meters* 
McWhirter  C1969)  reported  that  Ft  heterosis  for  forage  yield 
was  not  present  In  inter  or  Intra-specific  De smod j. up  crosses 
while  Rotar  and  Chow  (1971)  indicated  that  the  internode 
length  of  a  two-species  Fj  cross  was  intermediate  to  its  two 
parents*  These  data  suggest  that  Desytodj-upi  is  not  affected 
by  inbreeding  depression  for  forage  yield* 

Than  Aung  (  1970)  described  D*  gpndwj ppp ^e  as  a  postrate 
herb  whereas  Rotar  and  Chow  (1971)  observed  an  erect  growth 
habit*  The  stem  color  has  been  reported  by  Rotar  £t  g,l*  , 
(1967)  to  be  either  red— brown  or  green*  The  segregation  of 
this  character  is  3  red-brown  to  1  green  (Hutton  and  Gray, 
1967)*  Than  Aung  (  1970)  found  the  stems  to  be  green  with  red 
pigmentation  at  the  nodes*  The  stem  was  found  to  be  round 
(Rotar,  el  &l* ,  1967)  or  trisulcate  (Than  Aung,  1970)  in 
cross  section*  The  stems  are  pubescent  with  slightly  hispid 
hairs  (Than  Aung,  1970)  or  glandular  (Rotar  el  &.1  •  ,  1967)* 
Rotar  e_t  a^*  (  1967)  found  the  Internode  length  varied  from 
2*9  to  3*2  cm  which  was  shorter  than  £>•  uncinatum  or  EL* 
aoarl nes*  These  authors  believe  that  the  internode  length  is 
controlled  by  multiple  genes  since  the  interspecific  cross 
with  |>*  aparines  was  intermediate  to  the  parental 

phenotypes* 


The  leaves  of  Desmodluro  are  alternate. 
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trifoliate  with  the  terminal  leaflet  being  larger*  McWhirter 
(1969)  reported  a  dominant  polyphyllous  leaf  gene*  The  leaf 
form  of  D*  ssndfflcens^  is  lanceolate  with  acute  apices 
(Rotar  gjt  &l*  ,  1967)  or  ovate  and  cordate  at  the  base  (Than 
Aung)  1970)*  Rotar  ej:  al  *  (  1967)  found  the  length/width 
ratio  of  the  sixth  leaflet  from  the  flower  bud  to  range  from 
1*79  to  1*89  and  the  leaf  area  (length  x  width)  to  be  from  8 
to  12  sq  cm*  Than  Aung  (1970)  found  the  leaf  area  to  vary 
from  8  to  37  sq  cm  for  all  terminal  leaflets*  Hutton  and 
Gray  (  1967  )  reported  that  the  leaflets  of  gftpdwj  cense  and 
J2*  unci  natum  contain  a  silver-grey  marking  along  the  midrib* 
This  character  is  dominant  to  the  plain  leaf  and  segregated 
as  a  single  locus*  The  spotted  leaf  character  common  in  2* 
apafines  is  occasionally  found  in  gfindwjpftpge  (McWhirter, 
pers*  com*)*  Hutton  and  Gray  (1967)  reported  that  this 
character  was  dominant  to  the  plain  leaf  and  segregated  with 
a  3:1  ratio  in  the  F2  generation*  McWhirter  (pers*  com*  ), 
however,  suggested  that  the  spotted  leaf  character  is  a 
complementary  dominant  gene  system*  The  spotted  leaf  was 
found  to  be  closely  linked  with  the  stem  color  (Rotar  ej 
&1.  ,  1967). 

The  flowering  habit  of  D*  sandwicense  is  relatively 
indeterminant  and  differs  from  the  £>•  J±n£_i.na.!.!A“  and  D. 
aoarlnes  which  flower  in  response  to  short  days  (Rotar,  fit 
e,l .  *  1967)*  The  above  authors  also  report  that  ]>• 
sandwj cense  will  outcross  at  a  rate  of  18%*  McWhirter  (  pers* 


com* )  suggested  that 


is  highly  self  pollinated  In 
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■the  wild  and  estimates  the  outcrossing  to  be  about  20%*  He 
indicated  that  most  collections  give  an  impression  of  being 
rather  highly  inbred* 

The  flowers  of  Desmodium  are  borne  on  terminal  or 
axillary  racemes*  The  basal  flowers  of  the  raceme  are  first 
to  open  and  flowering  usually  continues  for  eight  days 
( Rotar  et  al*  *  1967)  with  about  five  flowers  opening  each 
day*  Preceding  the  opening  of  the  flowers?  the  raceme 
elongates  rapidly  separating  the  flowers*  The  petals 
elongate  beyond  the  sepals  the  day  before  the  flower  opens* 
The  anthers  dehise  at  this  time  and  the  pollen  is  mature* 
The  female  parts  of  the  flower  mature  the  following  morning 
when  the  standard  becomes  erect*  The  pistil  is  held  between 
the  keel  florets  and  when  tripped  it  halts  midway  to  the 
standard  without  striking  it*  The  pollen  is  discharged  in  a 
cloud  about  the  stigma*  Hutton  (1960)  reported  that  D* 
sandwicense  trips  spontaneously  while  Rotar  ejt  (1967) 
state  that  it  will  pollinate  without  tripping*  McWhirter 
( pers*  com*  )  reported  that  the  pistil  is  1  to  2  mm  longer 
than  the  stamens  and  that  pollination  will  not  occur  until 
tripped  or  disturbed*  Following  pollination,  the  flower 
wilts  rapidly  with  the  standard  folding  over  the  stigma. 

The  flower  color  in  Desmodium  is  polymorphic  varying 
from  white  to  purple  (Than  Aung,  1970)*  The  purple  flower 
color  is  dominant  to  white  and  this  character  shows  a 
complementary  dominant  gene  inheritance  in  crosses  of 
certain  white  flowered  species  (  D. 
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aparines )  (McWhirter,  pers.  com* )• 

The  seed  pods  of  Desmodlum  are  distinctly  articulated 
and  contain  an  average  of  6*3  seeds  per  pod  (  Rotar  £t.  a.1  •  ♦ 
1967)*  The  pods  are  covered  with  short  hairs  and  stick 
readily  to  clothing  or  fur*  The  seeds  are  kidney-shaped  and 
light  brown  in  color*  The  1000-kernel  weight  of  &• 
sandwi cense  was  3*5  gm  falling  between  I>*  uncinatum  (4*0  gm ) 
and  JO*  apftrlnes  (1*8  gm  )  (Rotar  and  Urata*  1966  )* 

The  germination  for  machine  harvested  seed  is  99%  but 
germination  is  reduced  when  hand  harvested*  Hand  harvested 
seeds  may  be  scarified  by  chipping  the  testa  (McWhirter, 
pers*  com*)  or  immersion  In  concentrated  sulfuric  acid  for 
15  min  (  B*  Groff,  pers*  com*  )•  Following  scarification, 
germination  is  normal* 

McWhirter  (1963)  first  reported  successful  crosses 
between  JD*  sandwlcens  e  and  JD*  ftPftrlngg*  Hutton  and  Gray 
(1967)  reported  successful  crosses  between  £•  sapdwi cgpge 
and  D.  unc  i  nat um  and  D.  uncinatum  and  £>•  aoarines*  They 
found  no  gross  cytological  barriers;  however, 

incompatab i li t ies  were  found  by  McWhirter  (pers*  com*)  when 
specific  lines  were  used*  McWhirter  (pers*  com*)  suggested 
that  JD*  aoarlnes  be  used  as  a  bridging  species* 

Rotar  and  Chow  (1971)  found  that  pod  formation  was 
highest  when  £•  sandwlc  ense  was  the  female  parent  but  much 
reduced  when  it  was  the  male  parent*  McWhirter  (pers*  com*) 
also  found  a  reduced  pod  set  when  fi*  S&U^ll££D.££  was  the 
male  parent*  Rotar  and  Chow  (1971)  found  crosses  between  D. 
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unclnatum  and  £. 


to  be  low  for  pod  set  and  they 


concluded  that  physiological  factors  may  be  influential* 
Pollen  abortion  was  highest  in  the  2*  sa pdwj pepge  x  D* 


crosses  and  was  lowest  in  the  D*  sandwi cense  x  D* 


uncina  turn  crosses  (Rotar  and  Chow,  1971).  Rotar  and  Chow 


(  1971  )  concluded  that  2* 


was  more  closely 


related  to  D«  uncj.  nft  tym  than  to  2*  gtpfirlP?S«  They  reported 
that  pod  set  was  negatively  correlated  with  percentage 
pollen  abortion  (r  =  *74,  df  =  21)*  McWhirter  (  pers  •  com*) 
reported  that  morphological  differentiation  was  correlated 
with  genetic  divergence*  He  concluded  that  D*  unc j na turn  was 
morphologically  distinct  from  D*  ^gindwi  cen  se  and  D* 
asanines . 

Chow  and  Crowder  (  1972)  report  a  cross  between  D*  £&Q.y.ta 
and  D*  unc  ina turn*  They  conclude  that  2*  c  g_n  urn  is  genetically 
separate  from  2*  sandwi cense «  2*  uncina turn  and  2*  apftrlnes  • 

B*  Genetic  studies  of  symbiotic  nitrogen  fixation* 

The  study  of  symbiotic  nitrogen  fixation  has  been  of 
considerable  interest  to  agriculture  for  numerous  years  and 
this  area  has  recently  been  comprehensively  reviewed 
(Qulspel,  1974;  Hardy  and  Gibson , 1 977 ) •  The  limitation  of 
genetlcal  studies  of  the  inheritance  of  symbiotic  nitrogen 
fixation  was  of  major  concern  (Dart,  1974)*  The  majority  of 
the  genetical  studies  Involved  the  study  of  major  genes 
conditioning  non-nodu la t ing  or  Ineffective  noduletJon 
interactions*  Few  genetlcal  studies  have  been  reported 
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giving  results  for  plant  variation  for  total  nitrogen  fixed* 

Major  genes  conditioning  strain-specific  nodulation 
have  been  reported  in  alfalfa  (  Medi pfigo  Sfttivfi  L*)  (Burton 
and  Wilson,  1939;  Gibson,  1962),  red  clover  (  Tr  1  f  pi  i  ynn 
pra tense  L*  )  (Nutman,  1954),  subterranian  clover  (  X* 
sub t erraneum  L*  )  (Gibson,  1964),  soybean  (  Glycine  mgiy  L* 
Merril  )  (Williams  and  Lynch,  1954;  Caldwell,  1966;  Vest, 
1970;  Vest  and  Caldwell,  1972),  end  peas  (  Pi  sum  g.a.iiyum  ) 
(Holl,  1975;  Lie  gjh  a,  l*  ,  1976  )  *  The  major  genes  appeared 
to  be  of  two  types:  (a)  recessive  genes  with  modifiers,  or 
(b)  dominant  genes*  The  investigations  in  Trlf ol i um 
suggested  that  s tra in— spec i f ic  ineffective  nodulation  was 
controlled  by  a  single  recessive  gene  modified  by  recessive 
suppressor  genes  which  restored  effectiveness  (Nutman,  1954; 
Gibson,  1964)*  Single  dominant  genes  (RJ2»  RJ3  and  RJ*  )  were 
reported  in  soybeans  conditioning  strain-specific 

ineffective  nodulation  (Caldwell,  1966;  Vest,  1970;  Vest  and 
Caldwell,  1972)  while  the  recessive  gene  (  r Ji  )*  conditioned 
a  universal  non— nodulat ion  interaction  in  soybeans  (Williams 
and  Lynch,  1954)* 

Vest  aJL*  (  1973)  classified  the  soybean  plants 

homozygous  for  the  RJ  alleles  (cultivers  -  Hardee,  Hill)  as 
nodulated-ineffective  since  the  plants  exhibited  small, 
white  cortical  proliferations  at  the  node  sites*  They 
concluded  that  infection  had  taken  place.  However,  nitrogen 


^Designated  no  by  Williams  and  Lynch  (1954). 
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was  not  being  fixed.  Vest  el  &l.  (1973)  reported  that  Rj2 
and  R  j 3  appeared  to  be  linked  yet  both  were  independent  of 
rji*  Although  RJ2  and  Rj3  were  similar  phen© typically , 
producing  numerous  small,  white  cortical  proliferations, 
they  appeared  to  differ  in  the  serological  properties  of  the 
ineffective  strains.  Rj2  was  ineffective  with  most  strains 
of  the  cl  and  122  serogroups  while  RJ3  was  ineffective  with 
only  one  particular  strain.  Caldwell  et  &.!•  (  1966)  performed 
reciprocal  grafts  of  the  cultivars  Hardee  (  R J2 ,  Rj2  )  and  Lee 
(rj2,  rj2 )  and  reported  that  the  root  genotype  dictated  the 
nodulation  interaction. 

The  ineffective  response  for  the  RJ4  gene  (Hill  x 
strain  61)  appeared  to  be  a  specific  response  in  the  sense 
of  the  RJ3  gene.  The  RJ*  gene  produced  only  a  few  small, 
white  cortical  proliferations  since  nodule  abortion  probably 
occurred  at  a  very  early  stage  of  development  (Vest  and 
Caldwell,  1972). 

A  universal  non— nodula t i on  gene  has  also  been  reported 
for  peas.  Holl  (1975)  and  Holl  and  LaRue  (1975)  studying 
crosses  of  an  Afghanistan  cultivar  (  non— nodulat ing )  reported 
two  genes  affecting  nodulation  and  fixation.  Holl  (1975) 
reported  that  the  Sym2  locus  (nodulation)  and  the  Sym3  locus 
(fixation)  interacted  to  give  four  phenotypes  for  fixation. 
Plants  with  the  genotype  Sytn2  — ,Syn»3  —  were  classed  as 
nodulat i ng  and  fixing  while  plants  with  the  genotype  sym2 
sym2,  —  —  were  non-nodulat ing  and  non— fixing.  However, 


when  the  Sym2  gene  was  combined  with  the 


homozygous 
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recessive  condition  for  the  syni3  gene  the  resulting 
combination  conditioned  differing  interactions.  If  the  Sym2 
gene  was  heterozygous  the  resulting  interaction  was 
nodulation  but  no  fixation*  whereas*  a  homozygous  dominant 
condition  for  the  Sym2  gene  resulted  in  nodulation  and  low 
fixation.  The  Sym2  locus  conditioned  nodulation  and 
segregated  3  nodulatlng  to  1  non— nodula ti ng  phenotypes. 

Lie  e_£  &l.  (  1976)  reported  a  strain-specific* 
tempe ra ture— sensi t 1 ve  non— nodulat ion  interaction  in  the  pea 
cultivar,  Iran.  The  authors  indicated  that  nodulation 
occurred  only  when  the  root  was  subjected  to  favorable 
nodulation  temperatures  (  26  C  )  for  the  second  and  third  day 
following  inoculation.  Adverse  temperatures  during  this 
period  reduced  nodulation  whereas  adverse  temperatures 
before  or  after  the  critical  period  had  no  effect.  Lie  £t 
al.  (1976)  reported  that  this  strain-specific*  temperature- 
sensitive  ineffective  nodulation  interaction  was  conditioned 
by  a  dominant  allele  at  a  single  locus. 

The  study  of  symbiotic  nitrogen  fixation  has  generally 
dealt  with  qualitative  interactions  relating  to  nodulation 
vs.  non— nodulati on •  However,  genetical  studies  have  been 
reported  for  quantitative  variation  for  estimates  of 
nitrogen  fixation.  Burton  and  Wilson  (1939)  reported  that 
varieties  of  alfalfa  differed  significantly  for  total 
nitrogen  fixed  when  inoculated  by  several  strains  of  R. 
meli loti .  Conversely  the  authors  also  noted  significant 
differences  when  strains  were  compared  as  well  as  specific 
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host— strain  interactions*  These  data  suggested  quantitative 
differences  among  hosts  and  among  strains*  Burton  and  Wilson 
(1939)  also  noted  that  the  results  tended  to  be  erratic  and 
attributed  these  deviations  to  environmental  factors*  Later 
work  in  alfalfa  by  Gibson  (1962)  confirmed  genetic  variation 
of  the  host— strain  interaction*  He  noted  that  the  Canadian 
cultivar  Rambler  was  inefficient  when  inoculated  with 
Australian  strains  of  f£*  me l iloti  yet  was  efficient  when 
inoculated  with  Canadian  strains*  Gibson  (  1962  )  suggested 
that  the  M*  falcata  parentage  of  Rambler  may  be,  in  part, 
responsible  for  this  inefficiency*  Recently  Seetin  and 
Barnes  (1977)  reported  genetic  variation  in  alfalfa  for 
quantitative  nitrogen  fixation  estimated  by  the  acetylene 
reduction  assay  and  dry  weight  accumulation*  These  authors 
indicated  that  the  alfalfa  clones  expressed  Fi  heterosis  and 
inbreeding  effects  for  the  estimates  of  nitrogen  fixation* 
Duhiga  et  &1*  (1978)  reported  that  F|  progeny  from  Mesilla 
alfalfa  plants  selected  for  high  acetylene  reduction  values 
expressed  significantly  greater  acetylene  reduction  values 
than  the  parental  population*  Significant  results  were  not 
reported  for  the  low  selection  line* 

In  recent  research  with  field  beans  (  Vj ci a  f aba  ),  El  — 
Sherbeeny  et  si*  (  1977a, b)  reported  genetic  variation  among 
Rhizoblum  leaumi nos a rum  strains  and  among  field  bean  lines, 
respectively*  Mytton  et  a l *  ( 1977)  indicated  that  the 
genetic  variation  for  lines  or  strains  was  found  to  be  minor 


in  comparison  to  the  line  by  strain  interactions*  Mytton  et 
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al«  (1977)  partitioned  the  genetic  variance  into  three 
components;  namely,  (a)  the  additive  genetic  effect  of  the 
host  genotype,  (b)  the  additive  genetic  effect  of  the  strain 
genotype,  and  (c)  the  non-additive  effect  of  the  host  by 
strain  Interactions.  They  Indicated  that  73%  of  the 
phenotypic  variation  could  be  directly  attributed  to  host 
strain  interactions  and  concluded  that  simultaneous 
selection  for  a  specific  host— strain  interaction  would  be 
most  effective* 

Mytton  (  1S75),  studying  white  clover  (  Tpi  f  olj  upn  ) 
symbiosis,  found  that  host  variation  for  nitrogen  fixation 
was  larger  than  for  differences  among  strains.  He  concluded 
that  only  effective  strains  had  been  collected  thereby 
reducing  genetic  variation*  However,  Mytton  (1975)  did 
report  significant  host  x  strain  interactions*  Strains  which 
were  selected  from  one  plant  generally  entered  into  a  more 
productive  interaction  with  relatives  of  that  plant  than 
with  non-related  plants  suggesting  selective  hosts. 

Minchin  al*  (1978)  also  found  genetic  variation 
among  cultlvars  of  cowpea  (  Vlana  UESEiO-liiaLt g,  )-U  )  feat.  c.&± 
among  Rhizobiutn  strains.  The  above  authors,  unlike  Mytton  et 
fi.1.  (  1977  ),  concluded  that  any  breeding  program  should 
concentrate  on  the  macrosymbiont. 

Quantitative  assessments  of  the  amounts  of  nitrogen 
fixed  had  been  measured  in  earlier  experiments  by  total  dry 
weight  accumulation  of  the  shoots  when  plants  were  grown  in 
nitrogen  free  containers  (Gibson,  1962;  Vincent,  1970).  The 


‘  1  ’ 
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1SN  method  (Burris  and  Wilson,  1957)  and,  recently,  the 
acetylene  reduction  assay  (Hardy  g.£  a.1*  ,  1968)  have  been 
developed  as  assays  of  the  activity  of  the  enzyme 
nltrogenase  (Figure  1  )•  Of  these  two  assays,  the  acetylene 
reduction  assay  has  been  most  popular  because  it  was  simple, 
rapi d  and  more  sensitive  than  the  1SN  method;  but  it  was  an 
indirect  approximation  (Burris,  1974)*  Since  its 
introduction,  the  acetylene  reduction  assay  has  appeared  in 
many  publications  measuring  nltrogenase  activity  in  Intact 
nodules,  nodule  prepared  cell— free  extracts,  field 
experiments,  as  well  as  marine  and  fresh-water  habitats 
(Quispel,  1974)* 

The  reduction  of  acetylene  to  ethylene  by  nltrogenase 
was  independently  discovered  by  Schollhorn  and  Burris  (1966) 
(reviewed  by  Burris,  1974)  and  Dilworth  (1966)*  Its 
usefulness  as  an  Index  of  nltrogenase  activity  was 
immediately  realized  by  Hardy  e t  §*!•  (  1968  )  who  developed  a 
method  of  measuring  nltrogenase  activity  for  field 
experiments*  The  relationship  between  acetylene  and  Ng 
reduct  ion  may  be  represented  ast 

HC  -  CH  +  2e~  +  2H+  =  H2C  -  CH2 

N  -  N  +  6e~  +  6H+  =  2NH3 

These  formulae  suggested  that  the  conversion  rate  for  N2  to 
NH3  was  one-third  that  of  C2H2  to  C2H*  since  acetylene 
reduction  required  two  electrons  while  N2  reduction  required 
six.  Burris  (1974)  reported  that  the  conversion  factor  of 
1/3  was  not  generally  found  but  was  much  closer  to  1/4  and 
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Figure  1.  A  diagramatic  representation  of  the  symbiotic  nitrogen  fixation 
pathway  in  the  host  nodule  cell. 
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Burris  suggested  a.  comparison  of  the  acetylene  reduction 
assay  with  the  1 5N  method  to  determine  an  exact  conversion 
factor  when  required  for  a  particular  experiment* 

Recent  studies  with  the  acetylene  reduction  assay  in 
soybeans  suggested  that  nodule  activity  (C2H4)  was 
determined  by  the  physiological  stage  of  the  plant  (Lawn  £t 
al* »  1974)*  Lawn  and  Brun  (1974)  found  that  the  nodule  fresh 
weight  and  acetylene  reduction  activity  for  soybeans  reached 
a  maximum  near  the  termination  of  flowering  and  declined 
rapidly  during  the  onset  of  the  pod-filling  stage*  The 
authors  attributed  this  decline  to  a  competition  by  the  pods 
for  available  assimilate*  Additional  studies  with  reciprocal 
grafts  in  soybeans  Indicated  there  were  significant 
differences  among  shoot  and  root  genotypes  as  well  as  for 
sampling  dates  as  measured  by  the  acetylene  reduction  assay 
(Lawn  and  Brun*  1974)*  Significant  shoot  genotype  x  sampling 
date  interactions  were  found  for  all  characters  studied 
while  root  genotype  x  sampling  date  interactions  were  not 
significant  with  the  exception  of  the  photosynthetic  rate* 
There  were  significant  shoot  genotype  x  root  genotype 
interactions  for  nodule  number  and  nodule  size*  Nltrogenase 
activity  (  C2H4)  for  the  shoot  genotype  means  was 
significantly  correlated  with  nodule  fresh  weight*  root  dry 
weight*  shoot:root  ratio*  photosynthesis  and  shoot  dry 
weight*  However  nitrogenase  activity  (C2H4)  for  root 
genotype  means  was  not  significantly  correlated  with  any  of 


the  above  characters* 


Lawn  gjb  &JL*  (  1974)  concluded  that  root 
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genotypes  may  differ  in  their  ability  to  efficiently  utilize 
the  phot osyn tha te  assimilated  during  symbiotic  nitrogen 
fixation*  They  also  suggested  that  the  root  genotype 
controlled  nodule  initiation  and  nodule  size  which  were 
found  to  be  significantly  negatively  correlated  with  the 
root  genotype  effect  on  nodule  number* 

Dart  (  1974)  reported  responses  to  selection  in  several 
Tr 1  foil um  species  for  number  of  nodules  per  plant,  time 
taken  to  form  the  first  nodule  and  effectiveness  of  the 
nodules*  These  characteristics  were  also  affected  by  the 
Rh izobium  strain  used*  Dart  (1974)  concluded  that  plant 
breeding  for  increased  nodulation  was  of  significant 
potential* 


C*  Objectives* 

It  was  the  objective  of  this  study  of  the  inhe 
of  symbiotic  nitrogen  fixation  in  Degcnodiutn  to  es 
whether  or  not  there  exists  genetic  variation 
processes  of  the  symbiotic  relationship  between  hos 
and  Rh Izobium  with  particular  emphasis  upon  effect 
ineffective  nodulation  and  the  quantitative  varia 
total  symbiotic  nitrogen  fixed  per  plant*  It  was  a 
interest  to  determine  the  genetic  basis  fo 
differences,  if  disclosed*  The  dry  weight  accumulatio 
bottle  jar  (Leonard  jar  assembly)  and  the  ac 


reduction  assay  were  assessed  to  determine  their  po 


r i tance 
tabl i sh 
in  t  he 
t  plant 
i  ve  vs 
tion  in 
Iso  of 
r  such 
n  in  a 
etylene 
ten ti al 


es  plant  breeding  tools* 


.  ,  X 
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MATERIALS  AND  METHODS 


A*  Outline# 

In  this  study  of  the  inheritance  of  symbiotic  nitrogen 
fixation  in  Deg,  modi  Utn  #  experimentation  was  conducted  at 
three  Levels#  These  levels^  in  increasing  order  of  interest, 
have  been  broadly  defined  as  methodology,  qualitative 
genetic  variation,  and  quantitative  genetic  variation  (Table 
1)#  The  methodological  studies  encompassed  the  background 


Table  1#  Levels  of  study  of  the  inheritance  of 

symbiotic  nitrogen  fixation  in  the  host# 


Level  1 

Methodology 

Comments  Z 

Preliminary  studies 
How?  What?  Where? 
Improvements?  When? 

Level  2 

Qualitative  Genetic 

Variat ion 

Host-strain  interactions* 

effective  ineffective 

Which? 

( Major  genes? ) 

Level  3 

Quantitive  genetic 

varl at i on 

Fixa  ti on 

ef f icient 

M  efficien cy" * 

inef f 1c  lent 

How  much? 

( Minor  genes? ) 

Ul t ima te 

goal 

Why? 

Select! on 

Super-fixing  inbred 
lines  or  hybrids 

*May  also  be  modified  by  strain  manipulation# 


studies  necessary  for  the  development  and  improvement  of  the 


technique  of  studying  symbiotic  nitrogen  fixation# 
Experiments  in  qualitative  genetic  variation  were  conducted 
to  determine  the  genetic  basis  for  ineffective/effective 
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nodu tat i  on  interactions*  The  experiments  involved  “the 
testing  of  the  observed  segregation  ratios  in  F2  and 
backcross  populations  for  a  goodness  of  fit  to  an  expected 
segregation  ratio* 

The  experiments  involving  quantitative  genetic 
variation  constituted  the  major  interest  in  this  study  and 
relate  to  the  Investigation  of  the  inheritance  of  total 
symbiotic  nitrogen  fixed  per  plant*  The  experimental  designs 
which  were  chosen  to  establish  and  examine  the  genetic 
variation  were:  (1)  a  five— parent  diallel  cross*  and  (2)  a 
cross  between  two  Inbred  lines*  The  diallel  cross  consisted 
of  inbred  lines  of  J)*  sandwi cense  which  were  not  selected 
for  the  characters  under  study*  The  cross  between  two  inbred 
lines  constituted  crosses  of  D*  san dyj pgpge  with  two  other 
species  of  Pesmodj  um  which  were  conspicuously  different  from 
£•  sandwlcense  as  well  as  a  within  £>•  sandwicense  cross* 


B*  Materials* 


1 «  Rhlzoblum  s t r a ins* 

To  study  one  com 
of  symbiotic  nitrogen  f 
other  component  constan 
population  created  from 
major  interest*  the 
experiment  has  been  hel 


ponent  of  the  host-strain  interaction 
ixatlon*  it  was  necessary  to  hold  the 
t*  Since  the  variation  within  the 
crosses  among  Desmodium  lines  was  of 
Rh  1  zobi  upi  component  within  each 
d  constant  and  was  considered  to  be 

pure  strains  of 


an 


environmental 


component 


The  two 
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Ekiz.ob.Lum  selected  were  strain  A  (CB627)*  and  strain  B 
(CB1789)#  which  were  isolated  from  D*  a pg§  and  D  » 
e ,  respectively.  These  strains  differed  in  growth 
properties.  Both  strains  are  classed  as  slow  growing; 
however,  strain  A  is  considered  to  be  relatively  faster 
growing  than  strain  B,  Both  strains  are  members  of  the  "Cow- 
pea"  group  of  Rhjjzo.klutq  (Vincent,  pers,  com,)  and  strain  A 
has  been  released  in  Australia  as  a  commercial  inoculum  for 
Desmodlum , 

2 •  Desmod i urn  lines. 

The  Desmod i um  lines  used  in  this  study  are  described  in 
Table  2,  J)«  san dw i c e n se  or  its  crosses  were  used  in  most 

experiments  because  it;  (a)  is  diploid,  self-fertile  and 
readily  crossed,  (b)  is  relatively  insensitive  to 
phot operi odism  and  is  easily  maintained  under  glasshouse 
conditions,  (c)  has  a  robust  seedling  (desirable  for  culture 
in  a  bottle  jar  assembly),  and  (d)  can  be  hybridized  easily 
with  other  conspicuously  variable  and  commercially  important 
species  of  Desmodlum  (Than  Aung,  1970),  Prolonged  inbreeding 
does  not  adversely  effect  plant  vigor  and  Inbreeding 
depression  is  generally  considered  to  be  absent  (McWhirter, 
1969;  Rotar  and  Chow,  1971),  These  characteristics  made  D, 
sandwl cense  an  excellent  model  organism  for  nitrogen 
fixation  studies. 


^Commonwealth  of  Australia  identification  number. 


■ 


Table  2#  Description  of  the  Pesmofl jum  lines  involved  in  this  study 
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Five  inbred  lines  of  £•  sandwlcepge  ( Dsl  to  Ds5 )  were 
crossed  in  all  possible  combinations  including  reciprocal 
crosses  to  create  a  dlallel  cross.  Although  the  inbred  lines 
were  not  selected  for  the  characters  studied,  the  lines  have 
been  considered  as  a  fixed  sample  of  D •  sandwj  ccpge  lines. 
The  six  populations  (Ptf  P2*  Fi ,  F2,  BC!  ,  and  BC2  )  derived 
from  a  cross  between  two  Inbred  lines  were  made  for  each  of 
the  following  crosses:  Dsl  x  Dal,  Ds2  x  Dsp,  and  Ds2  x  Ds5. 
The  diallel  cross  and  the  three  crosses  between  two  inbred 
lines  constituted  the  basic  material  to  determine  whether 
genetic  variation  among  differing  populations  existed  for 
symbiotic  nitrogen  fixation  in  Despaodj ytm, 

C,  Experimental  Procedures, 

1,  Crossing  technique. 

Two  crossing  methods  have  been  used  in  Desmodi urn . ,  the 
emasculation  method  (Rotar  ejt  al,  ,  1967)  and  the  tripping 

method  (McWhirter,  pers,  com,  ),  The  emasculation  method  is 
achieved  by  opening  the  flower  bud  the  night  before  maturity 
and  gently  removing  the  anthers.  The  racemes  with  the 
emasculated  flowers  are  then  covered  with  a  glassine  bag 
containing  a  moist  ball  of  cotton.  During  the  morning  hours, 
the  flowers  are  pollinated  with  the  desirable  pollen  using  a 
sandpaper— tipped  toothpick. 

In  the  tripping  method  the  mature  flower  is  gently 
tripped  and  the  stigma  is  examined  with  a  10X  hand  lens. 


Contaminated  flowers  are  discarded. 


The  remaining  flowers 
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ar*e  pollinated  and  following  pod  formation  glessine  bags  are 
placed  over  the  raceme  to  avoid  seed  contamination* 

Rotar  et  fil*  (  1967  )  found  the  emasculation  technique  to 
be  superior  but  more  difficult  than  the  tripping  method* 
They  reported  that  the  emasculation-wet  cotton  method 
produced  60%  crossed  seed  from  a  10%  flower  set  compared  to 
the  tripping  method  which  yielded  40%  crossed  seed  from  a 
13%  flower  set*  The  authors  only  brushed  the  stigma  in  the 
tripping  method  which  resulted  in  the  low  percentage  crossed 
seed*  Rotar  et  g/l*  (  1967)  and  Hutton  (  1960)  both  noted  that 
crossing  was  most  successful  when  the  temperature  was  cool 
and  the  relative  humidity  was  high*  Rotar  and  Chow  (1971) 
report  that  percentage  pod  set  and  the  monthly  average 
temperatures  were  negatively  correlated  (r  =  -*99*  df  =  3)* 
The  crossing  method  used  in  this  study  was  a 
modification  of  the  tripping  method*  Each  flower  was  grasped 
by  the  thumb  and  forefinger  at  the  calyx  and  tripped  by 
applying  pressure  at  the  base  of  the  standard*  Following 
trlpping»  the  anthers  and  pollen  were  brushed  away* 
Beginning  with  the  top  flowers*  the  stigmas  of  the  tripped 
flowers  were  thoroughly  washed  with  a  2*5%  (by  weight)  sugar 
solution  and  Immediately  pollinated  with  the  desired  pollen* 
The  pollen  was  either  taken  directly  from  the  pollinating 
flower  which  was  prepared  by  holding  the  flower  with  the 
thumb  and  forefinger  at  the  calyx  and  gertly  removing  the 
petals  without  disturbing  the  pollen  or  by  tripping  the 
pollinating  flower  directly  onto  a  fingertip*  The  finger  w as 
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dipped  into  alcohol  between  pollinations.  Following  pod 
formation  the  raceme  was  covered  with  a  glassine  bag  to 
avoid  seed  contamination.  Crosses  were  made  at  the  Colombian 
Institute  of  Agriculture,  Medellin,  Colombia  and  at  the 
University  of  Sydney,  Sydney,  Australia.  Study  of  the 
progeny  from  some  of  these  crosses  showed  that  95%  of  the  F| 
crosses  were  hybrids  for  the  spotted  leaf  character 
( Appendix  A  )  • 

2.  Bottle  jars. 

The  bottle  jars  (Leonard  jar  assembly),  for  nitrogen- 
free  growth  containers,  were  constructed  from  various 
inexpensive  and  readily  available  materials.  The  bottle  jar 
which  conformed  to  the  description  of  Vincent  (1970)  is 
illustrated  in  Figure  2.  During  this  study,  the  bottle  jar 
was  modified  to  facilitate  the  removal  of  the  plant  roots 
intact  for  the  acetylene  reduction  analysis.  The  modified 
bottle  jar  is  illustrated  in  Figure  3.  After  filling  the 
spirits  bottle  with  sand  and  compacting  the  sand  about  the 
wick,  the  sand  was  wetted  and  the  reservoir  was  filled  with 
a  nutrient  solution  composed  of  the  following: 


CaHPO*  1.0  g 

k2hpo4  .2  g 

MgS04.7H20  .2  g 

NaCl  .2  g 

FeC I3  .  1  g 

Tap  water  5.0  l 
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Glass  Cover 
(removed  after 
seedling  emergence) 

Dry  Gravel  Mulch 
(placed  after 
seedling  emergence) 

Water-Proof  Tape 


740  ml  Spirits  Bottle 

Well-Washed  River  Sand 
(4  mm  maximum) 

Cotton  Wick 


Water-Proof  Paper  Bag 


Nutrient  Solution 


Reservoir  Large  -  2  I 

Small  —  600  ml 


Cotton  Ball 


Figure  2.  Initial  bottle  jar  design  (Vincent,  1 970) 
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Glass  Cover 
(removed  after 
seedling  emergence) 

Dry  Rock  Mulch 
(placed  after 
seedling  emergence) 

Water-Proof  Tape 


740  ml  Spirits  Bottle 


Well-Washed  River  Sand 
(1.5  mm  Maximum) 


Water-Proof  Paper  Bag 


Nutrient  Solution 


Cotton  Wick 


Reservoir  Large  —  2  I 

Small  —  600  ml 

Cotton  Ball 


Figure  3.  Modified  bottle  jar  design. 
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bottle  jfi  r  was  covered  with  a  waiter*”  pro  of  paper  bag 
which  was  secured  with  water-proof  tape,  capped  with  half  of 
a  petri-dish,  and  autoclaved  for  2  hours  at  15  p.s.i..  After 
cooling,  the  bottle  jars  were  planted  with  surface  sterlized 
seeds.  The  seeds  were  sterilized  in  an  acidified  0.2%  HgCl3 
solution  (Vincent,  1970)  and  washed  with  six  changes  of 
sterile  water.  The  testa  of  the  seeds  had  been  previously 
chipped  to  facilitate  germination.  Following  germination  and 
inoculation  the  bottle  jars  were  covered  with  a  pebble  mulch 
(Vincent,  1970)  and  left  undisturbed  until  harvest. 

3.  Rhizobi um  culture. 

Stock  cultures  of  the  Rhizoblpm  strains  were  maintained 
as  agar  slopes  and  stored  at  2  C  in  McCartney  bottles. 
Inoculum  for  the  experiments  was  grown  as  broth  cultures  in 
350  ml  flats.  The  agar  slopes  and  broth  cultures  were  made 
from  yeast-mannitol  medium  (Vincent,  1970)  which  consisted 
of  the  following  ingredients: 


k2hpo4 

MgSO*. 7H20 
NaCl 
CaCC>3 
Manni tol 
Yeast  water 
Distilled  water 
Agar  (when  required) 
Fifty  ml  of  the 


0.5  g 

0.2  g 
0.  1  g 

3  m 

10  g 

1 00  ml 

900  ml 
15  g 

above  medium  was  placed 


i  n 


each  flat 


and  autoclaved  at  15  p.s.i.  (  120  C)  for  15  min.  The  sterile 
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flats  were  stored  at  2  C  until  required#  Yeast  water  was 
freshly  prepared  by  mixing  100  g  of  fresh  compressed  baker's 
yeast  in  one  liter  of  cool  water#  I  he  mixture  was  allowed  to 
stand  for  1  to  2  hours  at  room  temperature  before 
autoclaving  for  50  min*  The  suspended  material  was  allowed 
to  settle  and  the  clear  supernatant  was  used  as  the  yeast 
water* 


Depending  on  the  strain  of  Rfclzokium,  either  ten  to 
fourteen  days  prior  to  the  inoculation  time  of  the  plants,  a 
flat  was  removed  from  storage  and  inoculated  with  the 
desired  strain  using  observed  bacteriological  precautions* 
The  inoculated  flat  was  placed  in  a  growth  cabinet  at  26  C 


for  the  required 

time 

to  facilitate 

the 

growth  of 

Ejllz.obla.rn «  For 

the 

inoculation  of 

the 

experiments 

suspension  of  the  Rh i zoblum  was  made  by  filling  the  flat 
with  sterile  water  and  shaking  vigorously  to  disperse  the 
Rhi zoblum.  Five  ml  of  this  suspension  was  applied  to  the 
base  of  each  plant  3  to  4  days  following  emergence,  which 
was  approximately  1  week  after  planting* 


4*  Quantitative  variables* 

(a)  Acetylene  reduction  assay* 

The  acetylene  reduction  assay  of  Hardy  gt.  al*  (  1968  ) 
was  used  to  assess  the  ni trogenase  activity  of  the  nodulated 
roots*  Individual  plants  were  harvested  from  the  bottle  jars 
and  the  root  section  from  each  plant  was  placed  in  a  13*5  ml 
vial  which  was  plugged  with  an  air  tight  seal*  The  vials 
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were  immediately  placed  on  a.  vacuum  system,  refilled  with 
argon  and  Injected  with  2.8  cc  of  oxygen  and  3*5  cc  of 
acetylene  (ethylene  free).  After  60  min,  a  200  microliter 
sample  was  removed  from  each  vial  and  injected  into  a  F-ll 
Perkin— Elmer  gas  chromatograph  with  a  hydrogen  flame 
ionization  detector.  The  detector  was  attached  to  a  Tohshin 
electron  chart  recorder.  Nitrogen  was  used  as  the  carrier 
gas  and  the  column  was  a  stainless  steel  tube  of  30  cm  in 
length  and  .32  cm  in  diameter  and  filled  with  Porapack  T 
packing  material.  The  column  temperature  was  held  at  60  C 
giving  a  retention  time  of  15  sec  for  ethylene  and  30  sec 
for  acetylene.  At  periodic  intervals,  duplicate  samples  were 
processed  to  check  the  accuracy  of  sampling.  An  acetylene 
reduction  estimate  (  nanomolar( nM  )/pl/hr  )  was  calculated  from 
peak  height  using  the  following  equation: 

Es(  Ef  -  Eo) 

Acetylene  reduced  =  - X  Ac  X  St 

As(  Af  —  Ao  ) 

where  Ef  and  Af  were  the  final  ethylene  and  acetylene  peak 
heights,  respectively,  Eo  and  Ao  were  the  origin  values  for 
ethylene  and  acetylene,  and  Es  and  As  were  the  scales  of 
measurement  for  ethylene  and  acetylene.  Ac  was  the  acetylene 
constant  which  was  determined  as  the  mean  of  the  acetylene 
peak  heights.  St  was  the  correction  for  a  standard  known 
sample  of  ethylene. 

Acetylene  reduction  assay  conducted  in  Experiment  VI 
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generally  followed  the  procedure  outlined  above,  except  that 
a  35  ml  vial  was  used  with  no  argon*  The  ethylene  peak  was 
used  as  an  estimate  of  acetylene  reduced* 

(b)  Dry  weight  accumulation* 

Dry  weight  accumulation  data  were  collected  for 
individual  plants  grown  in  bottle  jar  assemblies.  The 
harvested  plants  were  clipped  at  the  cotyledonary  node  and 
the  shoot  portion  was  retained  as  the  dry  weight 
accumulation  sample*  The  leaves  dropped  during  the  course  of 
the  experiment  were  collected  and  included  in  the  sample* 
The  samples  were  oven-dried  at  105  C  for  24  hours  and  left 
at  room  temperature  for  at  least  one  week  before  weighing  to 
the  nearest  milligram* 

(c)  Shoot  length* 

The  shoot  length  was  measured  either  j_n  g.i_tu  or  on  the 
harvested  plants  before  oven  drying  with  measurements  being 
recorded  to  the  nearest  millimeter  for  the  primary  shoot 
length  and  the  two  secondary  shoots  originating  at  the 
cotyledonary  node*  Total  shoot  length  was  calculated  by 
summing  the  primary  and  the  two  secondary  shoot  lengths* 

D*  Experimental  designs  and  analyses* 

1*  Methodological  studies* 

(a)  Selection  for  leaf  size* 


Prior  to  the  commencement  of  the  nitrogen  fixation  studies, 


* 
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studies  were  conducted  at  the  Colombian  Agricultural 
Institute,  Medellin,  Colombia  to  become  familiar  with 
growing  in  its  natural  habitat*  A  preliminary 
selection  study  was  conducted  with  [)•  sandy  i  cense  for  larger 
and  smaller  eighth  node  terminal  leaflet  area.  Four  hundred 
plants  of  a  wild  D*  s_a ndw icense  population  were  selected  as 
the  base  population*  Several  characters  relating  to  the 
eighth  node  leaf  were  measured;  namely,  the  8th  leaflet 
length  and  width,  the  petiole  length,  the  internode  length 
(between  the  7th  and  8th  nodes),  and  the  shoot  length*  From 
these  characters  the  leaf  area,  the  length/width  ratio,  the 
petiole/ internode  ratio,  and  the  average  internode  length 
were  computed* 

Having  defined  the  selection  criterion  (8th  node 
terminal  leaflet  area)  and  the  related  characters,  the 
selection  was  initiated  by  selecting  five  plants  from  the 
high  and  low  'tails*  of  selection  criterion  distribution  and 
randomly  mating  the  plants  within  each  line*  For  the  first 
cycle  of  selection,  duplicate  data  for  leaf  size  were 


recorded  for  the  tenth  leaf  and  were  included  to  determi 


ne 


their  response  to  a  selection  pressure  applied  to  the  eighth 
node  terminal  leaflet*  One  hundred  plants  were  grown  for 
each  selection  cycle  and  selecting  five  plants  from  the 
'tails'  of  the  selection  criterion  distribution*  The  data 
for  all  characters  for  cycles  1  and  2  were  analysed  with  a 
one-way  analysis  of  variance  to  determine  whether  the 


characters  were  responding  to  selection* 
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(b)  Growth  patterns* 

To  determine  whether  the  bottle  jar  was  an  effective 
container  for  the  study  of  nitrogen  fixation  in  Destnodj  um  ,  a 
study  was  Initiated  to  observe  the  growth  of  Desmodlum  in  a 
bottle  jar*  The  fifteen  dlallel  genotypes  were  planted*  two 
plants  of  one  genotype  per  bottle  jar,  along  with  five 
(inbred  parents)  un inoculated  controls  to  determine  a  mean 
growth  curve*  The  experiment  was  conducted  during  the  winter 
season  and  the  primary  shoot  length  was  recorded  every 
twenty  days  for  100  days*  From  the  mean  growth  curve,  a 
sampling  point  for  the  nitrogen  fixation  studies  would  be 
determined • 

(c)  Linearity  of  acetylene  reduction* 

Hardy  ejt  &.1*  (  1968)  reported,  that  following  injection 
of  acetylene  into  a  vial  with  nodulated  soybean  roots, 
acetylene  was  reduced  to  ethylene  at  a  linear  rate  until 
approximately  90  minutes  following  injection*  To  verify  that 
acetylene  reduction  by  nodulated  Pe  smodj.  um  roots  occurred  at 
a  linear  rate,  four  Desmodium  species  ( Dsl ,  Dal,  Dp,  Du? 
Table  2)  of  plant  breeding  interest  were  inoculated  with 
strain  A  and  grown  for  8  weeks*  Two  plants  for  each  line 
were  bulk  harvested  and  prepared  for  the  acetylene  reduction 
assay*  Samples,  50  microlitres  each,  for  each  vial  were 
extracted  and  processed  at  10  min*  intervals* 

The  five  diallel  lines  were  also  prepared  for  the 
acetylene  reduction  assay  and  a  50  microlitre  sample  was 
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measured  every  10  min*  to  determine  if  fixation  was  linear# 

(d)  Bottle  jar  design* 

A  test  of  the  bottle  jar  design  was  initiated  by 
screening  river  sand  into  the  following  four  sizes: 

(a)  <  4.0  mm  (equivalent  to  that  of  Vincent  (1970)) 

(b)  1.0  —  3.0  mm 

(c)  1.0  —  1.5  mm 

( d )  <  1*0  mm 

The  sand  sizes  a  to  c  were  placed  in  the  initial  bottle  jar 
design  (Figure  2)  and  in  the  modified  design  (Figure  3)  for 
size  d*  The  sand  in  all  bottle  jars  was  compacted  about  the 
wick  to  ensure  capillary  movement  of  the  nutrient  solution* 
The  bottle  jars  were  arranged  into  a  randomized  block 
experiment  with  six  blocks*  planted  with  the  Dsl  line  and 
Inoculated  with  strain  B.  Three  seeds  were  planted  in  each 
bottle  jar  and  after  the  rock  mulch  was  applied*  each  cell 
was  thinned  to  the  two  most  vigorous  plants*  The  experiment 
was  harvested  after  eight  weeks  and  the  data  for  dry  weight 
and  shoot  length  were  recorded  for  each  cell  as:  (a)  the 
mean  of  the  two  plants*  and  (b)  the  difference  of  the  two 
plants  in  the  same  bottle  ( withln-cell  difference)*  The  cell 
mean  data  were  analysed  with  a  two-way  analysis  of  variance 
to  determine  significant  differences  for  dry  weight 
accumulation  among  sand  sizes* 
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2-  Strain-specific  ineffective  nodtulation  Interaction. 

S t r a i n— s pec i f i c  ineffective  nodulati op  interactions 
between  some  lines  of  I),  aga r  1  ne s  and  strain  B  were  found  in 
this  study.  To  determine  the  mode  of  inheritance  of  these 
interactions,  D.  §.a.adlLL£S.nS£.  (Dsl)  X  |).  apart  P  (Dal  and 
Da2  )  cross  plants  were  scored  for  green  or  yellow-chlorotic 
color  in  the  F2  and  backcross  (to  Dsl)  generations  to 
determine  the  segregation  ratios.  The  scored  results  were 
tested  with  the  chi-square  goodness  of  fit  to  an  expected 
segregation  ratio. 

Two  hundred  and  fifty  F2  seeds,  150  backcross  seeds  and 
10  seeds  each  for  parental  controls  were  planted  in  bottle 
Jars  and  inoculated  with  5  ml  of  strain  B  to  each  seed.  At 
eight  weeks  the  plants  were  scored  for  green  or  yellow  color 
and  a  random  sample  of  seventy  ineffectively  nodulated  F2 
plants  were  removed  from  the  bottle  Jars  and  inspected  for 
nodulation  characteristics.  Following  inspection  these 
plants  were  planted  in  sand  and  fertilized  with  a  source  of 
nitrogen  to  determine  if  nitrogen  was,  in  fact,  the  limiting 
criterion. 

Since  some  lines  of  D.  apar jnes  display  the  spotted 
leaf  character,  all  the  plants  were  also  scored  for  this 
character  and  tested  for  a  chi-square  goodness  of  fit  to  an 
expected  segregation  ratio.  Joint  segregation  tests  were 
determined  for  ineffective/effective  nodulation  and 


spotted/plain  leaves. 
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3*  Experiment  I* 

The  diallel  cross  was  arranged  in  Experiment  I  as  a 
randomized  block  design  with  four  blocks  (replicates)  and 
fifteen  bottle  jars  per  block#  The  bottle  jars  were  randomly 
assigned  the  fifteen  diallel  genotypes  and  were  randomized 
on  a  30  cm  grid  within  each  block.  Similar  bottle  jars 
(original  design,  Figure  2)  were  used  within  each  block  with 
blocks  1  and  4  having  similar  bottle  shapes  while  blocks  2 
and  3  were  different.  Each  bottle  jar  was  planted  with  two 
surface  sterilized  seeds  and  each  plant  was  inoculated  with 
a  5  ml  suspension  of  strain  B  following  emergence.  Blocks 
were  planted  at  differing  planting  dates  due  to  space 
limitations  and  to  facilitate  the  acetylene  reduction 
analysis.  Experiment  I  was  conducted  during  the  winter 
season  at  the  University  of  Sydney;  therefore,  the  average 
glasshouse  temperature  was  near  21  C. 

Each  block  was  harvested  after  a  growth  period  of  10 
weeks  and  data  relating  to  shoot  length,  dry  weight 
accumulation  and  the  acetylene  reduction  assay  were  recorded 
as  the  mean  of  the  two  plants  in  each  bottle  jar. 

Although  the  inbred  lines  were  not  selected  for  the 
characters  studied,  the  lines  have  been  considered  as  a 
fixed  sample  of  ]£•  sandwl  cense  lines.  The  data  were  analysed 
under  the  assumptions  of  the  model  I,  method  2  (Griffing, 
1956)  and  Analysis  II  (  Gardner  and  Eberhart,  1966  )  diallel 
design  to  establish  the  combining  ability  effects  (Griffing, 
1956)  and  genetic  constants  (Gardner  and  Eberhart,  1966  ), 


.. 
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respectively* 

4*  Experiment  II* 

The  Experiment  II  dialiel  cross  was  arranged  as  a 
randomized  bloclc  design  with  four  blocks  (  replicates  )  and 
fifteen  bottle  jars  per  block*  The  bottle  jars  were  randomly 
assigned  the  fifteen  dialiel  genotypes  and  were  randomized 
on  a  30  cm  grid  within  each  block*  Similar  bottle  jars 
(original  design,  Figure  2)  were  used  within  each  block  with 
blocks  1  and  2  having  one  bottle  shape,  whereas  blocks  3  and 
4  had  another*  Each  bottle  jar  was  planted  with  two  surface 
sterilized  seeds  and  each  plant  was  inoculated  with  a  5  mi 
suspension  of  strain  A  following  emergence*  Blocks  were 
planted  at  different  planting  dates  due  to  space  limitations 
and  +o  facilitate  the  acetylene  reduction  analysis* 
Experiment  II  was  conducted  during  the  spring  season  at  the 
University  of  Sydney;  therefore,  the  average  glasshouse 
temperature  was  near  24  C* 

Each  block  was  harvested  after  a  growth  period  of  10 
weeks  and  data  relating  to  dry  weight  accumulation  and  the 
acetylene  reduction  assay  were  recorded  for  each  plant*  The 
data  were  analysed  under  the  assumptions  of  the  model  I, 
method  2  (Griffing,  1956)  and  Analysis  II  (Gardner  and 
Eberhart,  1966)  dialiel  design  to  establish  the  combining 
ability  effects  and  genetic  constants* 
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5.  Experiment  III. 

Experiment  III  was  grown  as  a  factorial  experiment  with 
two  strains  (A  and  B),  two  blocks  for  each  strain,  and 
bottle  jars  per  block  with  three  plants  for  one 
genotype  planted  in  each  bottle  jar.  The  fifteen  diallel 
genotypes  were  randomly  assigned  to  bottle  jars  which  were 
randomized  on  a  30  cm  grid  within  each  block.  Each  block  was 
randomly  assigned  a  bench  location  and  blocks  inoculated 
with  different  strains  were  separated  on  the  bench  to  avoid 
contamination.  The  total  experiment  was  planted  on  one  day 
et  the  University  of  Sydney  during  the  summer  season.  The 
mean  glasshouse  temperature  for  this  experiment  was 
approximately  26  C. 

Similar  types  of  bottle  jars  (modified  design,  Figure 
3)  were  autoclaved  in  two  groups  such  that  autoclaving  has 
been  confounded  with  blocks. 

During  the  early  juvenile  leaf  stage.  Experiment  III 
was  Invaded  by  insect  larvae.  Since  leaf  loss  resulted,  the 
experiment  was  sprayed  with  an  insecticide  and  each  plant 
was  clipped  back  to  the  cotyledonary  node.  After  10  weeks, 
all  plants  were  harvested  and  measurements  for  the  acetylene 
reduction  assay  were  taken  for  strain  A  on  one  day  and  for 
strain  B  the  following  day.  One  block  for  each  strain  was 
analysed  during  the  late  morning  wi th  the  other  block 
analysed  during  the  early  afternoon.  Dry  weight  accumulation 
were  also  recorded  for  the  shoot  portion  of  each  plant.  Due 


to  the  death  of  several  plants,  the  two  plants  with  the 
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largest  dry  weight  in  each  bottle  Jar  were  retained  for 
diallel  analysis*  The  data  were  analysed  under  the 
assumptions  of  the  model  I,  method  2  (Griffing,  1956)  and 
Analysis  II  (Gardner  and  Eberhart,  1966)  diallel  design  to 
establish  the  combining  ability  effects  and  genetic 
constants* 

6*  Combined  analysis* 

The  same  diallel  cross  was  used  in  the  preceding 
diallel  experiments?  therefore,  the  data  may  be  pooled  for  a 
comprehensive  analysis*  Data  for  individual  plants  were 
available  for  block  4  in  Experiment  I  and  all  blocks  in 
Experiments  II  and  III  giving  nine  blocks  in  total*  Three  of 
these  blocks  were  inoculated  with  strain  B  while  the 
remaining  six  were  noculated  with  strain  A*  Since  strains 
have  been  considered  as  environmental  effects  for  this 
study,  strains  have  been  pooled  with  blocks  for  the  combined 
analysis*  The  data  were  analysed  under  the  assumptions  of 
the  model  I,  method  2  (Griffing,  1956)  and  Analysis  II 
(Gardner  and  Eberhart,  1966)  diallel  design  to  establish  the 
combining  ability  effects  end  genetic  constants* 

7*  Experiment  IV* 

A  cross  between  the  two  inbred  lines,  Dsl  x  Dal  (Table 
2)  was  made  creating  the  six  populations:  (Pi(Dsl),  P2<  Dal  ), 
Fit  t  BC i  (Pi  x  Fi  )  ,  and  BC2  ( P2  x  Ft)*  Experiment  IV  was 


grown  during  the  spring  season  at  the  University  of  Sydney 
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as  a  randomized  block  design  with  seven  blocks  and  twelve 
bottle  jars  per  block*  All  bottle  jars  were  of  the  original 
design  (Figure  2)  and  blocks  1  to  5  had  600  ml  reservoirs 
(Figure  2  )  while  blocks  6  and  7  had  2  1  reservoirs*  The  six 
populations  were  allotted  to  the  twelve  bottle  jars  such 
that  the  Plt  p2  and  Fi  populations  received  one  bottle  jar 
each,  the  F2  population  received  five,  and  the  backcross 
populations  received  two  each*  The  bottle  jars  were 
randomized  within  each  block  and  all  blocks  were  planted  on 
the  same  day  with  surface  sterilized  seeds*  After  emergence, 
each  plant  was  inoculated  with  5  ml  of  a  suspension  of 
strain  A*  The  experiment  was  harvested  after  seven  weeks 
growth  with  dry  weight  accumulation  and  the  acetylene 
reduction  assay  being  recorded  for  each  plant*  The  analysis 
of  variance  of  population  means  (Hayman,  1S58)  was  conducted 
as  a  two-way  analysis  of  variance  with  the  population  means 
within  each  block  constituting  the  data  which  was  analysed 
as  outlined  in  Section  E2*  Although  Hayman  did  not  partition 
eplstatic  effects,  a  study  was  initiated  to  investigate  the 
feasibility  of  such  a  partition  (see  Results  and  Discussion 
sec •  A  5  )  • 

8*  Experiment  V* 

The  six  populations  of  the  Ds2  x  Dsp  (Table  2)  cross 
were  grown  in  Experiment  V  as  a  randomized  block  design  with 
five  blocks  and  fifteen  bottle  jars  per  block*  All  bottle 
jars  were  of  the  modified  design  (Figure  3)  and  all  had  600 
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oil  reservoirs  (Figure  3).  The  six  populations  were  allotted 
to  the  fifteen  bottle  Jars  such  that  the  Pt ,  P2  and  Fi 
populations  received  one  bottle  jar  each,  the  F2  population 
received  six,  and  the  backcross  populations  received  three 
each*  The  bottle  jars  were  randomized  within  each  block  and 

blocks  were  planted  on  the  same  day*  After  emergance  of 
each  plan t  was  inoculat ed  with  a  5  ml  suspen sion  of  strain 
B*  The  experiment  was  harvested  after  ten  weeks  growth* 

Dry  weight  accumulation  and  the  acetylene  reduction 
assay  were  recorded  for  each  plant  and  the  analysis  of 
variance  of  population  means  ( Hayman ,  1958)  was  conducted  as 
a  two-way  analysis  of  variance  with  the  population  means 
within  each  block  constituting  the  data  which  were  analysed 
as  outlined  in  Section  E2* 

9*  Experiment  VI* 

Experiment  VI  consisted  of  the  six  populations  from  the 
cross  between  the  Inbred  lines  Ds2  and  Ds5  (Table  2)  and  was 
grown  in  a  growth  chamber  at  the  University  of  Alberta* 
Three  plants  of  one  population  were  grown  in  e  ten  cm  pot 
and  inoculated  at  planting  with  five  ml  of  a  suspension  of 
strain  A*  The  pots  were  arranged  into  four  blocks;  two  of 
which  had  sand  of  very  fine  (  *15ram)  particle  size,  one  with 
fine  (  *25mm)  sand,  while  the  remaining  block  was  a  50% 
mixture  of  each  sand  size*  Within  each  block  the  populations 
were  allotted  such  that  5  pots  each  were  allotted  to  the  Pi 
(Ds2),  P2  (  Ds5  ),  and  Fi  populations  while  20  pots  were 
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allotted  to  the  the  F2  population  and  10  pots  each  to  the 
backcross  populations*  All  pots  were  wetted  with  the 
nutrient  solution  prior  to  planting*  The  experiment  was 
maintained  at  a  temperature  of  23  C  and  harvested  after  11 
weeks*  Data  recorded  for  each  plant  were  the  acetylene 
reduction  assay*  shoot  fresh  and  dry  weight*  root  fresh 
weight,  the  cotyledon  length  and  width  (first  node)*  and  the 
width  of  the  juvenile  leaf  (second  node)  since  the  juvenile 
leaf  Is  nearly  circular*  The  data  were  analysed  on  an 
individual  plant  basis  as  outlined  in  Section  E2* 

E*  Genetical  analyses* 

1*  Dlallel  cross* 

A  diallel  cross  is  a  set  of  all  possible  crosses  among 
a  set  of  inbred  lines  and  the  diallel  fable  is  a  square 
matrix  of  dimension  equal  to  the  number  of  parents*  The  main 
diagonal  represents  the  inbred  lines  while  the  data  above 
and  below  the  main  diagonal  represents  the  Fj  and  reciprocal 
Fi  crosses,  respectively*  In  the  absence  of  reciprocal 
effects  the  matrix  is  symmetrical*  A  diallel  table  may  be 
formed  by  four  methods  from  the  components  of  the  square 
matrix  (Griffing,  1956)* 

The  diallel  analysis  of  variance  of  the  data  in  the 
diallel  table  has  been  reviewed  and  discussed  by  several 
workers  (Sprague  and  Tatum*  1942;  Hayman *  1954;  Griffing* 
1956;  Cockerham,  1963;  Gardner  and  Eberhert*  1966)  and  the 
relationship  among  the  four  methods  and  the  analyses  for 
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each  method  reviewed  by  Baker,  1978. 

The  interpretation  of  the  diallel  analysis  of  variance 
depends  upon  the  model  under  investigation.  Griff ing  (1956) 
reported  four  differing  models  underlying  the  analysis  of 
variance.  The  models  are: 

Model  I.  Fixed  treatment  and  block  effects 
Model  II.  Random  treatment  and  block  effects 
Model  A.  Random  treatment  and  fixed  block  effects 
Model  B.  Fixed  treatment  and  random  block  effects. 

The  basic  difference  between  these  models  relates  to  the 
difference  in  the  interpretation  of  the  fixed  and  random 
treatment  (genotype)  effects.  The  fixed  model  assumes  that 
the  inbred  lines  and  their  Ft  crosses  constitute  the  total 
population  about  which  inferences  are  to  be  made  and 
estimates  combining  ability  effects  and  variances  for  the 
inbred  lines.  The  random  model  assumes  that  the  inbred  lines 
are  a  random  sample  from  an  ancestral  population  about  which 
inferences  are  to  be  made.  The  random  model  estimates  the 
genetic  and  non-genetic  components  of  the  population 
variance • 

Eayman  (1954)  presented  several  assumptions  underlying 
a  diallel  analysis  which  must  be  satisfied  before  the 
conclusions  may  be  considered  reliable.  These  assumptions 
er  e : 

(a)  that  the  parents  are  homozygous, 

(b)  there  is  diploid  segregation, 


(c)  there  are  no  reciprocal  effects. 


* 


' 
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<d)  there  is  no  multiple  allelism, 

(e)  the  genes  are  independently  distributed,  and 

(f)  there  is  independent  action  of  non-allelic  genes 
(epistatic  effects)* 

^  ^  1956  )  re  po  rted  that  the  i  nde  penden  t  action  of 

the  non-allelic  genes  was  the  assumption  of  major 

importance*  The  dialiel  cross  among  Inbred  lines  of  D* 

g S,jy dw j- c £ o. sje  used  in  this  study  satisfied  assumptions  (a)  and 
(b)  while  assumptions  (c)  and  (d)  were  considered  to  be  of 
minor  importance*  Assumptions  (e)  and  (f)  concerning  linkage 
and  epistatic  effects,  respectively,  were  of  major  concern 
in  the  present  investigation  since  Despiod jL um  has  never  been 
tested  for  such  effects. 

The  analysis  of  the  acetylene  reduction  assay,  the  dry 
weight  accumulation  characters  and  other  characters  were 
based  on  a  Model  I,  method  2  dialiel  analysis  as  outlined  by 

Griffing  (1956)  and  modified  to  Include  Analysis  II  of 

Gardner  and  Eberhart  (1966)*  The  Model  I,  method  2  dialiel 
analysis  used  in  this  study  implies  a  statistical  rather 
than  a  genetic  interpretation  (Baker,  1978)  to  partition  the 
genetic  component  sum  of  squares  and  to  estimate  genetic 
constants  (Gardner  and  Eberha r t , 196 6 )  and  combining  ability 
effects  (Griffing,  1956)*  A  listing  of  the  Fortran  computer 
program  employed  in  this  study  for  the  dialiel  analysis  of 


the  data  is  presented  in  Appendix  B* 


*  - 
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2.  Cross  between  two  inhrect  lines. 

An  analysis  of  a  cross  between  two  inbred  lines  as 
outlined  by  Mather  and  Jinks  (1971)  may  consist  of  six  or 
more  populations.  The  six  populations  of  the  minimum  design 
are:  parent  1  (Pi),  parent  2  (P2),  the  initial  cross  (  Fx  ), 
the  segregating  generation  (F2),  and  the  backcrosses  of  the 
Fx  to  both  parents  ( B!  and  B2  )  .  For  this  fixed  model, 
estimates  of  the  type  of  gene  action  may  be  obtained  from 
first  degree  statistics  (Mather  and  Jinks,  1971). 

Estimates  of  the  magnitude  and  significance  of  the 
genetic  parameters  were  calculated  from  the  population  means 
with  the  formulae  as  presented  by  Mather  and  Jinks  (1971). 
The  equations  for  the  calculation  of  the  genetic  parameters 
are  : 


m  = 

.5P!  + 

*  5P2  +  4F2  -  2Bi 

-  2B2 

a  = 

•  SPi  - 

•  5  P2  t  Bi  —  B2 

d  = 

1.5Pi  - 

■  1 • 5P 2  —  F i  —  8F2 

t  6Bi  t  6B2 

ea 

II 

1 

4*. 

Til 

ro 

+  2Bi  +  2B2 

ad 

=  -Pi  + 

P2  +  2B !  -  2B2 

dd  =  Px  +  P2  +  2F i  +  4F 2  -  4 Bi  -  4B2 
whe  re 

m  =  F2  mean 
a  =  additive  effects 
d  =  dominance  effects 

aa  =  additive  X  additive  interaction 
ad  =  additive  X  dominance  interaction 
dd  =  dominance  X  dominance  interaction. 


' 
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These  equations  allow  for  an  estimate  of  the  magnitude 
of  the  parameters  to  be  made  with  the  significance  of  each 
being  tested  by  the  t-test  of  the  estimate  to  its  standard 
deviation  (Mather  and  Jinks,  1971  )•  The  variance  of  each 
parameter  is  determined  from  the  variances  of  the  six 
populations.  For  example,  an  estimate  of  the  variance  of  d 
is  obtained  from  the  following  equation! 

Vd  =  2 • 25VP !  +  2. 25VP2  +  VFt  +  64VF2  +  36VBt  +  36VB2 

where  Vd  =  the  variance  of  the  dominance  deviations,  VPt  = 

the  variance  of  Pj ,  VP2  =  the  variance  of  P2 ,  etc. 

An  alternate  analysis  of  population  means  is  the 
analysis  of  variance  of  the  population  means  (Hayman,  1958; 
Voigt  et  si • t  1966).  To  estimate  the  genetic  parameters  of 
the  model  a  coded  matrix  was  constructed  (Mather  and  Jinks, 


1971  ). 

The 

code  d 

- —  -  -  —  — —  J 

matrix  is  as  follows: 

Mean 

a 

d 

aa 

a  d 

dd 

Pi 

1 

0 

1 

0 

0 

P2 

1 

0 

1 

0 

0 

Fi 

0 

1 

0 

0 

1 

f2 

0 

1/2 

0 

0 

1/4 

Bi 

1/2 

1/2 

1/4 

1/4 

1/4 

I2 

1/2 

1/2 

1/4 

-1/4 

1/4 

The  matrix 

may 

be  weighted 

by  dividing  each  row  by  the 

standard  error  of  the  population  mean.  The  coded  matrix  was 
subjected  to  a  stepwise  least  squares  analysis  (Draper  and 
Smith,  1966)  which  has  been  restricted  in  this 


investigation.  The  genetic  parameters  of  the  model  were 


\r 

V 
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estimated  on  the  assumption  that  within  locus  effects  are  of 
more  importance  than  between  loci  interactions  in  a  genetic 
model,  i«e«  a  simple  model*  Therefore  the  variance  among  the 
population  means  was  first  analysed  for  additive  and 
dominance  effects  before  estimating  the  interaction 
parameters  (see  Results  and  Discussion,  Section  A5  )•  The 
significance  of  the  parameter  mean  squares  was  tested  by  the 
F-ratio  with  the  error  term  from  the  analysis  of  variance  of 
population  means  as  the  denominator  (Voigt  et  &l.  ,  1966  )* 

Properly  coding  the  genotypes,  replications,  etc.  prior 
to  data  entry  allows  the  usage  of  the  SPSS  statistical 
analysis  package  (Nle,  et  al, ,  1975)  for  analysing  the  data 
for  a  cross  between  two  inbred  lines.  SPSS  is  particular  1 ly 
suited  to  the  recoding  of  the  data  or  the  creation  of  new 
variables  to  assist  in  the  complete  analysis  and 
interpretation  of  the  data.  The  partitioning  of  the 
population  sum  of  squares  into  its  genetic  components  is  not 
available  in  SPSS;  however,  an  interactive  Fortran  computer 
program,  listed  in  Appendix  C,  was  developed  for  this 
analysis.  Program  HAYMAN  follows  the  analytical  procedure 
for  the  analysis  of  variance  of  population  means  as  defined 
by  Hayman  (1958)  and  earlier  discussed.  HAYMAN  is  suited  for 
genetic  analysis  of  any  type  of  statistical  design  whether  a 
one-way  or  a  multi— factor! al  analysis  of  variance  or 
covariance  since  the  population  sum  of  squares  is 
partitioned  into  genetic  components  by  the  correlation  form 


of  the  stepwise  regression  analysis. 


This 


regress ion 
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analysis  is  independent  of  scale  or  design  (Draper  and 
Smithy  1966  )• 

(3)  Principal  component  transformations. 

Total  nitrogen  fixation  has  been  estimated  in  this 
study  by  many  measures  and  correlations  among  the  measures 
would  be  expected*  When  variables  are  highly  correlated  the 
communal  variation  within  the  correlation  matrix  can  be 
summarized  with  a  principal  components  transformation 
(Winer,  1971)  to  construct  "variables"  which  estimate  the 
variation  along  the  axes  of  principal  variation.  This 
transformation  is  also  called  principal  components  analysis 
(Overall  and  Klett,  1972)  where  the  primary  emphasis  is  to 
define  "composite  or  factor-score  variates  that  have  certain 
desirable  statistical  properties".  Principal  components 
analysis  may  be  divided  into  three  stages  of  analysis* 
Firstly  the  analysis  begins  with  the  calculation  of  a 
correlation  or  a  covariance  matrix*  Then  the  correlation  or 
covariance  matrix  is  mathematically  reduced  to  eigenvectors 
(factor  loadings)  and  eigenvalues  (the  amount  of  variance 
accounted  by  each  eigenvector)*  Generally  a  correlation 
matrix  with  unity  (total  variances)  along  the  main  diagonal 
is  structured  because  this  matrix  is  scale  independent 
(Overall  and  Klett,  1972)*  Finally,  factor  score 
coefficients  are  calculated  from  the  factor  loading  matrix 
with  a  straight-forward  mathematical  procedure  analogous  to 


a  regression  of  an  unknown  dependent  variable  on  the 
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original  variables*  A  composite  variable  is  a  linear 
combination  of  factor  score  coefficients  of  the  variables. 
The  first  composite  variable  accounts  for  greatest  amount  of 
common  variation  in  the  data*  The  second  composite  variable 
is  orthogonal  to  the  first  and  accounts  for  the  next  largest 
portion  of  the  common  residual  variation.  The  remaining 
composite  variable  accounts  for  the  remaining  variation. 
Principal  components  analysi s»  being  primarily  concerned 
with  construction  of  composite  variable  scores,  normalizes 
the  eigenvectors  by  converting  the  eigenvector  matrix  to  one 
of  correlations  among  the  original  variables  and  the  derived 
factors.  Such  a  normalization  is  of  interest  in  determining 
what  the  principal  axes  represent.  From  the  normalized 
eigenvectors,  composite  variables  are  defined  which  maximize 
existing  differences  among  individual  cases. 

Principal  components  analysis  has  been  used  in  agronomy 
by  Denis  and  Adams  (  1978)  to  relate  plant  variables  to  yield 
of  Phase ol us  vulgaris  L.  (common  bean)  and  by  Hussaini,  g.t 
&l.  (1977)  to  define  major  groupings  within  the  world 
collection  of  Eleusi ne  coracana  L.  (finger  millet).  However, 
principal  components  analysis  has  not  been  employed  as  a 
method  of  analysing  differences  among  the  population  means. 
Since  the  principal  axis  of  variation  is  expected  to  be 
nitrogen  fixation,  it  would  be  pertinent  to  define  this  axis 
and  investigate  its  variation. 

The  principal  axis  of  variation  in  this  study  was 


scaled  to  a  distribution  with  a  mean  and  standard  deviation 
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of  1/3  that  of  the  acetylene  reduction  assay.  This  scale  iras 
used  because  nitrogenasp  theoretically  reduces  three  times 
the  amount  of  acetylene  than  N2  (Burris,  1974).  Following 
scaling,  the  composite  N2  fixation  estimate  was  treated  as 
another  character  for  the  analyses.  The  interpretation  of 
the  principal  axis  of  variation  will  be  restricted  to  the 
investigation  of  the  differences  among  population  means 
rather  than  attempting  to  estimate  the  amount  of  nitrogen 
fixed. 


RESULTS  AND  DISCUSSION 

A»  Methodology 

To  achieve  the  ultimate  goal  of  lines  or  hybrids  with 
greater  total  nitrogen  fixation  it  is  necessary  to  attain  an 
understanding  of  the  genetic  mechanisms  operating  in  the 
host*  In  this  study,  the  area  of  genetic  mechanisms  of 
nitrogen  fixation  has  been  divided  into  three  levels  (Table 
^  ^ this  generalized  section  of  methodology  several 
studies  were  undertaken  both  initially  and  during 
experimentation  to  establish  or  to  investigate  various 
aspects  of  methodology* 

1*  Selection  for  leaf  size* 

A  preliminary  study  was  initiated  in  Des medium  to 
become  familiar  with  this  species  growing  in  its  natural 
habitat*  This  study  involved  applying  e  selection  pressure 
to  the  eighth  node  terminal  leaflet  area  to  determine  its 
response  to  selection*  The  frequency  distribution  of  the 
eighth  node  terminal  leaflet  sizes  for  the  base  population 
are  presented  in  Figure  4*  The  distribution  indicated  that 
considerable  variation  was  present  in  the  wild  population* 
Principal  components  analysis  of  the  base  population 
characters  suggested  that  most  characters  were  related  to 
one  factor  i*e*  a  leaf  size  factor  (Factor  1,  Table  3)  while 
the  ratio  characters  are  individually  related  to  other 
f ac  tors* 

Data  collected  for  the  tenth  node  leaf  for  a  Dsp  (Table 
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Figure  4.  Frequency  distribution  of  the  eighth  node  terminal  leaflet 
area  for  the  base  population. 
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Table  3*  Eigenvectors  for  factors  from  the  base  population. 


Factor 


Character 

1 

2 

3 

Terminal  leaflet  length! 8*  ) 

0.94 

0.07 

0.01 

Terminal  leaflet  width!  8) 

0.91 

0.05 

-.34 

Terminal  leaflet  area! 8  ) 

0.94 

0.06 

-.  17 

Length/width  ratio!  8) 

-.04 

0.02 

0.95 

Petiole  length! 8) 

0.85 

0.28 

0.  11 

Internode  length!  7-8) 

0.74 

-.54 

0.11 

Petiole/internode!  8) 

0.  10 

0.96 

0.  04 

Shoot  length 

0.81 

-.33 

0.25 

Average  inter node  length 

0.81 

-.41 

0.  25 

^node  number. 


_ 


- 
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2)  population  indicated  that  leaf  measurements  were  related 
to  one  factor  while  plant  weights  related  to  another  (Table 

4) *  Ratio  characters  again  were  related  each  to  a  single 

f ac  tor* 

Principal  components  analysis  of  the  cycle  1  data  still 
suggested  leaf  size  characters  were  related  to  the  first 
factor  while  ratio  characters  were  related  to  another  (Table 

5) .  The  means  and  standard  deviations  for  each  character  for 
cycle  1  are  presented  in  Table  6.  The  base  population  was 
found  to  have  a  significantly  larger  eighth  terminal  leaflet 
area  than  either  the  high  or  low  selection  lines  with  the 
exception  of  the  ratio  characters  where  the  reverse  was 
true*  Means  and  standard  deviations  for  the  factor  score 
variables  indicated  the  same  patterns  as  the  selection 
criterion  (Table  7)* 

Cycle  2  data  also  indicated  the  same  pattern  between 
the  characters  and  the  factors  (Table  8)*  The  response  to 
selection  by  the  high  line  was  found  to  be  negative  for 
cycle  1;  however,  the  high  line  mean  was  found  to  be 
significantly  larger  than  the  base  mean  suggesting  a 

response  to  selection  (Table  9)*  The  low  line  mean  was  found 

to  be  significantly  lower  than  the  base  mean  while  the 

factor  score  means  also  indicated  that  the  response  to 
selection  away  from  the  base  population  was  significant 

(Table  10)*  The  analyses  suggested  that  the  leaf  size 
characters  were  highly  related  as  indicated  by  their  factor 
loadings  to  factor  1  and  similarity  with  which  they 


. 
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Table  4*  Eigenvectors  tor  factors  from  the  Dsp  population* 


Charac  ter 

1 

Factor 

2  3 

4 

Terminal  leaflet  leng th( 1 0*  ) 

-•04 

0.94 

-.07 

0.28 

Terminal  leaflet  length( 10  ) 

-.  06 

0.97 

-.07 

-.12 

Terminal  leaflet  area(lO) 

-•06 

0.97 

-.10 

0.10 

Length/width  ratlo(lO) 

0*05 

0.18 

0.00 

0.96 

Petiole  length  (10) 

0*  14 

0.88 

0.15 

0.  11 

Internode  length  (9-10) 

-.06 

0.73 

-.57 

0.07 

Pet  i  ole/ inte  mode  ratio(lO) 

0*19 

-.06 

0.93 

0.  01 

Shoot  length 

0*  91 

0.03 

-.21 

0.  17 

SQRT(  no*  of  nodes  ) 

0.76 

-.36 

0.  16 

-.10 

Average  internode 

0*73 

0.26 

-.32 

0.28 

SQRT(  no*  of  tillers) 

0*76 

-.12 

0.40 

.00 

Shoot  weight 

0*95 

0.  01 

0.17 

-.03 

Leaf  weight 

0.88 

-.07 

0.19 

0.01 

Total  weight 

0*93 

-.04 

0.19 

0.  00 

Leaf/shoot  ratio 

-.72 

-.  14 

0.01 

0.11 

*node  number* 

SORT  —  square  root* 


- 
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Table  5*  Eigenvectors  for 

f ac  tors 

from  cycle 

1  . 

Factor 

Charac  ter 

1 

2 

3 

4 

Terminal  leaflet  len?th( 8*  ) 

0.88 

-.27 

-.01 

0.22 

Terminal  leaflet  width(8) 

0.89 

-.26 

-.04 

-.20 

Terminal  leaflet  area( 8  ) 

0.91 

—  .  28 

-.04 

0.02 

Length/width  ratlo(8) 

-.04 

-.03 

0.07 

0.84 

Petiole  length  (8) 

0.  74 

0.42 

0.18 

0.  1 1 

Internode  length  (7-8) 

0.36 

-.76 

0.28 

0.10 

Pet iole/ internode  ratio(8  ) 

0.14 

0.92 

-.  15 

-.03 

Terminal  leaflet  length! 10) 

0.91 

0.17 

0.12 

0.  17 

Terminal  leaflet  width!  10) 

0.88 

0.20 

0.05 

-.31 

Terminal  leaflet  area(10) 

0.93 

0.20 

0.  08 

-.07 

Length/width  ratlo(lQ) 

0.01 

-.  05 

0.10 

0.  85 

Petiole  length!  10  ) 

0.56 

0.33 

-.39 

-.01 

Internode  length(  9  — 10) 

0.46 

0.02 

0.  75 

0.09 

Petiole/ internode  ratio(lO) 

0.10 

0.27 

-.93 

-.07 

*node  number# 


>  •- 


Table  6*  Means  and  standard  deviations  from  a  one-way  analysis  of  variance  of  the 
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Table  7.  Means  and  standard  deviations  from  a  one-way  analysis  of 
variance  of  the  factors  for  cycle  1# 
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Table  8*  Eigenvectors  for 

factors  from 

cycle  2. 

Factor 

Cha  rac ter 

1 

2 

3 

Terminal  leaflet  length( 8*  ) 

0.95 

-.12 

0.  05 

Terminal  leaflet  length(8) 

0.  92 

-.09 

-.24 

Terminal  leaflet  area( 8  ) 

0.94 

-.12 

-.10 

Length/width  ratio(8) 

-.17 

-.03 

0.  84 

Petiole  length  (8) 

0.  91 

-.06 

-.11 

Internode  length  (  7-8  ) 

0.56 

-.75 

-.11 

Pet  iole/ in te  mode  ratio(  8  ) 

0.16 

0.89 

0.04 

Terminal  leaflet  length(lO) 

0.96 

0.02 

0.  06 

Terminal  leaflet  width(  10) 

0.93 

0.01 

-.30 

Terminal  leaflet  area(10) 

0.96 

0.00 

-.13 

Length/width  ratio(lO) 

-.14 

0.05 

0.  86 

Petiole  length  (  10  ) 

0.  84 

0.21 

-.20 

Internode  length  (9—10) 

0.64 

-.63 

-.10 

Petiole/lnternode  ratio(lO) 

0.11 

0.90 

-.11 

*node  number* 
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Table  9.  Means  and  standard  deviations  from  a  one-way  analysis  of  variance  of  the 
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Table  10. 


Means  and  standard  deviations  from  a  one-way 
analysis  of  variance  of  the  factors  for  cycle  2. 


Factor 


Selection  line 

n 

1 

2 

3 

High 

95 

0 #74±.98a* 

-0.47±0. 97b 

—0 • 20±0. 80b 

Bas  e 

99 

0.00±. 79b 

0. 28± 1 • 12a 

-0. 23±0. 97b 

Low 

95 

-0.74±.59c 

0. 17±0. 70a 

0. 44±1. 06a 

*means  followed  by  the  same  letter  are  not  significantly 

level  for  the  new  multiple  range  test 
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responded  to  a  selection  pressure  upon  one  of  them.  The 
selection  response  for  the  eighth  node  terminal  leaflet  and 

the  first  principal  factor  over  two  cycles  has  been 
summarized  in  Figure  5. 

2#  Growth  patterns. 

An  initial  study  was  conducted  with  the  crosses  among 
£•  SanMceflse  lines  inoculated  with  strain  B  to  observe  the 
growth  of  Pefimodjj^n  in  a  bottle  Jar  and  to  observe  the 
general  growth  pattern  over  a  period  of  time.  The  mean 
growth  patterns  of  the  inoculated  and  uninoculated  plants 
are  presented  in  Figure  6.  During  the  winter  season,  the 
transition  from  lag  to  log  growth  phase  occurred  at 
approximately  the  10th  week  for  the  inoculated  plants.  Since 
this  transition  occurred  later  in  the  winter  season,  the 
10th  week  was  chosen  as  an  appropriate  sampling  time  for  a 
study  of  the  relationship  between  dry  weight  accumulation 
and  the  acetylene  reduction  assay.  The  uninoculated  controls 
were  yellow  and  chlorotic  and  their  height  did  not  Increase 
beyond  the  height  at  juvenile  stage.  This  chlorotic 
condition  suggested  that  the  bottle  Jar  was  an  effective 
assembly  for  the  study  of  symbiotic  nitrogen  fixation  in 
Desmodlum . 

Diallel  tables  for  several  characters  measured  on  the 
inoculated  plants  of  the  above  study  were  constructed  and 
partloned  into  genetic  component  mean  squares  which  are 
expressed  as  a  percentage  of  the  genotype  mean  square  and 
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Figure  5.  Selection  response  of  the  high  and  low  lines  for  the  selection 
variable  and  the  first  principal  factor. 
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Figure  6.  Mean  growth  curve  for  the  fifteen  genotypes  of  D.  sandwicense 
inoculated  with  Rhizobium  strain  B  and  the  uninoculated  controls. 


65 


presented  in  Table  11.  The  pattern  which  emerged  was  the 
consistency  with  which  general  effects  and  average  heterosis 
values  exceeded  the  mean  square  for  genotypes  while  the 
other  components  were  less.  These  data  suggested  the 
presence  of  heterosis  and  additive  genetic  variance. 

It  was  observed  in  the  above  study  that  the  parental 
lines  appeared  to  differ  in  growth  before  the  initiation  of 
the  log  growth  phase  and  remained  different  until  the  study 
terminated.  To  confirm  these  results  the  primary  shoot 
length  was  recorded  weekly  for  the  parental  plants  in  block 
1  of  Experiment  I.  The  growth  curves  for  the  five  inbred 
lines  are  presented  in  Figure  7.  The  differences  between 

these  lines  became  evident  as  early  as  the  sixth  week  after 
plan  ti ng. 

Hybrid  vigor  for  forage  yield  has  been  generally 
considered  to  be  absent  in  Desffiodliim  when  adequately 
supplied  with  nitrogen( McWhir ter ,  1969;  Fotar  and  Chow, 
1971).  The  heterosis  present  in  the  estimates  of  symbiotic 
nitrogen  fixation  must  be  due  solely  to  the  processes  of 
nitrogen  fixation  because  the  non-nodula t ed  controls  did  not 
accumulate  dry  matter.  Therefore,  the  bottle  jar  is  an 
effective  tool  for  creating  a  nitrogen  limiting  environment 
for  expressing  differences  in  nitrogen  fixing  potential 
among  population  genotypes# 

3#  Linearity  of  acetylene  reduction. 


Acetylene  reduction  samples  were  prepared  to  determine 


66 


v 

0)  pJ 


c 

1 

* 

a) 

1 

4» 

© 

© 

1 

r*4 

£ 

E 

£ 

d 

>> 

do 

4* 

4* 

•H 

© 

0 

*0 

© 

a 

C 

H 

* 

>> 

0 

4- 

44 

0 

£ 

© 

c 

(0 

4* 

tJ 

© 

4> 

O 

0 

GL 

C 

vH 

c 

& 

© 

c 

£ 

0 

cd 

A 

& 

4- 

a 

' 

d 

v 

13 

s 

r 

d 

00 

0 

«H 

0 

(4 

c 

0 

0 

1 

4- 

u 

H 

JS 

© 

G 

•H 

d 

4- 

*0 

© 

4» 

d 

1/5 

0 

u 

© 

c 

u 

c 

© 

© 

a 

ob 

*0 

© 

c 

*0 

cd 

© 

CM 

0 

r -» 

c 

(0 

a 

d 

t) 

4* 

o 

1/5 

co 

o 

^4 

-t 

00 

r* 

^4 

r- 

o 

CM 

<0 

00 

CO 

* 

o 

05 

O 

o 

\0 

▼H 

^4 

CM 

o 

05 

nO 

iH 

'O 

05 

o 

O 

05 

co 

05 

00 

1-t 

O 

p* 

r* 

CO 

o 

X 

O 

o 

05 

05 

^4 

X 

rH 

▼H 

4t 

o 

'f 

4T 

NO 

1/5 

NO 

o 

NO 

t- 

05 

X 

r"4 

tH 

CM 

10  ® 


o  >> 


05 


co  ^ 

4 

o 

ti 

o 

NO 

f* 

p* 

05 

©  l-> 

4-H 

V 

CM 

U  <0 

U 

a  m 

q 

X  tJ 

Of 

o 

>> 

o 

NO 

X 

CM 

CO 

nO 

© 

• 

X 

d 

o 

© 

1/5 

1/5 

P* 

CM 

© 

i 

*0 

cm 

0)  JC 

tJ 

4i 

©  4» 

3 

01 

G 

4* 

01 

o 

>> 

o 

o 

X 

1/5 

* 

d  44 

10 

4 

NO 

d 

o 

X 

nO 

05 

X 

CO 

3  0 

cd 

t5 

^4 

CM 

0* 

1 

(Ti  © 

ti 

>1 

U 

•H 

u 

o 

>» 

o 

1/5 

o 

CM 

TT 

X 

G  d 

d 

d 

o 

cm 

05 

* 

U5 

X 

ti  3 

•H 

§ 

*0 

rH 

CM 

©  C 

e 

•4 

3  cc 

H 

aj 

© 

0 

• 

4* 

© 

•H 

© 

•*4 

▼H 

c 

© 

r-> 

CO 

(30 

© 

a 

d 

0 

d 

•H 

c 

>> 

c- 

fc 

u 

© 

0 

© 

0 

4> 

© 

« 

« 

C 

© 

«-> 

a 

0 

c 

4* 

> 

•H 

a 

E 

c 

© 

© 

«< 

•j 

w 

d 

0 

© 

0 

w 

H 

o 

o 

. 


■ 


Ds3 


67 


(tuo)  q;6uei  wajs  AjEiuud 


Figure  7.  Growth  curves  for  the  D.  sandwicense  dialiel  lines  inoculated  with 
Rhizobium  strain  B. 
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if  the  reduction  of  acetylene  by  was  linear.  The 
acetylene  reduction  diagrams  for  four  lines  are  presented  in 
Figure  8.  All  samples  indicate  linear  reduction  and  a  sample 
taken  at  one  hour  should  adequately  represent  the  amount  of 
acetylene  reduced  in  a  one  hour  period.  It  should  be 
cautioned  that  Figure  8  is  intended  to  refer  to  the  linear 
rate  of  acetylene  reduction  only  and  any  comparison  among 
genotypes  should  be  avoided.  The  acetylene  reduction  curves 
for  the  five  diatlel  lines  are  presented  in  Figure  9.  These 
lines  also  present  a  linear  reduction  rate. 

The  above  results  suggested  that  symbiotic  nitrogen 
fixation  resulted  from  nodulated  roots;  however,  the  studies 
did  not  indicate  whether  the  nodule  was  the  site  of  nitrogen 
fixation.  Nodules  were  removed  from  plant  roots  and  placed 
separately  in  one  vial  and  the  roots  in  another.  The  pairs 
of  vials  were  prepared  and  analysed  by  the  acetylene 
reduction  assay.  Only  the  nodule  sample  indicated  any 
acetylene  reduction.  It  may  be  concluded  that  the  nodules 
are  the  site  for  nitrogen  fixation  in  E££modjLufii. 

4.  Bottle  jar  design. 

During  Experiment  II  in  the  quantitative  genetic 
variation  section  (Level  3,  Table  1),  two  short-comings  of 
the  bottle  jar  design  became  apparent.  Firstly,  the  plant 
roots  would  grow  into  the  wick  and  it  was  dificult  to 
separate  the  nodulated  roots  intact.  It  was  also  observed 
that  the  environmental  variation  in  the  quantitative  genetic 


« 


. 


: 


Acetylene  Reduced  (nM/pl/hr) 


69 


Figure  8.  Acetylene  reduction  curves  for  several  conspicuously  different 
Desmodium  species  inoculated  with  Rhizobium  strain  A. 
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studies  was  excessively  large  and  appeared  to  be  confined 
largely  to  within  blocks  which  suggested  that  the  bottle  jar 
design  might  be  manipulated  in  an  effort  to  reduce  this 
variation.  It  was  suspected  that  the  sand  particle  size  and 
the  distribution  of  the  particles  within  the  bottle  jar  may 
contribute  to  this  environmental  variation.  Within  bottle 

jar  variation  might  be  reduced  if  the  sand  particle  size  was 
more  constant. 

The  analysis  of  variance  for  sand  sizes  is  presented  in 
Table  12.  Significant  differences  were  recorded  for  sand 
size  means.  The  ceil  means  and  within-cell  differences 
indicated  that  some  means  and  differences  had  larger 
variances  than  others  (Table  13).  The  <  4.0  mm  sand  size 
means  were  found  to  be  less  for  dry  weight  and  shoot  length 
while  the  <  1.0  mm  size  had  the  lowest  difference  and 
variance  for  the  within— cell  difference  between  the  two 
plants.  Sand  sizes  with  large  inter-particle  spaces  appeared 
to  have  larger  within  cell  differences  suggesting  that 
perhaps  small  sand  grains  assist  capillary  movement  of  the 
nutrient  solution. 

The  germination  of  the  plants  in  the  above  experiment 
was  of  considerable  interest.  It  was  noted  that  the  <  1.0  mm 
size  was  the  quickest  to  germinate  while  the  1.0  -  3.0  tntn 
size  emerged  2—3  days  later.  The  other  two  sizes  were 
intermediate  with  the  <  4.0  mm  size  being  slightly  earlier. 
The  sand  surface  of  the  1.0  -  3.0  mm  and  1.0  -1.5  mm  sizes 
dried  before  the  rock  mulch  could  be  applied  while  the  plant 
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Table  12.  Mean  squares  from  the  analysis  of  variance  of  the 
sand  size  experiment. 


Source 

df 

Dry  weight 
(  mg/ pi  ) 

Shoot  length 
( cm/pl  ) 

Blocks 

5 

455.80 

2.57 

Sand  size 

3 

1718. 92# 

35. 15*# 

Interaction 

15 

1191.83 

16 • 99 * 

Erro  r 

23 

568.70 

5.62 

*signlf leant  at 

the 

5%  level  of  probability. 

^^significant  at 

the 

1%  level  of  probability. 

Table  13.  Means  and  standard  deviations  and  within-cell 

differences  and  standard  deviations  for  the  sand 
size  experiment. 


Sand 

Size 
(  mm  ) 

Dry  wt. 

( mg/ pi) 

Shoot  length 
( cm/pl  ) 

Mean  Wcd+ 

Mean  Wcd+ 

< 

4.0 

42.4±1 1.3 

15 

1. 

0-3.0 

63 • 4±41 .7 

37 

1. 

o 

l 

• 

a\ 

67. 7±26. 5 

26 

< 

1.0 

62. 4± 13.7 

7 

2±  9.4 

1 1 .9± 1.3 

1 . 5± 1.1 

2±43. 3 

15 • 0±4 • 4 

3 • 6±3. 9 

0±  9.5 

15.5±3. 0 

2. 8± 1. 1 

2±  4.9 

14 • 2± 1 .9 

0.  9±0. 5 

+wl thin-cell  difference. 
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roots  of  the  1.0  -  3.0  mm  size  were  found  to  be 
predominantly  confined  to  the  wick  area.  This  information 
suggested  that  seeds  planted  in  the  finer  sand  sizes  imbibed 
water  more  readily.  These  results  may  be  due  to  a  uniform 

capillary  movement  of  the  nutrient  solution  within  the  <  1.0 
mm  sand  size. 

The  reduction  of  the  within  bottle  Jar  variation  would 
reduce  environmental  variation?  therefore  the  modified 
bottle  jar  was  accepted  as  a  less  variable  container. 

5.  Analysis  of  variance  of  population  means  for  individual 
epistatic  effects. 

Voigt  ei  &l.  (1966)  stated  that  an  analysis  of 
population  means  from  a  cross  between  two  inbred  lines 
"explains  as  much  of  the  variance  among  means  as  possible  in 
terms  of  additive  gene  action”.  Meredith  and  Laster  (1975) 
report  that  "the  estimates  of  A  and  D  are  independent  and 
uncorrelated  only  if  E  is  not  significant".  To  investigate 
these  seemingly  contradictory  statements  a  coded  matrix  for 
the  genetic  parameters  (Mather  and  Jinks,  1971)  was 
constructed  and  three  types  of  least  squares  analysis  of 
this  matrix  were  Investigated.  The  three  types  of  analyses, 
differing  in  order  of  parameter  input,  were  (a)  the  foward 
solution  least  squares  analysis  which  assumed  a  simple  model 
and  which  explained  the  maximum  amount  of  the  epistatic 
effects  in  terms  of  additive  gene  action,  (b)  a  restricted 
stepwise  least  squares  analysis  which  assumed  a  simple  model 
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by  forcing  into  the  analysis  the  additive  and  dominance 
effect  vectors  and  then  explained  the  remaining  variance 
among  the  population  means  by  entering  the  coded  interaction 
vector  with  the  largest  partial  correlation.  The  remaining 
vectors  were  entered  by  the  same  criterion  until  the 
variance  of  the  population  means  was  explained.  And  (c)  a 
free  stepwise  least  squares  analysis  where  the  coded 
parameter  vectors  were  entered  in  terms  of  the  largest 
partial  correlation.  This  analysis  assumed  no  model. 

The  three  analyses  of  variance  of  the  population  means 
are  presented  in  Table  14.  For  this  study  the  restricted 
stepwise  least  squares  analysis  has  been  assumed  to  be  the 
superior  method  due  to  the  free  partition  of  the  epistatic 
effects.  If  the  epistatic  effects  were  not  partitioned,  the 
restricted  stepwise  least  squares  analysis  was  identical  to 
the  forward  solution  least  squares  analysis  for  order  of 
parameter  input. 

According  to  Meredith  and  Laster  (1975)  the 
significance  of  the  epistatic  effects  caused  the  loss  of 
independence  between  the  additive  and  dominance  effects.  To 
present  the  relationships  between  the  genetic  parameters  and 
the  population  means,  path  diagrams  (Li,  1975)  were 
constructed  for  the  above  analysis.  The  path  diagrams  are 
presented  in  Figure  10.  The  additive  and  dominance  effects 
remain  orthogonal  for  both  the  direct  and  the  indirect 
pathways.  The  non-orthogonality  was  between  the  main  effects 
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coefficients  of  determination  (  r2 )  between  the  main  effects 
and  the  interaction  effects  the  analysis  of  variance  of 
population  means  should  be  Interpreted  with  caution,  since  a 
large  proportion  of  the  main  effects  may  be  explained  as 
interaction  effects.  However,  the  analysis  is  of  some  value 
because  the  main  effects  were  orthogonal  even  in  the 
presence  of  significant  eplstatlc  effects  and  also  provided 

an  estimate  of  the  significance  of  the  Interaction 
parameters. 


6.  Strain-specific  ineffective  nodulation  interaction. 

Several  lines  of  conspicuously  different  fieamsdiars 
species  were  tested  with  the  selected  strains  to  determine 
the  type  of  host-strain  interaction.  The  interactions  are 
presented  in  Table  15.  D.  &2a£irje.s  lines  Dal,  Da2  and  Da3 
were  found  to  be  ineffectively  nodulated  by  strain  B.  These 
lines  also  exhibited  small  seed  size  and  all  originated  from 
Costa  Rica  or  Guatemala  (Imrie,  1973).  The  D.  aperjpeg  lines 
Da4  and  Da5  exhibited  effective  nodulation  interactions  with 
strain  B  as  well  as  large  seeds.  These  lines  were  originally 
collected  from  Colombia  and  Mexico,  respectively. 


E.  Qualitative  genetic  variation. 

The  area  of  study  of  the  inheritance  of  nitrogen 
fixation  which  has  commanded  the  greatest  interest  has  been 
the  study  of  non-nodulatlng  or  strain-specific  ineffective 
nodulation.  These  conditions  were  found  to  be  produced  by 
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Table  15,  Strain-specific  interactions  of  Pesmodl um 
with  EiliiotiiLlin  strains* 


Lines 

—r— 

Shizobi um 

- 

A 

Strain 

B 

Dsl 

E 

E 

Ds2 

E 

E 

Ds3 

E 

E 

Ds4 

E 

E 

Ds5 

E 

E 

Dal 

E 

I 

Da  2 

E 

I 

Da3 

E 

I 

Da4 

E 

E 

Da  5 

E 

E 

Dp 

E 

E 

Du 

E 

E 

Dsp 

E 

E* 

Dse 

I? 

E* 

E  ~  effective  nodulation  (green  plant)* 

~  ineffective  nodulation  (yellow-chlorotic  plant)* 
-  not  verified* 


1  ines 


appears  to  be  pH  sensitive* 


. 
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major  genes*  In  preliminary  experiments  of  this  study  in 
fiesisod ium  it  was  noted  that  a  strain-specific  ineffective 
nodula  t  ion  interaction  was  found  between  some  £.  ao^r  j_Re& 
lines  when  inoculated  with  strain  B*  The  numbers  of 
observations  for  each  of  the  scores  of  the  F2  plants  for 
nodula  t i on  interaction  and  spotted  leaf  character  are 
presented  in  Table  16*  The  statistical  analyses  of  the  F2 
plants  suggested  that  inheritance  of  this  strain-specific 
ineffective  nodulatlon  interaction  was  conditioned  by  a 
dominant  allele  at  a  single  locus* 

The  analysis  of  the  leaf  flecking  character  indicated 
the  presence  of  complementary  gene  interaction  in  the  Dsl  X 
Da2  cross.  The  Dsl  X  Dal  cross  appeared  to  be  segregating  at 
a  single  locus  as  the  data  conformed  to  a  3:1  segregation 
ratio*  The  joint  segregation  of  the  ineffective/effective 
nodulation  and  spotted/ pla in  leaf  character  for  the  Dsl  X  Da 
cross  suggested  that  the  two  characters  were  independent  in 
i nheri tance • 

test  the  hypothesis  that  a  single  dominant  gene 
conditioned  the  strain-specific  ineffective  nodulation 
interaction  between  some  D*  aflarlnes  lines  and  strain  B,  the 
backcrosses  to  Dsl  were  grown  and  the  plants  scored  for 
green  vs  yellow-chlorotic  color  and  for  spotted  vs  plain 
leaves*  The  character  scores  for  the  backcross  plants  are 
presented  in  Table  17*  These  data  for  nodulation  conformed 
to  the  expectations  of  a  single  gene  hypothesis  thereby 
suggesting  the  single  gene  hypothesis  for  the  strain- 
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Table  17,  Observed  nodulation  and  leaf  fleck  scores  and  Chi-square  probabilities  for 
the  backcross  plants  of  the  Iifii&n&dXiim  crosses  inoculated  with  strain  B« 
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specific  ineffective  modulation  in  £e£Ei£M2iLim.  These  results 
have  been  published  by  Pinchbeck  (1978)  (Appendix  D)  and  are 
in  accord  with  the  published  results  for  soybean  where  three 
single  dominant  genes  were  reported  for  strain-specific 
ineffective  interactions  (Caldwell,  1966;  Vest,  1970;  Vest 
and  Caldwell,  1972  ). 

Seventy  randomly  selected  ineffectively  nodulated  F2 
plants  were  inspected  for  nodulation  characteristics  and  re¬ 
planted  in  pots  and  supplied  with  nitrogen.  Ail 
ineffectively  nodulated  plants  responded  to  nitrogen 
fertilizer.  The  roots  of  the  ineffectively  nodulated  plants 
exhibited  various  degrees  of  nodule  formation.  The  majority 
of  nodules  were  small  whereas  a  few  were  near  normal  size. 
All  nodules  lacked  pink  pigmentation,  a  sign  of  an  active 
nodule  (Vincent,  1970),  These  nodules  appeared  to  be  similar 
to  those  reported  by  Caldwell  (1966)  and  Vest  (1970)  for  the 
i?J2  and  Rj3  genes  in  soybean. 

The  symbols  Nj/nx  are  suggested  for  the  designation  of 
the  !>•  aparl  nes  strain-specific  ineffective  gene.  D. 
gS-Xldwj.cec^e  has  the  genotype  nltm  while  the  ineffectively 
nodulated  D,  a  pari  neg,  lines  have  the  genotype  Nt  »  Ni  . 

Similar  phenotypes  for  ineffective  nodulation 
interaction  were  exhibited  by  Dal  and  Da2  lines  •  However, 
it  is  not  known  whether  these  phenotypes  were  conditioned  by 
similar  or  different  genes.  To  confirm  this  question  a 
standard  test  of  functional  allelism  (Srb,  Owen  and  Edgar, 
1965)  would  need  to  be  applied  to  the  F2  generation  of  a  Dal 
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X  Da2  cross* 

Experimentation  for  the  relationship  between  the  Nt 
gene  and  seed  size  was  not  pursued  in  this  study  as  a  single 
gene  influencing  seed  size  has  not  been  reported.  This 
avenue  would  be  best  examined  by  exploring  the  fertile  D. 
aparlnes  X  £.  unci  net  cross  since  both  species  exhibited 
the  extremes  for  seed  size  and  have  different  genotypes  for 
the  N*  gene* 

C*  Quantitative  genetic  variation. 

A  great  need  in  Interna t ion a 1  a griculture  has  been  to 
increase  total  symbiotic  nitrogen  fixation*  The  traditional 
method  of  increasing  total  fixation  has  been  manipulation  of 
the  EM^obium  strains  used  for  inoculum  (Vincent,  1974).  An 
alternate  research  avenue  is  to  establish  and  exploit  the 
variability  within  the  host.  It  was  an  objective  of  this 
study  to  investigate  the  hypothesis  that  there  is  variation 
within  Deg modi um  for  total  symbiotic  nitrogen  fixation  and 
that  this  variation  has  a  heritable  basis.  The  diallel  cross 
and  a  cross  between  two  inbred  lines  have  been  employed  as 
ouantitative  genetic  designs  to  establish  this  genetic 
variation. 

1.  Diallel  cross, 
a.  Experiment  I. 


The  analysis  of  variance*  for  two  characters  from 


. 


- 
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Experiment  I  Is  presented  in  Table  18.  Tukey's  test  for  non¬ 
additivity  of  the  data  (Snedecor  and  Cochran,  1967) 
indicated  that  a  transformation  of  the  data  was  necessary. 
Vincent  (1970)  suggested  that  dry  weight  accumulation  data 
should  be  analysed  in  the  logrithmetic  scale.  Non- 
significance  for  genotypes  was  recorded  for  dry  weight  and 
the  cause  appeared  to  be  block  2  which  was  significantly 
different  from  the  remainder  (Table  19).  This  block  may  have 
differed  for  three  reasons:  (a)  due  to  the  mechanical 
failure  of  the  gas  chromatograph  this  block  was  harvested 
one  week  later  than  the  remainder,  (b)  block  2  may  have  been 
inoculated  with  a  contaminated  broth  of  Rhizobiup.  since  the 
values  of  parent  4  and  some  of  its  crosses  are  greatly 
increased,  or  (c)  that  the  bottle  type  used  in  block  2 
created  a  different  environment  thereby  affecting  the  dry 
weight  accumulation  results.  Block  2  was  deleted  from 
further  analyses. 

Principal  components  analysis  of  the  collected 
characters  indicated  that  the  characters  were  highly 
correlated  and  related  to  one  principal  factor  (Table  20). 

The  genotype  means  for  the  characters  measured 
(transformed  to  logl0)  are  presented  in  Table  21.  The 
diallel  genotypes  revealed  that  most  crosses  exhibited 
heterosis  suggesting  that  heterosis  may  be  a  characteristic 


*Due  to  a  mechanical  failure  of  the  gas  chromatograph,  the 
acetylene  reduction  data  were  not  collected  for  blocks  2  and 
3  of  this  study.  Therefore  acetylene  reduction  character  was 
not  reported  for  Experiment  I. 
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Table  18*  Mean  squares  of  two  original  characters  for 
Experiment  I* 


Source 

d.  f  • 

Total 
shoo  t 
length 
( cm/pl ) 

Tota  l 

dry 

weight 

(  mg/ pi  ) 

Blocks 

3 

2084*# 

48750** 

Genotypes 

14 

109* 

2542 

Error 

42 

44 

1873 

Non— add i t ivi ty 

1 

6  10** 

23439** 

Res idual 

41 

31 

1347 

i gni f 1  cant  af  the  5%  level  of  probability# 
^significant  at  the  1%  level  of  probability# 
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Table  19.  Dry  weight  ( mg/pl  )  of  the  diallel  genotypes  for 
Experiment  I. 


Block 


Genotype 

1 

2 

3 

4 

Dsl 

70 

136 

51 

111 

Ds2 

41 

63 

31 

35 

Ds3 

91 

118 

43 

39 

Ds4 

11 

281 

51 

37 

Ds5 

45 

142 

49 

65 

Dsl  x  Ds2 

64 

109 

107 

111 

Dsl  x  Ds3 

115 

163 

74 

102 

Dsl  x  Ds4 

48 

316 

117 

120 

Dsl  x  Ds5 

60 

137 

81 

112 

Ds2  x  Ds3 

73 

194 

123 

54 

Ds2  x  Ds4 

42 

157 

34 

112 

Ds2  x  Ds5 

33 

214 

66 

122 

Ds3  x  Ds4 

52 

300 

65 

96 

Ds3  x  Ds5 

55 

220 

48 

141 

Ds4  x  Ds5 

40 

170 

51 

136 

Block  means 

56 

181 

66 

93 

-I 


Table  20.  Correlations  and  principal  eigenvector  loadings  for  Experiment 


87 


88 


o 

* 

+- 


c 

9 

(V 

! * 


V 

A 

to 

O 

C 

9 

u 

9 


a  a 

«J  >>  -4  M 
01  (.  11  E 
•J  tJ  *  w 


I 

J 

I 

I 

I 

I 

I 

i 

1 


05 

tn 

00 

co 

ao 

in 

05 

CO 

05 

00 

sO 

r* 

sO 

05 

t" 

05 

r» 

00 

cs 

CD 

CO 

O' 

00 

V© 

Oj 

tH 

t" 

CO 

CO 

co 

05 

CO 

*H 

05 

tH 

^4 

tH 

▼H 

tH 

v4 

o 

V 

i— 4 

A  f* 

*  a 

to 

«j 

0 

SB  N. 

X. 

v 

c 

C  E 

■o 

0 

A 

9  0 

H 

01 

*  w 

<H 

C 

' 

t, 

>> 

0 

(4 

u 

<d 

b 

t) 

£  p-> 

to 

c 

M 

a 

0 

0 

u  ^ 

*0 

0 

0 

C  E 

t* 

• 

9 

A 

9  0 

<d 

M 

W 

01 

*  w 

*0 

c 

-*■» 

\ 

a) 

c 

+- 

a 

(0 

5 

•H 

“0 

t 

£  f» 

c 

9 

+> 

v  a 

0 

a 

c 

0 

V. 

x 

*t4 

0 

e  6 

0) 

w 

<d 

A 

9  0 

c 

3» 

(Tj 

^  w 

d 

c. 

to 

0 

E 

<w 

I 

I 

I 


o 

* 

C" 


00 

o 

'<r 


o 

m 0 


co 

n 


00 

in 


in 

Gj 

CO 

10 

sD 

00 

uo 

05 

V 

tH 

05 

lO 

CO 

00 

tH 

05 

o 

00 

00 

00 

r~ 

o* 

43 

4; 

*4 

tH 

^4 

tH 

tH 

tH 

tH 

w4 

tH 

tH 

I 

I 

l 

l 

l 

> 


00 

o* 

00 

"0 

1/5 

r* 

05 

tH 

05 

o 

'/) 

05 

45 

vO 

O' 

CO 

05 

CO 

U5 

CO 

O 

CO 

C5 

00 

0* 

O' 

O' 

1/5 

co 

in 

43 

1/5 

00 

05 

43 

vC 

O' 

CO 

vO 

O' 

CO 

1/5 

1/5 

• 

o 

h 

l 


o 

O' 

05 

£ 

vO 

00 

© 

05 

o 

1/5 

05 

05 

rr 

05 

^4 

O 

05 

O 

00 

o 

^4 

i-4 

o 

4-1 

O 

tH 

o 

tH 

i-4 

^4 

tH 

co 

05 

in 

vC 

0* 

05 

1/5 

00 

o 

00 

1/5 

05 

O 

a 

o 

tH 

05 

o 

O 

tH 

e 

o 

▼H 

O 

1-4 

tH 

o 

o  co 

-  c 

M  «J 
0  9 
J  E 


05 

9 

r-4 

<Q 

0 

H 


& 

a 

>> 


o 


c 

05 

CO 

't 

9 

co 

05 

01 

(0 

O 

a 

Q 

a 

Q 

05 

co 

in 

CO 

CD 

01 

01 

(0 

(0 

Q 

O 

Q 

Q 

a 

X 

X 

X 

X 

X 

U5 

1-1 

H 

^4 

05 

m 

(0 

01 

co 

(0 

CO 

Q 

Q 

Q 

Q 

Q 

Q 

in 

Tf 

in 

in 

(0 

(0 

CO 

01 

01 

• 

o 

o 

Q 

C5 

o 

*0 

X 

X 

X 

X 

X 

<H 

0 

05 

05 

co 

CO 

to 

CO 

CD 

01 

01 

9 

Q 

Q 

Q 

Q 

Q 

(0 

l 


89 


c 

c  = 

• 

03 

0  l- 

M 

tafi 

•h  0 

sO 

CO 

so 

03 

Tf 

in 

03 

03 

00 

s© 

pH 

pH 

03 

't 

JZ 

0 

p  p 

© 

00 

o 

o 

S© 

00 

03 

0* 

pH 

CO 

O 

03 

03 

o 

00 

-H 

■H 

u 

S3  0 

co 

o 

Tf 

r- 

o 

O' 

in 

't 

M 

pH 

o 

03 

o 

in 

c 

p 

X  81 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

<H 

03 

*H 

*H  p 

o 

o 

o 

pH 

o 

o 

o 

o 

o 

o 

O 

o 

o 

O 

© 

o 

0 

s 

z 

p  r 

1 

1 

1 

1 

1 

•H 

u 

u 

0 

a. 

u 

a 

u 

X 

Q> 

u 

tJ 

u 

t. 

0 

o 

<w 

p  r-> 

tJ 

p 

4:  a 

NO 

a» 

00 

03 

03 

o 

o 

CO 

03 

in 

in 

00 

s0 

00 

c 

US 

c 

tf 

sO 

03 

co 

pH 

sc 

00 

pH 

03 

CO 

o 

CO 

in 

pH 

S3 

c 

p 

>s  p  &£ 

00 

in 

r- 

"t 

O' 

03 

03 

03 

03 

00 

O' 

00 

00 

00 

00 

pH 

-P 

<3 

0 

U  03  S 

• 

(A 

p 

H 

T)  2  w 

pH 

pH 

pH 

P< 

pH 

pH 

pH 

pH 

pH 

pH 

pH 

pH 

pH 

o 

t! 

C 

0 

to 

* 

P  V 

x 

c 

c 

0 

0) 

<D 

03 

a 

* 

P  P 

■P 

P 

43  a 

03 

o  - 

00 

sO 

© 

V© 

so 

03 

o 

s© 

03 

0* 

in 

in 

o 

03 

0 

tjf  V. 

V© 

pH 

00 

Tj* 

O' 

03 

CO 

03 

o 

CO 

03 

sO 

o 

H 

43 

© 

^5  *»3  01) 

03 

o 

o 

00 

o 

co 

03 

03 

03 

03 

pH 

03 

pH 

pH 

pH 

»-< 

tt 

£ 

U  03  s 

0 

03 

W 

*0  *  w 

fp 

pH 

rH 

o 

pH 

pH 

pH 

pH 

pH 

H 

pH 

pH 

p* 

pH 

pH 

o 

P 

t> 

C 

• 

03 

U 

•P 

03 

c 

0 

0 

•H 

w 

V 

03 

CO 

Tf 

in 

co 

Tf 

in 

o- 

in 

in 

pH 

03 

(A 

(A 

x 

X 

(A 

(A 

(A 

(A 

(A 

(A 

• 

03 

a 

>» 

p 

p 

P 

p 

p 

P 

p 

P 

P 

P 

P 

V 

03 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

0 

0 

43 

c 

PH 

03 

CO 

in 

pH 

pH 

03 

03 

03 

CO 

CO 

<3 

0) 

03 

US 

(0 

(A 

U3 

us 

<A 

X 

X 

(A 

(A 

X 

(A 

(A 

(A 

© 

H 

o 

P 

P 

p 

p 

P 

p 

P 

p 

P 

P 

P 

P 

P 

P 

P 

US 

90 


In  Desnodium  for  total  symbiotic  nitrogen  fixed. 

The  diallel  analyses  of  variance  for  the  loglo 
transformed  characters  in  Experiment  I  are  presented  in 
Table  22.  All  characters  including  the  "Nitrogen  fixation" 
factor  were  noted  to  have  a  similar  analysis  pattern.  The 
significant  mean  squares  for  genotypes  suggested  that 
genetic  variation  for  the  characters  existed  within  the  D. 
^anrfwlce.nse  population.  Significant  general  effects 
indicated  that  a  portion  of  this  variation  was  heritable. 
Although  heterosis  was  not  significant,  average  heterosis 
*as  shown  to  be  highly  significant  indicating  that  the  Fi 
progeny  means  were  significantly  different  from  the  parental 
means • 

The  estimates  of  the  G.C.A.  effects  and  variety 
constants  are  presented  in  Table  23.  Line  Dsl  exhibited 
positive  G.C.A.  effects  and  variety  constants  for  total  dry 
weight  while  line  Ds4  was  found  to  have  a  negative  G.C.A. 
effect  and  variety  constant. 

Estimates  of  the  heterosis  constants  revealed  that,  in 
general,  the  line  and  specific  heterosis  constants  were 
small  in  magnitude  in  comparison  to  the  average  heterosis 
constant  (Table  24).  These  results  were  expected  since  line 
and  specific  heterosis  effects  were  not  significant  in  the 
diallel  analysis  of  variance. 

The  estimates  of  the  S.C.A.  effects  failed  to  define 
»ny  hybrid  which  was  superior  to  the  remainder  (Table  25). 


The  negative  S.C.A. 


estimates  for  inbred  lines  suggested 


■ 


a*  *tntx4utn  -  r  ?  f  r  '"  v  *l<  ^  .!*-*•**» 


Table  22.  Mean  squares  from  the  dlallel  analysis  of  variance  for  Experiment 
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Table  25.  Estimates  of  Griffing's  S.C.A.  effects  for 
Experiment  I. 


S.C.A. 

effect 

Ma  in 

shoo  t 
length 
( cm/pl ) 

Sec  onda  ry 
shoot 
length 
( cm/ pi ) 

Total 

shoot 
length 
( cm/pl ) 

Leaf 
dry 
weight 
(  mg/ pi  ) 

s(  Ds  1  ) 

-.  059 

-.311 

-.090 

-.136 

s(  Ds2  ) 

-.086 

-.521 

-.140 

-.171 

s(  Ds3  ) 

-.064 

-.393 

-.094 

-.127 

s(  Ds4  ) 

-.109 

-.525 

-.  139 

-.139 

s(  Ds5  ) 

-.038 

-.247 

-.081 

-.092 

se[  s(  i  )  ] 

.  041 

.  188 

.051 

.069 

se[  s (  i  )-s(  j  )  ] 

.055 

.252 

.  068 

.093 

s(  Ds  1  f  Ds2  ) 

.  076 

.  146 

.095 

.097 

s(  Ds  1  ,  Ds3  ) 

-.  004 

.  195 

.023 

.047 

s(  Ds  1 1  Ds4  ) 

.068 

.283 

.  071 

.113 

s(  Dsl,Ds5) 

-.022 

-.002 

-.009 

.015 

s(  Ds2 1  Ds3  ) 

.048 

.  331 

.  081 

.103 

s(  Ds2  r Ds4  ) 

.021 

.390 

.055 

.087 

s(  Ds2  y  Ds5  ) 

.028 

.176 

.050 

.054 

s(  Ds3  T  Ds4  ) 

.  072 

.  158 

.  058 

.086 

s(  Ds3 y  Ds5  ) 

.013 

.  102 

.  026 

.018 

s(  Ds4, Ds5  ) 

.  057 

.218 

.  095 

.097 

se[  s(  i,  j)] 

.  052 

.238 

.  064 

.088 

se[s(  if  j)-s(  i  »  k  )  ] 

.  078 

.356 

.  096 

.132 

se[  s(  i,  j  )-s(  k,l)] 

.  071 

.325 

.  088 

.120 

.1  ►«-*" 


•  - 

— 
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Table  25 


(coni:).  Estimates  of  Grlfflng's  S.C.A.  effects  for 
Experiment  I. 


S  •  C  •  A  • 

effect 

Shoo  t 
dry 
we Ight 
( mg/ pi  ) 

Total 
dry 
weight 
( mg/pl ) 

Ni trogen 
fixation 
"factor" 

s( Dsl  ) 

132 

-.  140 

-.549 

s(  Ds2  ) 

-.167 

-.162 

-.747 

s(  Ds3  ) 

-.  152 

-.134 

-.  580 

s(  Ds4  ) 

-.192 

-.190 

-.840 

s(  Ds5  ) 

-.065 

-.088 

-.374 

se[  s(  i  )] 

.  079 

.072 

.  284 

se[  s (  i  )-s(  J  )  ] 

.106 

.  097 

.  381 

s(  Dsl,  Ds2  ) 

.113 

.  102 

.446 

s(  Ds  1 ,  Ds 3  ) 

.044 

.046 

.  178 

s(  Dsl  ,  Ds4  ) 

.  154 

.126 

.  518 

s( Dsl , Ds5  ) 

-.  047 

•  006 

-.044 

s( Ds2, Ds3  ) 

.115 

.  108 

.  468 

s( Ds2,Ds4 ) 

.  066 

.065 

.  346 

s(  Ds2,Ds5  ) 

.  039 

.  049 

.234 

s(  Ds3, Ds4  ) 

.102 

.092 

.386 

s(  Ds3  ,  Ds5  ) 

.  043 

.023 

.  127 

s(  Ds4, Ds5  ) 

.  096 

.  097 

.  431 

se[  s(  i  ,  j  )  ] 

.100 

.091 

.359 

se[  s(  i,  j)-s(  i  ,  k  )  ] 

.  149 

.  137 

.  538 

se[  s(  i  ,  J  )-s(  k  ,  1  )  ] 

.  136 

.125 

.  491 
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that  homozygous  lines  may  be  exhibiting  inbreeding 
depression  for  total  symbiotic  nitrogen  fixed. 


b.  Experiment  II. 


The 


data  from  Experiment  I  indicated  that  total 
symbiotic  nitrogen  fixed  appeared  to  be  heritable.  To 
further  investigate  the  hypothesis  that  there  exists  genetic 
variation  for  nitrogen  fixation,  the  diallet  cross  was  grown 
with  £h_Uobiym  strain  A  to  determine  if  the  mode  of 
inheritance  was  consistent  for  a  different  strain. 

Principal  components  analysis  of  the  characters 
measured  for  Experiment  II  revealed  that  all  characters  were 
highly  correlated  and  related  to  one  factor  (Table  26). 

The  genotype  means  from  the  log10  transformed 
characters  are  presented  in  Table  27.  An  examination  of  the 
individual  genotype  means  revealed  that  every  cross 
exhibited  he t erobe 1 1 i os  is  except  cross  Ds2  x  Ds3  which 
portrayed  only  complete  heterosis  for  the  acetylene 
reduction  assay. 

The  diallel  analyses  of  variance  for  the  log1Q 
transformed  characters  for  Experiment  II  are  given  in  Table 
28.  There  were  significant  differences  among  genotypes  which 
confirmed  that  genetic  variation  was  present  in  the  diallel 
population.  General  effects,  heterosis  and  average  heterosis 
components  were  significant  while  line  and  specific 
heterosis  were  non-significant.  The  genetic  variance  was 
confined  to  among  parental  lines  and  between  the  progeny  and 


Table  26.  Correlations  coefficients  and  principal  eigenvector 
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Table  27.  Loe10  means  and  standard  error  of  the  difference 
between  two  means  for  Experiment  II. 


Genotype 

Leaf 
dry 
we ight 
( mg/pl ) 

Shoo  t 
dry 

weight 

( mg/pl ) 

Total 

dry 

weight 

( mg/pl  ) 

Dsl 

1*627 

1.056 

1.732 

Ds2 

1 . 660 

1.  147 

1.777 

Ds3 

1.838 

1.355 

1.962 

Ds4 

1.749 

1.127 

1.842 

Ds5 

1.795 

1.272 

1.909 

Dsl  x  Ds2 

1.842 

1. 308 

1.954 

Dsl  x  Ds3 

1.948 

1.452 

2.  070 

Dsl  x  Ds4 

2.  003 

1.466 

2.115 

Dsl  x  Ds5 

1.950 

1.457 

2.073 

Ds2  x  Ds3 

1.914 

1.425 

2.037 

Ds2  x  Ds4 

2.  035 

1.494 

2.  145 

Ds2  x  Ds5 

2.036 

1.525 

2.153 

Ds3  x  Ds4 

2.214 

1.683 

2.327 

Ds3  x  Ds5 

2.  141 

1 . 640 

2.  261 

Ds4  x  Ds5 

1.  941 

1.346 

2.041 

se  of  d. 

0.164 

0.  199 

0.172 

e^r .  r 


t 
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Table 

27 

( cont ) • 
di f  fere 

Log10  means  and  standard 
nee  between  two  means  for 

error  of 
Experl  me 

Genotype 

Ac  etylene 
reduced 

( nmoles/ pi /hr ) 

Ni trogen 
f Ixa ti on 
*'  f  ac  tor" 

Dsl 

1.648 

1.831 

Ds  2 

1.488 

1. 831 

Bs3 

1.864 

1  •  847 

Ds4 

1.50  6 

1.835 

Ds5 

1.726 

1.842 

Ds  1  x 

Ds2 

1.759 

1. 845 

Dsl  x 

Ds3 

1 . 880 

1.853 

Ds  1  x 

Ds4 

1.835 

1.855 

Ds  1  x 

Ds5 

1.996 

1.855 

Ds2  x 

Ds3 

1 . 86  3 

1.851 

Ds2  x 

Ds4 

1.898 

1 . 857 

Ds2  x 

Ds5 

2.  122 

1. 861 

Ds3  x 

Ds  4 

2.  083 

1.87  0 

Ds3  x 

Ds5 

2.  030 

1.866 

Ds4  x 

Ds5 

1.732 

1.849 

se  of 

d  • 

0.  206 

0.460 
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parental  means.  Block  x  genotype  interactions  were  not 
significant  for  dry  weight  characters  but  significant  for 
the  acetylene  reduction  assay.  This  component  was  distorted 
because  of  a  problem  of  the  separation  of  the  nodulated 
roots  from  the  lamp  wick  (Figure  2)  in  the  preparation  of 
the  plants  for  the  assay.  As  a  result  of  this  problem  the 
bottle  jar  design  was  modified  (Figure  3). 

The  estimates  of  the  G.C.A.  effects  and  variety 
constants  are  presented  in  Table  29  and  line  Ds3  was  shown 
to  have  the  greatest  G.C.A.  effect  and  variety  constant  for 
all  characters-  Line  Dsl ,  however,  exhibited  a  negative 
G-C-A,  effect  and  variety  constant  which  is  the  reverse  of 
Experiment  I.  These  results  were  of  interest  since  line  Ds3 
was  a  Dsl  backcross  line  from  an  original  Dsl  x  D.  a oaring^ 
(Dal)  cross*  Since  strain  A  was  isolated  from  D.  aag.ill.neg,  a 
strain  x  genotype  interaction  may  be  suggested* 

The  estimates  of  the  heterosis  constants  again  revealed 
that  average  heterosis  was  much  larger  than  either  the 

estimates  of  the  line  or  specific  heterosis  constants  (Table 

30) . 

The  S • C  *  A •  effect  estimates  revealed  inbreeding 
depression  in  the  Inbred  lines  for  all  characters  (Table 

31) .  The  cross  Ds3  x  Ds4  showed  considerable  S.C.A.  effect 
for  the  dry  weight  characters.  However,  the  cross  Ds2  x  Ds5 
exhibited  a  superior  S.C.A.  effect  for  the  acetylene 
reduction  assay.  The  cross  Ds2  x  Ds5  also  had  a  larger 
acetylene  reduction  assay  S.C.A.  value  relative  to  its  dry 
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weight  accumulation  S.C.A.  value,  whereas,  the  cross  Ds3  x 
Ds4  showed  equal  dry  weight  accumulation  and  acetylene 
reduction  S.C.A.  values.  These  results  suggested  that  dry 
weight  accumulation  and  the  acetylene  reduction  assay, 
although  assumed  to  measure  symbiotic  nitrogen  fixation,  may 
measure  closely  related  components  of  a  complex  character. 

c.  Experiment  III, 

Since  both  strain  and/or  environmental  effects  appeared 
to  be  of  importance  in  the  preceding  studies,  Experiment  III 
was  inititated  to  determine  if  strain  X  genotype 
interactions  were  of  significant  importance  in  a  D. 
sajidw  1  Qgn£e  population, 

A  review  of  the  original  and  l og j q  transformed  means 
for  the  within-s tral n  and  overall  analyses  for  Experiment 
III  illustrated  that  heterosis  was  of  the  approximate 
magnitude  of  1.6  (Tables  32,  33  and  34;  respectively).  The 
low  mean  values  for  the  acetylene  reduction  analysis  in 
strain  A  (Table  32)  appeared  to  be  due  to  the  nodules  of 
that  strain  being  in  a  state  of  senescence.  Nodule 
senescence  may  have  been  due  to  the  harsh  "environmental" 
conditions  imposed  upon  this  experiment;  yet  strain  B 
appeared  unaffected.  Since  strain  A  was  faster  growing  than 
strain  B  the  nodules  of  strain  A  may  have  been  sampled  at  a 
later  physiological  stage.  The  compatability  of  the  dry 
weight  means  may  be  explained  by  the  releasing  of  nitrogen 


from 


the 


decaying  nodules  for  plant  growth. 


Further 


. 


Table  32.  Means  and  standard  errors  of  the  difference  between  two 
means  for  the  &hlzpfc> j,mp  strain  A  inoculated  blocks  for 
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Table  34.  Means  and  standard  error  of  the  difference  between  two 
means  from  the  overall  analysis  for  Experiment  III. 
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experimentation  would  be  necessary  to  clarify  these  results. 

The  within-st rain  diallel  analyses  of  variance  for  the 
normal  characters  are  presented  in  Table  35  and  for 
comparison  the  log10  transformed  diallel  analyses  of 
variance  is  given  in  Table  36.  The  similarity  of  the  two 
scales  suggested  that  scale  effects  had  been  removed  from 

this  experiment  through  improved  experimental  and  bottle  jar 
des ign. 

Significant  variation  was  again  found  among  the 
genotypes  and  the  majority  of  this  variation  was  due  to 
average  heterosis.  Block  X  genotype  interactions  were  also 
found  to  be  significant  for  must  genetic  components 
suggesting  genetic  component  X  environment  interactions. 
These  results  were  expected  since  Vincent  (1965)  reported 
that  the  environment  Influenced  symbiotic  nitrogen  fixation. 

Factorial  diallel  analyses  of  variance  for  original  and 
1-O@io  transformed  characters  are  presented  in  Table  37. 
Strains  were  not  significantly  different  for  dry  weight 
accumulation  yet  highly  significant  for  the  acetylene 
reduction  assay  for  both  scales  due  to  nodule  senescence  in 
Strain  A.  Genotypes  were  highly  significant  for  all 
characters  and  scales  suggesting  that  genetic  variation  was 
present  in  this  population  and  that  this  variation  was  scale 
independent • 


A  partition  of  the  genotype  variance  suggested  that 


a 


predominant  portion  of  this  variance  was  due  to  dominance 
deviations  since  general  effects  were  significant  for  the 
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logl°  transformed  acetylene  reduction  assay  only.  The 
additive  effects  which  were  present  in  the  wi thi n-stra i n 
analyses  may  be  cancelled  in  a  factorial  analysis  (Comstock, 
fi.E.,  pers.  comm.).  Average  heterosis  was  highly  significant 
for  all  analyses  while  the  significance  of  line  and  specific 
heterosis  was  character  and  scale  dependent. 

Strain  by  genotype  interactions  were  found  to  be 
significant  in  this  experiment.  Strain  X  heterosis  and 
strain  X  specific  heterosis  Interactions  were  significant 
for  all  analyses  while  there  appeared  to  be  little  strain  X 
average  heterosis  interaction. 

A  limitation  to  this  analysis  is  the  limited  number  of 
degrees  of  freedom  for  the  interaction  terms.  Since  strains 
have  been  considered  to  be  an  environmental  component  in 
this  study,  the  strains  and  blocks  were  confounded  and 
analysed  as  a  four  block  experiment  so  that  direct 
comparisons  could  be  made  with  Experiment  II  (Table  28).  The 
overall  dialiel  analyses  of  variance  for  the  original  and 
lo§io  transformed  characters  are  presented  In  Table  38. 
Genotypes  were  found  to  be  significantly  different  for  both 
scales  as  were  all  components  with  the  exception  of  specific 
heterosis  for  the  logio  transformed  acetylene  reduction 
assay.  All  interaction  terms  were  significant  with  the 


exception  of  blocks/strai ns  X  general  which  wai 


non¬ 


significant  for  the  logt o  transformed  dry  weight 
accumulation.  The  significance  of  the  interaction  terms 
suggested  that  genetic  studies  for  total  symbiotic  nitrogen 
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fixed  should  be  tested  with  different  strains  and  in 

different  environments  if  the  experimental  results  are  to  be 
re  liable* 

The  error  mean  squares  for  dry  weight  accumulation  and 
the  acetylene  reduction  assay  in  Experiment  III  (Table  38) 
have  been  reduced  by  approximately  one-half  those  of 
Experiment  II  (Table  28).  This  reduction  is  most  likely  due 
to  the  modification  of  the  bottle  jar  design  (Figure  3) 
which  reduced  withln-cell  variation*  Once  the  variation  due 
to  bottle  jar  design  was  removed,  the  log10  transformation 
was  not  necessary  as  the  original  analyses  were  similar  to 
the  l l o  analyses*  These  data  suggest  that  further 
investigation  of  the  design  of  bottle  jars  or  sand  size 
would  be  beneficial* 

The  estimates  of  the  variety  constants  and  G*C*A* 
effects  for  the  original  and  logio  transformed  characters 
from  all  analyses  are  presented  in  Tables  39  and  40, 
respectively*  The  estimates  of  the  G.C.A*  effects  for  the 
logio  transformed  characters  for  strain  B  (Table  39)  were 
found  to  be  similar  to  the  results  of  Experiment  I*  Line  Dsl 
was  found  to  have  the  largest  G*C*A*  effect  for  both 
characters  while  Ds4  showed  the  largest  negative  G*C,A. 
effect*  Line  Ds2  showed  the  largest  heterosis  constant  for 
dry  weight  but  Dsl  had  the  larger  constant  for  the  acetylene 
reduction  assay*  The  G*C*A*  effects  for  strain  A  (Table  39) 
were  not  similar  to  those  reported  for  Experiment  II*  As 
earlier  discussed  the  nodules  of  the  plants  of  strain  A  were 
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Table  40,  Estimates  of  G.C.A.  effects  and  variety  constants  for  the  log.n 
transformed  characters  for  Experiment  III* 
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in  a  state  of  senescence  and  this  physiological  state  may 
have  affected  the  results  for  the  acetylene  reduction 
analysis.  However  the  G.C.A.  effects  for  dry  weight 
accumulation  were  similar,  therefore  nodule  senescence  may 
not  explain  the  divergence  from  the  G.C.A.  effects  of 
Experiment  II.  Line  Ds5  presented  the  largest  G.C.A.  effect 
for  strain  A.  Its  crosses  had  larger  means  than  the  other 
crosses  (Table  32).  Line  Ds5  and  its  crosses  may  have  been 
resistant  to  nodule  senescense  or  responded  more  quickly  to 
clipping  of  the  juvenile  plants.  Line  Ds5  also  has  the 

larges!  G.C.A.  effect  for  both  characters  for  the  overall 
analysis  (Table  39). 

The  estimates  of  the  heterosis  constants  and  the  S.C.A. 
effects  derived  from  the  original  and  the  log1Q  characters 
in  the  various  analyses  are  presented  in  Tables  41  to  44, 
respectively. 

d.  Combined  analysis. 

The  preceding  experiments  were  conducted  with  the  same 
diallel  cross  but  differed  in  environmental  and  strain 
effects.  Since  strains  have  been  initially  considered  as  an 
environmental  effect,  it  would  be  pertinent  to  pool  the 
above  experiments  for  a  combined  analysis.  Blocks  1,2  and  3 
of  Experiment  I  have  been  deleted  from  the  combined  analysis 
because  subsamples  within  cells  were  not  available  for  block 
1  and  the  acetylene  reduction  assay  was  not  recorded  for 
blocks  2  and  3.  The  remaining  nine  blocks,  six  inoculated 
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^1  lh  strain  A  and  thrcs  c  t  «•<>  >  „  n  , 

mree  with  strain  B,  were  pooled  such  that 

block  and  strain  effects  were  confounded* 

The  original  and  logl0  means  confirmed  that  the  mean  of 
the  F i  genotypes  was  much  larger  than  the  means  of  their 
inbred  parents  (Table  45).  The  mean  heterosis  for  all 
characters  appeared  to  be  approximately  1.8. 

The  di all el  analyses  of  variance  for  Gardner's  Analysis 
II  of  the  original  and  log10  transformed  characters  of  the 
combined  analysis  are  presented  in  Table  46  and  for 
comparison  Gardner's  Analysis  III  is  provided  in  Table  47. 
Both  analyses  are  compatible;  however,  the  general  mean 
square  was  larger  for  Analysis  III  than  for  Analysis  II 
suggesting  that  heritable  effects  are  given  more  weighting 
by  the  Analysis  III.  Gardner  and  Eberhart  (1966)  report  that 
Analysis  III  may  be  superior  when  average  heterosis  is 
highly  significant  as  in  the  present  data;  however,  either 
analysis  is  applicable. 

Significant  differences  among  genotypes  were  found  for 
all  characters  except  the  original  acetylene  reduction 
assay.  These  results  were  not  unexpected  since  the  acetylene 
reduction  assay  values  in  the  earlier  experiments  were 
affected  by  the  difficulty  of  removing  the  root  intact  from 
the  bottle  jar  wick.  Highly  significant  heterosis  and 
average  heterosis  effects  were  also  recorded.  All 
interaction  effects  for  the  logj  Q  transformed  characters 
were  significant  with  the  exception  of  block/ strains  x 


average 


heterosis  for  the  logio  transformed  acetylene 


"•  ■* 


I  ■  |  v  ,  I  ,<  ,-T 


Table  45.  Means  and  standard  error  of  the  difference  between  two 
means  for  the  combined  diallel  analysis. 
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reduction  assay.  These  results  suggested  that  genotype  x 
environment  interactions  were  of  significant  importance  in 
the  study  of  total  symbiotic  nitrogen  fixed.  Significant 
general  effects  for  the  log10  transformed  characters  for 
both  dial lei  analyses  suggested  that  heritable  genetic 
variation  was  present. 

The  combining  ability  analysis  for  G.C.A.  effects  and 
variety  constants  are  presented  in  Table  48.  As  indicated  by 
the  genotype  means,  line  Ds5  was  shown  to  have  the  largest 
G.C.A.  effect  and  variety  constant  for  the  loglo  transformed 
characters.  Lines  Ds2  and  Ds4,  on  the  other  hand,  had 
negative  G.C.A.  effects.  Heterosis  constants  and  S.C.A. 
effects  for  the  two  scales  are  presented  in  Tables  49  and 
50,  respectively.  The  average  heterosis  constant  was  much 
larger  that  the  other  heterosis  constants  suggesting  that 
the  major  portion  of  the  genetic  variance  is  due  to  Ft 
heterosis  (Table  49).  On  the  other  hand,  all  inbred  lines 
exhibited  negative  S.C.A.  effects  suggesting  inbreeding 
depression  (Table  50).  Further  studies  would  have  to  be 
undertaken  to  determine  whether  the  large  average  heterosis 
effects  were  due  to  Ft  heterosis  or  inbreeding  depression 
exibited  by  the  inbred  lines  when  grown  in  a  nitrogen 
limiting  environment. 

The  combined  analysis  indicated  that  total  symbiotic 
nitrogen  fixed  had  both  heritable  and  non-herl table 
components  of  the  genetic  variance.  The  larger  proportion  of 


var i ance 


was  contributed  by  the 


non— her i tab le  portion 


- 


Table  48*  Estimates  of  G# C*A,  effects  and  variety  constants  for  the 
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Table  50.  Estimates  of  Grlffing's  S.C.A.  effects  for  the 
combined  analysis* 


S  •  C  •  A 

effec t 

original 

logi  0 

Dry 

weight 

( mg/ d l  ) 

Acetylene 
reduced 
( nm/pl/hr  ) 

Dry 

weight 

( mg/ pi  ) 

Ac e  ty le ne 
reduced 
( nm/pl/hr  ) 

s(  Dsl  ) 

-68 

-  18.  3 

-.169 

-.172 

s<  Ds2  ) 

-20 

-18.3 

-.097 

-.140 

s(  Ds3  ) 

-70 

-21.4 

-.181 

-.  148 

s(  Ds4  ) 

-66 

-26.  3 

-.184 

-.170 

s( Ds5  ) 

-91 

-45.  1 

-.204 

-.219 

se[  s(  i  )  ] 

15 

12.2 

.032 

.  040 

se[  s(  i  )— s(  J  )  ] 

20 

16.3 

.043 

.0  54 

s(  Dsl  t  Ds2  ) 

9 

.5 

.059 

.045 

s(  Ds  1 1  Ds  3  ) 

71 

9.4 

.  135 

.  110 

s(  Dsl  t  Ds4  ) 

15 

4.6 

.055 

.053 

s( Dsl y  Ds5  ) 

42 

22.  2 

.  089 

.136 

s( Ds2r  Ds3  ) 

-16 

-2.4 

-.056 

-.  053 

s(  Ds2  y  Ds4  ) 

19 

1.5 

.090 

.125 

©(  Ds2, Ds5  ) 

27 

38.  0 

.100 

.163 

sC  Ds3  j Ds4  ) 

35 

26.2 

.  142 

.  131 

s(  Ds3  y  Ds5  ) 

51 

9.7 

•  139 

.109 

s(  Ds4t  Ds5  ) 

62 

20.  4 

.081 

.031 

se[  s( 1  ,  j  )  ] 

19 

15.4 

.041 

.051 

se[  s(  1 ,  J  )  —  s(  i  ,  k)  ] 

28 

23.1 

.061 

.  077 

se[  s(  1  ,  j  )-s(  k  ,  l  )  ] 

25 

21.1 

.056 

.  070 

-  - 

>  f 
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because  of  the  expression  of  average  heterosis.  Crosses  Ds4 
x  Ds5r  Dsl  x  Ds3  and/or  Ds2  x  Ds5  would  be  of  significant 
interest  for  an  intensive  study  since  these  crosses 
exhibited  high  mean  yield  and  specific  combining  ability 
effects  for  at  least  one  of  the  two  characters. 


e.  Correlation  between  dry  weight  accumulation  and  the 
acetylene  reduction  assay. 

A  secondary  objective  of  this  study  was  to  relate  dry 
weight  accumulation  and  the  acetylene  reduction  assay  as 
quantitative  measure  for  total  symbiotic  nitrogen  fixation. 
The  correlation  coefficients  for  these  characters  for  the 
diallel  experiments  are  presented  in  Table  51.  The 
correlation  coefficients  were  found  to  be  highly  significant 
with  the  exception  of  the  unadjusted  overall  correlation. 
These  data  suggest  that  the  acetylene  reduction  assay  is 
compatible  with  dry  weight  accumulation  for  juvenile  plants 
grown  in  bottle  jars.  A  similarly  large  correlation  (r  = 
.89)  was  found  by  Lawn  et  al.  (1974)  for  total  nodule 
activity  and  shoot  dry  weight  for  shoot  genotype  means. 

It  may  be  concluded  that  the  acetylene  reduction  assay 
would  be  suitable  as  a  plant  breeding  tool  for  screening  of 
juvenile  genotypes  into  various  classes  for  total  nitrogen 
fixed.  Since  the  acetylene  reduction  assay  does  not  require 
nitrogen-free  growth  conditions  this  assay  would  be 
effective  and  reliable  for  screening  plants  grown  under 


field  conditions. 


. 
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Table  51.  Unadjusted  and  block  adjusted  correlation 

coefficients  between  dry  weight  ( mg/pl )  and 
acetylene  reduced  (nM/pl/hr)  for  the  diallel 
exp  eriments • 


Exper imen t 

Unadjusted  ,,r” 

Bloc 
Adj  ust 

I 

•  92 

.86 

II 

.92 

. 

00 

a 

III 

.76 

•  83 

Combined 

.43 

.85 

*all  block  means  in  each  character  are  adjusted  to  the  grand 
mean  for  that  character* 


. 
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2,  Cross  between  two  inbred  lines. 

A  second  approach  to  the  study  of  quantitative  genetic 
variance  is  provided  by  the  analysis  of  populations  produced 
by  a  cross  between  two  inbred  lines.  Through  the  analysis  of 
first  degree  statistics,  the  cross  between  two  inbred  lines 
will  provide  information  about  epistatic  interactions  which 
cannot  be  obtained  from  a  diallel  cross.  In  this  study  of 
+he  inheritance  of  symbiotic  nitrogen  fixation  the  cross 
between  two  inbred  lines  has  been  employed  as  an  alternate 
quantitative  genetic  design  for  the  establishment  and  the 
partition  of  the  genetic  variance  and  to  investigate  the 
influence  of  eplstasis  in  the  expression  of  heterosis.  Three 
experiments  were  conducted  In  this  section. 

a.  Experiment  IV. 

Seven  blocks  were  grown  in  this  experiment,  three 
blocks  were  deleted  due  to  their  uniqueness.  Blocks  7  and  8 
were  grown  in  large  bottom  (21)  bottle  Jars  while  the 
remainder  had  small  bottles  (600ml).  Block  5  was  missing 
many  of  the  Bj  population  plants  and  was  not  considered  to 
be  complete.  Because  of  poor  germination  within  the  Bt 
population,  six  Bx  plants  from  the  deleted  blocks  were 
fitted  into  the  remaining  blocks  so  as  not  to  appreciably 
change  the  within  block  mean. 

The  original  and  logxo  means  revealed  Fx  heterosis  with 
the  F 2  mean  falling  very  near  the  midparent  value  suggesting 


the  random  segregation  of  alleles  (Table  52).  The  backcross 


1  «  -  > 


jr 


Table  52*  Means,  standard  deviations  and  coefficients  of  variation  for 
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moans  do  not  appear  to  fit  any  pattern. 

The  mean  squares  of  the  analysis  of  variance  of 
population  means  for  the  original  and  logto  transformed  data 
are  presented  in  Table  53.  The  analysis  of  variance  of  the 
population  means  indicated  the  blocks  were  not  significantly 
different  while  populations  were  highly  significantly 
different.  Additive,  dominance  and  epistatic  effects  were 
also  significant  with  dominance  by  dominance  interactions 
accounting  for  the  major  epistatic  variation.  The  test  for 
non-additivity  revealed  that  the  log10  scale  was  the  most 
appropriate  due  to  the  effects  of  the  bottle  jar  design 
(discussed  earlier  in  the  diallel  section). 

b.  Experiment  V. 

The  means  of  the  six  populations  created  from  a  cross 
between  Ds2  ans  Dsp  (Table  1)  revealed  that  the  Ft  mean  was 
larger  than  the  other  means  indicating  the  presence  of  Fj 
heterosis  in  this  population  (Table  54).  The  backcross 
population  means  tended  to  exceed  the  Fi  mean  which  is 
unexpected  for  an  additive  model.  These  results  suggest 
epistatic  interactions.  The  mean  for  Dsp  was  found  to  be 
much  less  than  expected.  Observing  the  F2  population 
distribution  for  the  original  and  loglo  dry  weight 
accumulation  it  was  noted  that  the  log10  data  was  skewed 
toward  the  upper  end  of  the  distribution  (Figure  11).  These 
distributions  may  be  interpreted  as  a  confounding  of  two 


distributions  (Figure  12).  The  phenomenon  producing  the  low 


.  r  •  ib 


Table  53.  Mean  squares  from  the  genetic  analysis  of  variance  of 
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Figure  11.  Normal  and  log10  distributions  of  the  F2  dry  weight 
values  for  Experiment  V. 
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Figure  12.  Normal  and  log10  distributions  suggesting  the  distributions  of 
segregating  factors  A  and  B  for  Experiment  V. 
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Dsp,  F 2  and  B2  scores  appeared  to  be  a  gene  segregating  for 
pH  effects  since  the  Dsp  line  grew  normally  at  pH  6.2±.2. 
Further  investigation  would  be  necessary  to  determine  the 
exact  nature  of  the  low  Dsp  scores* 

The  six  populations  of  a  cross  between  Ds2  and  Dsp  when 
analysed  for  the  analysis  of  variance  of  population  means 
showed  highly  significant  population  differences  (Table  55)* 
Additive,  dominance  and  epistatic  effects  were  highly 
significant.  Dominance  X  dominance  and  additive  X  dominance 
interactions  contributed  to  epistatic  variation.  These 
results  were  compatible  with  those  of  Experiment  IV. 

c.  Experiment  VI. 

The  combined  analysis  of  the  diallel  suggested  that 
three  crosses  showed  high  mean  yield  and  specific  combining 
ability  for  the  characters  studied.  Cross  Ds2  x  Ds5  was  the 
cross  which  had  the  highest  mean  yield  for  the  acetylene 
reduction  assay  so  in  Experiment  VI  the  six  populations  of 
the  Ds 2  x  Ds5  cross  were  created  for  intensive  study. 

The  correlation  cofficients  for  the  characters  studied 
are  presented  in  Table  56.  The  estimates  of  N2  fixation  were 
found  to  be  highly  correlated  among  themselves  and  partially 
correlated  with  the  leaf  characters.  Applying  the  principal 
components  transformation  to  the  fixation  estimate 
characters  yielded  one  factor  (Table  57).  The  estimate 
variables  were  highly  correlated  with  the  first  factor 


indicating  that  all  estimates  measure  the  principal  axis  of 


' 


Table  55*  Mean  squares  from  a  genetic  analysis  of  variance  of  population 
means  for  Experiment  V. 
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Table  56*  Correlations  for  Experiment  VI 
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Table  57.  Eigenvectors  for  factors  fro.  Experiment  VI. 


Factor 


Charac  ter 

1 

2 

3 

4 

Acetylene  reduced 

0.90 

-.35 

0.  24 

0.00 

Top  fresh  weight 

0.  97 

-.01 

-.22 

0.11 

Root  fresh  weight 

0.90 

0.38 

0.20 

0.00 

Top  dry  weight 

0.97 

-.02 

-.20 

-.11 

Ei genvalue 

3.51 

0.27 

0.19 

0.02 

%  variation  explained 

88 

7 

5 

0 

' 


(\ .  I 


144 


variation;  namely  N2  fixation.  Shoot  fresh  and  dry  weight 
appeared  to  be  the  better  estimates  of  N2  fixation  since 
they  are  most  highly  correlated  with  the  principal  axis. 
Since  scores  for  the  principal  axis  can  be  defined  with  a 
principal  components  transformation,  the  "factor"  scores 
were  calculated  for  the  principal  axis  from  the 
eigenvectors.  The  new  derived  •character*  has  been  termed 
the  "Nitrogen  fixation"  factor  and  differences  among 
population  means  are  of  importance  rather  than  the  amount  of 
N2  fixed.  The  "Nitrogen  fixation"  factor  has  been  given  a 
distribution  of  1/3  the  mean  and  standard  deviation  of  the 
acetylene  reduction  assay  since  nitrogenase  is  assumed  to 
theoretically  reduce  1/3  the  amount  of  N2  than  acetylene 
(Burris  1974).  The  "Nitrogen  fixation"  factor  was  added  to 
the  data  set  and  treated  as  another  character  for  the 
remaining  analyses  in  Experiment  VI. 

The  means  and  standard  deviations  for  the  original  and 
^ogl0  transformed  data  are  presented  in  Tables  58  and  59 T 
respectively.  The  Fi  mean  was  found  to  be  greater  than 
either  parent  for  all  scales  and  characters  while  the  F2  and 
backcross  means  were  not  significantly  different  from  the 
parents.  The  Ds2  mean  was  found  to  be  significantly  smaller 
than  the  Bi  mean  for  the  cotyledon  characters  (Table  58). 
Both  the  Ds2  and  Bt  means  were  significantly  less  than  the 
other  population  means  suggesting  a  genetic  basis  to  the 
difference.  A  similar  pattern  has  also  been  shown  for  the 


juvenile  leaf  character 


..  .1  »onl*  '«**•'  V  . . 

\  »  ,  -  n't  y  J  I J  3ft  ^  o« 


145 


d 

c 

•H 

•H 

0 


c  c  = 
0>  0  (* 
Mi  «H  O 
0  4*  4* 
(<  d  O 
+*  K  4J 
-4  *H  <H 
2  <H  = 


u 

£ 

0 

d 

0 

£ 

0 

£ 

£ 

d 

«■» 

PS 

co 

¥■4 

O 

r* 

in 

£ 

rr 

Tf 

in 

in 

in 

w 

w 

w 

w 

w 

CS 

r* 

o 

c 

l-t 

CS 

P“l 

pH 

¥■4 

¥■4 

44 

44 

44 

+1 

•H 

+1 

03 

f- 

O’ 

o 

03 

CS 

pH 

pH 

CS 

(N 

CS 

9 

£ 

4-» 

U 

0 

<H 

C 

0 


d 

•H 

u 

a) 

> 

4H 

O 

CD 

V 

c 

9 


<w 

o 

o 

o 

*0 

c 

(0  • 

M 

*  > 

w 

c  4- 

0  c 

-H  A) 


—  —  O'  •**  -M  JJ 

»  ^  ««* 


~  I 


•H 

o 

0 

£ 


V  r-» 

£  a 

m  ^ 

K*)  *H  U) 

bos 


W  X)  Z  W 


-H  *■* 

£  £  t-> 

+  w  an  a 

0  &  **4  v 
0  k  M 
»  'H  *  w 


+>  £  £  a 

0  03  tHi  V. 

o  a-  «h  ot 

£  t  9  S 
W  <*  *  w 


w  w  w  w  w  w 


£ 

£ 

d 

£ 

£ 

£ 

<■* 

03 

r* 

in 

03 

in 

CO 

co 

in 

’'t 

w 

w 

w 

w 

w 

w 

00 

r- 

cs 

o 

CO 

r* 

't 

£ 

£ 

r* 

+i 

+i 

+1 

-H 

-H 

+i 

r* 

CO 

00 

00 

03 

co 

cs 

CO 

CO 

trH 

tH 

1-H 

H 

tJ 

co 

£ 

SP 

£ 

00 

o 

-H 

r* 

CO 


"t 

U) 

w 

■cf 

00 

CS 

-H 

IT) 

<N 


t ) 

I 

£ 

oc 

CO 

hp 

CO 

>0 

o 

44 

£ 

£ 


£ 

CO 

w 

o 

44 

l£ 

t" 


£ 

d 

<»» 

co 

TT 

w 

CS 

00 

• 

o 

+1 

CM 

a 


<o 

£ 

CO 

w 

CS 

CO 

CS 

4H 

03 

If) 


0 

£ 

03 

w 

in 

00 

o 

-H 

CO 

O' 


CD 

pH 

CS 

44 

r- 

co 


t > 
o 

PS 

If) 

in 

w 

OJ 

00 

• 

o 

-H 

03 


r» 

TT 

sp 

co 

CS 

44 

CS 

CS 


d 

^ s 

1/3 

in 


+i 

£ 

O 

CS 


£ 

w 

£ 

£ 

CS 

-H 

r* 


E 

in 

Tf 

£ 

in 

in 

in 

d 

•*4 

•H 

u 

> 

9 

9 

a 

V 

x 

£ 

£ 

w 

d 

£ 

d 

£ 

£ 

d 

t3 

9 

C 

P'S 

5- 

C 

£ 

CO 

in 

CS 

CS 

£ 

d 

0 

9 

*0 

X. 

in 

in 

in 

in 

tJ 

pJ 

9 

w 

w 

w 

w 

w 

w 

C 

>> 

0 

a 

r* 

cs 

o 

o 

o 

in 

d 

CD 

V 

3 

\ 

co 

cs 

CO 

co 

CO 

co 

V 

9 

tJ 

ss 

-H 

■H 

-H 

•n 

•H 

-H 

(0 

• 

0 

9 

c 

in 

CS 

03 

00 

CO 

+- 

<5 

£ 

in 

in 

in 

£ 

*¥ 

0 

CD 

d 

c 

h 

d 

d 

c 

00 

in 

0) 

0) 

00 

PH 

9 

£ 

in 

in 

in 

cs 

o 

pH 

0 

cs 

pH 

pH 

c 

• 

0 

00 

•H 

in 

+* 

d 

cs 

in 

9 

r-» 

CD 

CD 

p-* 

3 

O 

O 

£ 

a 

w 

w 

d 

0 

*4 

IV 

•4 

(V 

*H 

(V 

H 

cu 

Oh 

(X 

b 

lx 

GQ 

GQ 

1 

c 

«j 

s 

* 

0) 

2 
I 

V 
I  C 
I  9 

I  -0 
I  3 
I  * 

I  C/3 

I 

I  CD 

4- 

U 

0 

<H 


C 

a/ 

ti 

& 

<H 

tH 

*C 

>> 

pH 

4-* 

c 
<0 
0 
I  *H 
!  'H 


•H 

c 

QL 

•H 

to 


4-> 

0 

c 


9 
$- 
I  d 
I 

tl 

0/ 

V 

4- 

9 


9 

s 

<0 

CD 

9 


£ 

>> 

£ 

• 

*0 

4- 

9 

(0 

» 

9 

0 

4H 

r-> 

pH 

9 

0 

dO 

94 

c 

d 

CD 

X 

C 

d 

9 

0 

E 

9 

9 

3d 

Table  58  (cont)*  Means*  standard  deviations*  and  coefficients  of 
variation  for  the  original  characters  for  Experiment  VI 
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The  analyses  of  variance  of  the  population  means  for 
the  original  and  the  logl0  transformed  characters  are  given 
m  Tables  60  and  61,  respectively.  Populations  and  dominance 
effects  were  highly  significant  for  all  characters  in  both 
scales;  however,  additive  effects  were  significant  for  root 
fresh  weight.  Epistasls  appeared  to  be  absent  for  most 
characters.  These  data  suggest  a  completely  non-heri table 

model  for  the  estimates  of  total  symbiotic  nitrogen 
fixation. 

The  analyses  of  variance  of  the  population  means  for 
the  cotyledon  and  juvenile  leaf  characters  indicated  that  an 
additive-dominance  model  adequately  explains  the 
distribution  of  the  means.  These  data  suggest  that  large 
cotyledons  and  juvenile  leaves  appear  to  be  conditioned  by  a 
dominant  allele  at  a  single  locus  or  a  block  of  loci 
segregating  as  a  single  unit.  To  determine  whether  the 
cotyledon  and  juvenile  leaf  size  expressed  pleiotropism  the 
cotyledon  area  and  juvenile  leaf  width  characters  were 
subdivided  into  groups.  The  crosstabulation  of  the  groups 
indicated  that  plants  with  small  cotyledons  and  large 
juvenile  leaves  and  Y.X.Q.S.  versa  were  found  implying  that  the 
cotyledon  and  juvenile  leaf  sizes  are  controlled  by 
different  genes  (Table  62). 

Having  defined  the  cotyledon  area  and  juvenile  leaf 
width  groups,  it  would  be  of  interest  to  determine  if  plants 
with  larger  cotyledons  and  juvenile  leaves  would  fix  more 


nitrogen.  The  means,  standard  deviations, 


and  coefficients 


' 


*  : 


Table  60.  Mean  squares  for  the  original  characters  for  Experiment  VI 
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Table  60  (cont).  Mean  squares  for  the  original  characters  for 
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Table  61*  Mean  squares  for  the  logl0  transformed  characters  for  Experiment  VI 
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Table  62* 


Cotyledon 

area 


<  32  mm 
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39—42  mm 
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total 


Numbers  of  plants  for  cotyledon  area  by  juvenil 
leaf  width  cr oss tabula ti on  for  Experiment  VI* 
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of  variation  for  the  nitrogen  fixation  estimates  are 
presented  in  Table  63.  Plants  with  larger  cotyledons  and 
juvenile  leaves  fixed  significantly  more  nitrogen.  It  was 
also  noted  that  plants  growing  in  fine  sand  also  fixed  more 
nitrogen  than  did  plants  growing  in  very  fine  sand  (Table 
63).  Combining  the  above  groupings  with  populations,  a 
factorial  analysis  of  variance  was  performed  on  the 

estimates  of  nitrogen  fixation  to  determine  whether 
population  interactions  existed.  These  analyses  of  the 
original  and  logio  transformed  data  are  presented  in  Tables 
64  and  65,  respectively.  All  main  effects  were  highly 
significant  for  all  estimates  except  root  fresh  weight  for 
populations  when  adjusted  for  the  other  main  effects.  The 
only  interaction  of  significance  was  population  by  sand  size 
indicating  that  changing  sand  size  in  an  experiment  may 
affect  the  estimates  of  nitrogen  fixation.  Since  cotyledon 
area  and  juvenile  leaf  width  groups  show  significant 
differences  for  the  total  fixation  estimates,  it  would  be  of 
interest  to  treat  cotyledon  area  and  juvenile  leaf  width  as 
covariates  to  adjust  the  estimates  for  effects  due  to 
cotyledon  or  juvenile  leaf  size  (juvenile  plant  vigor) 
because  these  leaves  are  functional  during  the  onset  of 
nodulation.  Because  plants  with  larger  cotyledons  and 
juvenile  leaves  may  fix  more  nitrogen  (Mytton,  1973),  this 
advantage  should  be  removed  from  the  analysis  to  understand 
the  mechanisms  of  N2  fixation.  The  analysis  of  covariance 


for  the  estimates  for  the  original  and  logio  transformed 


,  -  si  I 


t  ! 


♦ 


Table  63*  Means*  standard  deviations  end  coefficients  of  variati 
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data  are  presented  in  Tables  66  and  67,  respectively.  As 
expected  the  covarlates  were  highly  significant  as  were 
differences  among  the  populations  (see  Table  66).  However, 
dominance  deviations  were  no  longer  significant  whereas 
additive  effects  have  become  highly  significant.  These  data 
suggest  that  the  "efficiency"  (amount  of  fixation  following 
nodule  formation)  of  nitrogen  fixation  appears  to  be 
heritable.  Examining  the  population  means  following 
adjustment  for  the  covariates  revealed  that  the  Ds2  mean  was 

adjusted  upward  (Table  68)  as  expected  from  earlier  data 
(Table  60). 

Since  Ds2  was  shown  to  be  similar  to  Ds5  for  the 
unadjusted  estimate  means  and  significantly  different  for 
cotyledon  and  juvenile  leaf  size  it  may  be  assumed  that  Ds2 
is  more  "efficient"  for  total  fixation  per  unit  area  of 
cotyledon  or  juvenile  leaf  than  is  Ds5. 

A  method  to  investigate  this  "efficiency"  of  fixation 
would  be  to  divide  the  nitrogen  fixation  estimates  by  either 
the  cotyledon  area  or  juvenile  leaf  width  to  create 
characters  of  ratios  of  fixation  per  unit  of  cotyledon  area 
or  juvenile  leaf  width.  The  population  means  for  the 
estimate  ratios  have  been  given  in  Tables  69  and  70  for 
cotyledon  area  and  juvenile  leaf  width  adjustment, 
respectively.  For  the  estimate  ratios,  Ds2  was  again  found 
to  be  significantly  more  "efficient"  than  Ds5  for  ell 
estimates  except  root  fresh  weight.  The  analysis  of  variance 


for  the  estimate  ratios  revealed  that  total  genetic  effects 
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Table  67*  Mean  squares  from  an  analysis  of  covariance  of  the  log»o 
transformed  characters  for  Experiment  VI. 
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Table  68.  Adjusted  means  from  an  analysis  of  covariance  for 
Experiment  VI. 


Popula  tl on 

n 

Acetylene 
r educed 
( nm/ pi/ hr ) 

Shoot 

fresh 

weight 

(  mg/ pi  ) 

Foot 

fresh 
weight 
( g/pl) 

Shoot 

dry 

weight 

( mg/ pi  ) 

Nitrogen 

fixation 

"factor" 

origina 1 

Pl(  Ds2  ) 

58 

78.0 

657 

i 

* 

00 

00 

1 

1 

173 

24.9 

P2(  Ds5  ) 

55 

47.2 

432 

1  .47 

106 

15.  4 

Fi 

59 

68.  7 

617 

1.78 

164 

23.0 

f2 

229 

57.  1 

519 

1  .67 

134 

19.0 

Bi 

108 

65.0 

596 

1 .77 

158 

22.  1 

b2 

111 

58.  1 

521 

1 .85 

135 

20.  3 

lof!l  O 

Pl( Ds2  ) 

58 

1.85 

2.79 

0.21 

2.20 

1.35 

P2(  Ds5  ) 

55 

1.64 

2.  61 

0.  14 

2.00 

1.15 

F  i 

59 

1.78 

2.75 

0.20 

2.17 

1.30 

f2 

229 

1.70 

2.  68 

0.  17 

2.08 

1.22 

Bl 

108 

1.75 

2.  75 

0.  19 

2.15 

1.29 

b2 

111 

1.68 

2.  66 

0.  19 

2.06 

1.20 

. 


00  | 

Table  69.  Means*  standard  deviations*  and  coefficients  of  variation  for  N  fixation 
estimates/cotyledon  area  for  Experiment  VI. 
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were  additive  for  cotyledon  adjusted  date  (Table  71)  whereas 
dominance  effects  were  shown  to  be  present  for  shoot  weights 
and  the  "Nitrogen  fixation”  factor  for  the  juvenile  leaf 
adjusted  data  (Table  72).  Grouping  the  shoot  dry  weight 
estimate  ratios  and  crosstabulatlng  with  either  cotyledon 
area  or  juvenile  leaf  groups,  revealed  the  presence  of 
plants  with  large  cotyledons  or  juvenile  leaves  with 
"efficiency"  of  fixation  and  yicg,  y£££a  (Tables  73  and  74). 
It  would  appear  to  be  possible  to  select  plants  with 
vigorous  seedlings  which  were  "efficient"  for  fixation. 

D.  Host  genes  for  symbiotic  nitrogen  fixation  in  Desmodrlpm. 

The  genes  controlling  the  formation  of  nitrogenase  are 
found  on  the  Elilzo^ium  genome  (Bergerson,  1974).  Since  each 
experiment  was  inoculated  with  what  was  assumed  to  be  a 
homogeneous  culture  then  the  nitrogenase  enzyme 
should  be  uniform  within  each  experiment.  Significant 
genotype  effects  may  indicate  that  the  host  plants  differed 
in  their  ability  to  support  nitrogen  fixation 
physiologically.  Nitrogenase  would  not  be  expected  to  be  the 
limiting  factor  since  it  was  expected  to  be  the  same  for  the 
"efficient"  as  for  the  "inefficient"  genotypes.  It  must  be 
concluded  that  there  are  genes  within  Desmodium  which  may 
indirectly  regulate  nitrogen  fixation.  Increasing  the 
fixation  rate  of  the  Desmodium  -  Eki^okiiim  interaction  may 
not  be  completely  effective  by  solely  manipulating  the 


Rh izob jum  strains  which  has  been  the  traditional  practice. 


. 


Table  71.  Mean  squares  for  the  N  fixation  estimetes/cotyledon  area  for 
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Table  73. 


Numbers  of  plants  for  shoot  dry  weight/cotyledon 
area  ratio  group  by  cotyledon  area  or  Juvenile 
leaf  arid  th  c  r  oss  tabu  lat  i  on  for  Experiment  VI. 


Shoot  dry 

weight /cotyled 

on 

area 

ratio 

< 

3.5 

3. 5-7.0 

> 

7.0 

Cotyledon  area 

Row 

total 

I 

I 

I 

I 

<  32  mm2 

I 

66 

I 

67 

I 

24 

I 

157 

I 

I 

I 

I 

I 

I 

I 

I 

32—39  mm2 

I 

74 

I 

52 

I 

4 

I 

130 

I 

I 

I 

I 

I 

I 

I 

I 

39—42  mm2 

I 

95 

I 

73 

I 

7 

I 

175 

I 

I 

I 

I 

I 

I 

I 

I 

>  42  mm2 

I 

82 

I 

76 

I 

4 

I 

162 

I 

I 

I 

I 

Juvenile  leaf 

widt  h 

Row 

tota  1 

I 

I 

I 

I 

<  6  mm 

I 

103 

I 

73 

I 

17 

I 

193 

I 

I 

I 

I 

I 

I 

I 

I 

6—7  mm 

I 

133 

I 

103 

I 

16 

I 

252 

I 

I 

I 

I 

I 

I 

I 

I 

>  7  mm 

I 

80 

I 

92 

I 

6 

I 

178 

I 

I 

I 

I 

Column 

total 


317 


26  8 


39 


624 
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Table  74« 


Numbers  of  plants  for  shoot  dry  weight/ juvenile 
leaf  width  ratio  group  by  cotyledon  area  or 

juvenile  leaf  width  c ross tabulati on  for  Experiment 
VI  • 


Shoot  dry  we ight/ juvenile  leaf  width  ratio 


< 

20 

20-40 

>  40 

Cotyledon  area 

Row 

to  ta  l 

I 

I 

I 

I 

<  32  mm2 

I 

88 

I 

57 

I 

10 

I 

155 

I 

I 

I 

I 

I 

I 

I 

I 

32—39  mm2 

I 

88 

I 

37 

I 

4 

I 

129 

I 

I 

I 

I 

I 

I 

I 

I 

39—42  mm2 

I 

91 

I 

77 

I 

7 

I 

175 

I 

I 

I 

I 

I 

I 

I 

I 

>  42  mm2 

I 

74 

I 

83 

I 

5 

I 

162 

I 

I 

I 

I 

Juvenile  leaf 

width 

Row 

tota  l 

I 

T 

I 

I 

<  6  mm 

I 

101 

I 

80 

I 

11 

I 

192 

I 

I 

I 

I 

I 

I 

I 

I 

6—7  mm 

I 

143 

I 

98 

I 

11 

I 

252 

I 

I 

I 

I 

I 

I 

I 

I 

>  7  mm 

I 

98 

I 

76 

I 

4 

I 

178 

I 

I 

I 

I 

Column 

total 


342 


254 


26 


622 
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All  analyses  of  the  unadjusted  estimates  of  nitrogen 
fixation  revealed  that  all  Ft  crosses  exhibited  heterosis 
when  the  plants  were  grown  in  a  nitrogen  limiting 
environment.  However,  when  nitrogen  is  supplied  in  adequate 
amounts,  there  appears  to  be  no  difference  between  the  Fi  or 
parental  plants  (McWhirter,  1969;  Rotar  and  Chow,  1971).  The 
heterosis  in  this  study  may  be  simply  the  expression  of 
parental  inbreeding  depression  when  the  plants  are  grown  in 
a  nitrogen  limiting  environment.  These  conclusions  are 
supported  by  the  negative  estimates  of  the  parental  S.C.A. 
effects  from  the  diallel  studies  (Table  50). 

E.  Breeding  procedures  for  nitrogen  fixation  in  Desmodiym. 

The  analyses  of  this  study  revealed  two  major  points 
concerning  the  inheritance  of  symbiotic  nitrogen  fixation  in 
Desmodjurn.  These  hypotheses  are:(a)  that  there  was  a  major 
dominant  gene  conditioning  strain-specific  Ineffective 
nodulat i on ,  and  (b)  that  there  was  significant  genetic 
variation  for  total  symbiotic  nitrogen  fixation  among 
effectively  nodulated  Be s mod?  urn  genotypes  and  that  this 
variation  appears  to  have  a  heritable  basis.  Partitioning 
this  genetic  variance  suggested  that  additive  genetic 
variation  and  dominance  variance  were  present  with  dominance 
variance  being  much  larger  when  the  estimates  of  nitrogen 
are  not  adjusted  for  juvenile  plant  vigor.  Adjusting  the 
estimates  of  nitrogen  fixation  for  juvenile  plant  vigor 
revealed  an  almost  totally  additive  (heritable)  basis  to  the 
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genetic  variation. 


1.  Strain-specific  ineffective  nodulation  genes. 

The  presence  of  strain-specific  ineffective  nodulation 
genes  may  be  construed  as  a  disadvantage  when  breeding  for 
increased  nitrogen  fixation  in  Dg_s.jn2dium.  However,  these 
genes  may  also  be  used  to  advantage.  If  a  Degmod^im  line  or 
hybrid  were  to  be  released  as  a  commercial  cultivar  it  would 
be  necessary  to  eliminate  major  genes  causing  ineffective 


nodulation. 


rapid  screening  program  with  recurrent 


selection  for  the  recessive  genotype  (effective  nodulation) 
during  early  generation  breeding  trials  would  quickly  remove 
ineffective  nodulation  genes  from  the  potential  cultivars. 

The  ineffective  nodulation  genes  may  also  be  used  by  a 
microbiologist  as  a  means  of  selecting  for  nodulation 
competitiveness  among  strains  of  Rhizoblum  (Vincent,  1974). 
Having  selected  a  series  of  Desniogti um  tester  lines  which 
carry  at  least  one  strain-specific  ineffective  nodulation 
gene,  a  microbiologist  may  initiate  competition  studies  to 
rank  the  strains  from  most  to  least  competitive  for  the 
nodulation  sites.  The  most  competitive  strain  if  desirable 
in  all  other  characteristics  would  be  released  as  the 
commercial  inoculum.  The  plant  breeder  having  available  a 
desirable  strain  of  ghizoblum  would  begin  to  select  or  breed 
for  a  Desmodium  hybrid  or  inbred  line  which  would  maximize 


the  total  fixation  by  the 


line  -  commercial  strain 


interact  ion. 
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2*  Increasing  -total  symbiotic  nitrogen  fixation* 

The  outstanding  feature  of  the  data  from  the 
quantitative  experiments  was  the  expression  of  heterosis* 
Those  data  clearly  indicated  that  heterozygous  genotypes  of 
HSS.JSSldi.ujD  were  significantly  superior  to  inbred  lines  for 
total  symbiotic  nitrogen  fixation.  A  hybrid  breeding  program 
exploiting  this  heterosis  would  be  desirable* 

The  reporting  of  cytoplasmic-genetic  male-sterility  in 
MsmUdi  (McWhlrter,  1969)  greatly  assists  the  initiation 
of  a  hybrid  breeding  program.  The  male-sterile  factor  was 
found  in  line  Dsl  and  fertility  was  partially  restored  in 
crosses  with  line  Dal*  Male— sterile  factors  were  also  found 
in  lines  Ds2  and  Da2.  However  fertility  restoration  was  poor 
(McWhirter,  pers.  com.)*  Before  Initiating  a  hybrid  breeding 
program  for  total  fixation*  complete  fertility  restoration 
would  have  to  be  established* 

However,  with  the  presence  of  significant  average 
heterosis  the  accumulation  of  efficient  inbred  lines  should 
be  accompanied  by  a  hybrid  testing  program*  Such  a  breeding 
procedure  has  been  referred  to  as  reciprocal  recurrent 
selection  (Comstock  git  al.  ,  1949).  Base  populations  for  use 
in  a  reciprocal  recurrent  selection  program  could  involve 
the  phe no typically  distinct  yet  fully  fertile  D*  sandwicep^e 
by  H*  or  D.  ag.ar.lna.g  by  £.  uncjnatum  crosses. 
Inbred  lines  necessary  for  this  program  could  be  developed 
either  by  mass  selection,  l.e*  preservation  of  additive 


effects  or  by  single  seed  descent  (Goulden,  1939)* 


. 


"  '  ' 
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3*  Nil: rogen  fixation  components* 

Heterosis  or  dominance  genetic  effects  appear  to  be  the  rule 
when  the  nitrogen  fixation  estimates  are  analysed  in  their 
totality*  However,  when  the  total  varietion  for  each 
estimate  is  partitioned  into  components  the  dominance 
effects  disappear*  The  genetic  variation  of  each  part  is 
almost  exclusively  additive  (heritable)*  In  Experiment  VI, 
two  components  of  total  nitrogen  fixed  were  defined  as;  (a) 
Juvenile  plant  vigor;  estimated  by  cotyledon  area  and 
Juvenile  leaf  width,  and  (b)  "efficiency”  of  nitrogen 
fixation  as  estimated  by  the  ratio  of  total  nitrogen  fixed 
to  Juvenile  plant  vigor*  Large  cotyledon  and  Juvenile  leaf 
size  appeared  to  be  conditioned  by  a  dominant  allele  at  a 
single  locus*  Plants  with  large  cotyledons  were  shown  to  fix 
significantly  more  nitrogen  so  it  would  be  advantageous  at 
the  early  stages  of  a  breeding  program  to  select  for  large 
cotyledons  or  Juvenile  leaves*  These  plants  would  probably 
be  more  vigorous  and  therefore  nodulate  more  quickly  or  form 
nodules  at  more  sites*  These  results  ere  similar  to  the 
results  reported  by  Mytton  (  1973)  where  large  seed  weight 
assisted  nodulation*  It  would  be  expected  that  selection  for 
large  cotyledons  and  Juvenile  leaves  in  Desmod j  igp  would 
rapidly  respond  to  selection  since  a  single  gene  appears  to 


control  these  characters*  Also  achieved  with  this  selection 
would  be  plants  with  increased  seedling  vigor  which  would 
establish  a  sward  and  begin  fixing  nitrogen  and  dry  matter 
more  quickly*  Having  selected  a  population  of  seedlings  with 


. 
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large  cotyledons  or  juvenile  leaves  the  plants  would  enter 
the  second  phase  of  the  breeding  program;  namely  selection 
for  "efficiency"  of  nitrogen  fixation.  Highly  efficient 
plants  would  by  tested  against  the  selected  Rhizchf um 
strains  (discussed  in  part  1  of  this  section)  to  determine 
the  most  favorable  interaction  which  would  be  competitive 
for  noduiation  sites  under  field  conditions.  Such  superior 
interactions  could  be  selected  for  superior  general 
adaptability  or  superior  specific  adaptability  to  a 
particular  location  which  exhibits  unique  environmental 
conditions.  If  an  interaction  was  selected  for  general 
adaptability  this  interaction  would  have  to  be  tested  under 
many  locations*  sites*  and/or  years  since  environmental 
conditions  greatly  effect  total  nitrogen  fixation. 

At  this  early  stage  in  the  development  of  the  study  of 
host  variation  for  nitrogen  fixation,  the  best  avenue  to 
pursue  to  achieve  maximum  return  for  research  input  is  not 
clear.  A  hybrid  breeding  program  working  directly  with 
heterosis  may  be  more  effective  where  capital  and  facilities 
are  not  limiting  since  the  heritable  variation  is  maintained 
under  reciprocal  recurrent  selection.  The  limitation  of  the 
breeding  program  would  be  the  numbers  of  individual  plants 


which  would  have  to  be  maintained  to  insure  that  to  progra 


m 


was  effective. 


Initiati ng 


breeding  program  emphasising  the 


components  of  total  nitrogen  fixation  would  be  more 
advantageous  for  laboratories  with  limiting  resources.  After 


: 


, 
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selection  for  vigorous  juvenile  plants  the  populations  could 
be  effectively  reduced  by  the  elimination  of  the  non- 
vigorous  plants*  The  researcher  could  concentrate  his 
resources  to  a  smaller  number  of  plants  because  the  genetic 
basis  of  the  components  appears  to  be  heritable  and  a 
selection  response  for  greater  "efficiency"  would  be 
expected* 

It  may  be  easy  to  speculate  about  the  expected  results 
of  a  breeding  program  employing  either  approach  but  the 
ultimate  criterion  judging  the  success  of  any  research  for 
increased  nitrogen  fixation  would  be  a  selection  response 
culminating  in  a  inbred  variety  or  the  creation  of  a 
superior  hybrid  which  fixes  a  maximum  amount  of  nitrogen 
under  field  conditions*  Only  when  the  products  of  research 
achieve  practical  results  under  agricultural  conditions  is 
the  research  program  of  merit  or  achieving  its  goals* 

F*  Relevance  to  tropical  agriculture* 

The  importance  of  Dggmpdjum  in  tropical  agriculture  is 
two-fold.  Firstly,  Degm odium  being  a  forage  legume  would 
provide  increased  forage  protein  for  animal  consumption 
thereby  effectively  increasing  availability  of  protein  for 
human  consumption*  Secondly,  Pesrood 1 um  as  a  legume,  fixes 
atmospheric  nitrogen  which  assists  n on- legumi nous  plants  to 
produce  increased  yield*  If  the  total  fixation  could  be 
Increased  Pesmodi um  would  be  valuable  in  the  tropics  as  a 


green  manure  crop  and  a  soil  conditioner  as  well  as  a  forage 


. 
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legume  • 


The  pending  population  Increases  in  the  tropics  will 
dictate  that  the  food  chains  be  made  es  efficient  as 
possible,  i,e  direct  plant  to  human  pathways  for  protein.  D. 
^ajidwi.cease  may  be  used  in  this  context  as  a  model  organism 
for  legume  studies  because  of  its  many  attributes  for 
glasshouse  growth.  Principles  of  inheritance  found  in  D. 
Saa^l.Sgnae  may  be  directly  applicable  to  other  self-fertile 
grain  legumes  which  are  used  for  human  consumption.  The 
future  of  the  need  for  adequate  amounts  of  protein  for  human 
consumption  cannot  be  underes ti mated  and  it  is  in  this  area 
where  legumes  may  make  a  great  contribution. 

The  secondary  contribution  of  legumes  would  be  as  a 
rotational  crop  with  cereals  to  buffer  against 
catas trophi es .  In  short,  it  would  be  optimal  to  develop 
annual  grain  legumes  which  fix  more  total  nitrogen  and 
increases  protein  production  and  rotate  these  legumes  with 
the  cereal  crops  which  have  been  bred  to  produce  maximum 
energy. 


' 


SUMMARY  AND  CONCLUSIONS 


The  genetic  stucfy  of  symbiotic  nitrogen  fixation  in  the 
tropical  legume  Deg.modi.uin  was  divided  into  three  sections; 
namely  ,  methodology,  qualitative  genetic  variation,  and 
quantitative  genetic  variation.  The  initial  studies  in  the 
methodology  section  suggested  that  the  bottle  jar  was 
effective  for  studying  nitrogen  fixation  in  Desmo^j  VMn. 
Uninoculated  or  ineffectively  nodulated  control  plants  were 
yellow  and  chlorotic  suggesting  the  lack  of  nitrogen.  Green 
plants  were  taken  from  the  bottle  jars  and  analysed  with  the 
acetylene  reduction  assay.  Dry  weight  accumulation  was  also 
recorded.  Reduction  of  acetylene  to  ethylene  preceeded  at  a 
linear  rate  for  approximately  90  minutes.  Linear  reduction 
indicated  that  senescence  of  harvested  Desmod^u^  plants  was 
not  a  limiting  factor  for  the  acetylene  reduction  analysis. 

The  w i t h i n— bot t l e  jar  variation  among  similar  genotypes 
suggested  that  the  design  of  the  bottle  jar  could  be 
modified  to  reduce  wi thin-bot t le  jar  variation. 
Modifications  to  the  design  and  sand  particle  size  greatly 
reduced  this  variation. 

Initial  nodulat ion  studies  indicated  that  some  D. 
aparlnes  lines  were  ineffectively  nodulated  by  strain  B. 
Studies  of  F2  and  backcross  populations  from  £•  sftndwlgeqge 
x  D.  aparlnes  crosses  implied  that  the  strain-specific 
ineffective  nodulatlon  interaction  between  some  J>.  aparj  peg 
lines  and  strain  B  was  conditioned  by  a  dominant  allele  at  a 
single  locus. 
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The  analysis  of  quantitative  genetic  variation 
suggested  that  several  points  were  of  importance.  Firstly, 
the  correlation  between  dry  weight  accumulation  and  the 
acetylene  reduction  assay  was  highly  significant  and  there 
appeared  to  be  a  consistant  pattern  for  both  characters 
within  an  experiment.  The  acetylene  reduction  may  be  used  as 
a  plant  breeding  tool  to  screen  juvenile  Desmodium  plants. 

The  most  significant  aspect  illustrated  by  all 
quantitative  experiments  was  the  expression  of  Ft  heterosis. 
Fl  hybrids  were  significantly  superior  to  inbred  lines  for 
total  nitrogen  fixed.  A  breeding  program  incorporating 
reciprocal  recurrent  selection  would  maximize  the  effect  of 
Fl  heterosis.  Although  Fi  heterosis  was  significant,  a 
portion  of  the  genetic  variance  also  appeared  to  be 
heritable  and  analysis  of  components  of  total  nitrogen 
fixation  predominantly  expressed  additive  (heritable) 
genetic  variation.  Selecting  inbred  lines  for  these 
heritable  effects  would  be  beneficial  since  these  lines 
could  become  the  basis  for  recurrent  selection  programs  or 
released  as  culivars. 

Since  genes  controlling  the  formation  of  nitrogenase 
are  found  on  the  Ehizobium  genome,  quantitative  variation  in 
total  nitrogen  fixed  may  be  due  to  physiological  differences 
among  host  plants.  The  host  plant  may  indirectly  control 
symbiotic  nitrogen  fixation. 

Genotype  by  environment  interactions  were  of 


significant  importance  suggesting  that  studies  concerning 


■ 
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quantitative  variation  for  symbiotic  nitrogen  fixation 

should  be  conducted  under  differing  environments,  locations, 
and/or  years. 

Breeding  programs  utilizing  qualitative  and 
quantitative  genetic  variation  were  outlined.  Such  a  program 
would  be  useful  for  the  production  of  highly  efficient 
varieties.  The  relevance  of  such  varieties  to  tropical 
agriculture  was  also  discussed. 

The  preceding  analyses  suggested  that  the  null 
hypothesis  is  to  be  rejected.  The  null  hypothesis  was  that 
there  was  no  genetic  variation  in  the  processes  of  the 
symbiotic  relationship  between  Destnodlupi  and  Khizoh^um.  Two 
hypotheses  may  be  postulated.  These  are:  (a)  that  strain- 
specific  ineffective  nodulation  interaction  between  strain  B 
and  some  D.  aa&rineg  lines  was  contitloned  by  a  dominant 
allele  at  a  single  locus,  and  (b)  that  there  was  genetic 
variation  in  symbiotic  nitrogen  fixation  in  Desmodlum . 
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Table  1. 


Hybridlzat ion  and  emasculation  observations  using 
various  washing  solutions  and  flower  preparations* 


Number 

of 

Percent 

Av*  no* 

f  l  owe 

rs 

pod 

seeds 

Treatment  manipulated 

set 

per  pod 

Distilled  water  -  partial  flower* 

53 

50 

1.7 

~  entire  flower 

32 

84 

3.9 

-  emasculated 

45 

0 

2*5%  sugar  solution  -  partial  flower 

28 

60 

3.4 

-  entire  flower 

45 

85 

5.8 

—  emasculated 

63 

5 

2.3 

5*0%  sugar  solution  -  partial  flower 

43 

46 

1.8 

—  emasculated 

18 

0 

+the  tip  half  of  the  petals  removed* 

Table  2.  The  percent 
seeds  per  pod 

pod  set 
obtained 

and  the  average  number  of 
during  a  5x5  diallel  cross. 

Number  of 

Percent 

A  v  •  no. 

N urn b e r  of 

flowers 

pod 

seeds 

seeds 

Pe  rcen  t 

Method 

mani pula  ted 

set 

per  pod 

tested 

hybrids 

Pol linated 

387 

8.2 

4.3 

104 

95 

Emasculat  ed 

60 

6.6 

2.5 

Table  3*  The  correlations  between  environmental  factors  and 
hybridizat ion  results* 


Pe  rcent 
pod  set 

Av.  no. 

of  seeds 

Max.  daily  temp. 

0.  15 

-.39 

Min.  daily  temp. 

0.  30 

-.19 

Mean  dally  temp. 

0.  24 

-.41 

Temp,  at  pollination 

0.  07 

-.  17 

SQRT( Dally  precipitation) 

-.24 

0.19 

SQRT 


square  root 


. 

-* 
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Appendix  B*  Diallel  analysis  computer  program  listing* 

1:  DIALL  subroutine  instructions*  187 
2:  DIALL  subroutine  listing*  191 
3:  DIALL  sample  output*  206 


Acknowledgment!  A  model  I,  method  1  diallel  program  was 
obtained  from  Mr*  P*  Hsu,  Dept*  of  Plant  Science,  The  Univ* 
of  Alberta*  This  program  was  extensively  edited,  expanded 
and  corrected  by  the  author  to  produce  an  all  model  and 
method  diallel  program* 
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Subroutine  name: 
Subroutine  title: 
Programmer : 


DIALL 

Diallel  analysis 
B.  Pinchbeck 


Program  Description 

Source  Language:  Fortran  IV 

Purpose:  The  analysis  of  a  diallel  cross  for  models  I,  II, 

and  A;  methods  1  to  4. 

Description:  DIALL  is  a  algorithm  which  computes  the  diallel 

analysis  of  variance,  components  of  genetic 
variance,  Griffing’s  combining  ability  effects, 
or  Gardner-Eberhart ’ s  constants.  The  subroutine 
generally  follows  the  diallel  analysis  as  outlined 
by  Griffing  (1956)  yet  is  modified  to  include  the 
analysis  of  variance  of  Hayman  (1954)  and  Gardner 
and  Eberhart  (1966)  for  methods  1  and  2,  respect¬ 
ively.  The  diallel  table  (p  X  p)  or  within-block 
diallel  tables  (p  X  p  X  b)  are  calculated  from  the 
formulae : 


X.  . 
ij 


SQRT(bs)  Z  Zijkl 


(for  a  two-way  randomized  block  design  with 
interactions) 


or 


X.  . 
ij 


1 

SQRT(b) 


I  Z 


ijk 


(for  a  two-way  randomized  block  design) 


and 


(Vk 


SQRT(s) 


(Vk 


k=l,b 


and 


S 

g 


where 


p  =  number  of  diallel  parents 

b  =  number  of  blocks  (replications,  locations,  years,  etc.) 
s  =  number  of  subsamples  per  cell 
X  =  diallel  table  entry 
Z  =  individual  cell  entries 

X^  =  the  ij th  entry  in  the  overall  diallel  table 
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(Vk 

Sbg 

(Vk 


g 


=  the  ijth  entry  in  the  diallel  table  of  the  kth  block. 

=  the  sum  of  squares  of  the  block  x  genetic  component 
interaction 

=  the  sum  of  squares  for  the  genetic  component  in  the  kth 
block. 

-  the  sum  of  squares  for  the  genetic  component  for  the  overall 
diallel  table. 


The  above  adjustments  were  incorporated  so  that  the  sum  of  the  sum  of 
squares  of  the  main  effect  genetic  components  would  equal  to  the 
treatment  sum  of  squares  for  the  two-way  analysis  of  variance. 

Griff ing  s  combining  ability  effects  and  Gardner-Eberhart ' s  constants 
are  calculated  from  a  diallel  table  of  genotype  means  as  outlined  by 
Griffing  (1956),  and  Gardner  and  Eberhart  (1966),  respectively . 

Components  of  genetic  variance  are  calculated  from  the  diallel  analysis 
of  variance  such  that : 

/N  /\ 

=  2Vg  (inbred  parents) 

A  A 

=  4V^  (noninbred  parents) 

AAA 

(methods  1  and  3;  inbred  parents) 


VD  =  4V^  +  (methods  1  and  3;  non-inbred  parents) 


V  =  V  (methods  2  and  4;  inbred  parents) 
h  h 

/\  /\ 

=  4V^  (methods  2  and  4;  non-inbred  parents) 


VG  '  VA  +  VD 


V_  =  V  +  V_  (model  A) 

P  G  E 

AAA  /N 

vp  =  VG  +  VBG  +  VE  (model  II) 


and  where  the  estimated  variance  of  the  components  are  adjusted  for 
block  and  subsamples.  (for  example) 

~  i 

V  =  - 7 - r-  (M  -  H,  )  (for  model  II,  method  2) 

g  bs (p+2)  g  h 

whereas  Griffing  (1956)  presents  the  formula: 


A 

V  = 
g 


FTT  (Mg-V 


. 

' 


. 
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and  where 


A 

V  = 
g 

the  general  (G.C.A.)  variance  component 

A 

vh  = 

the  heterosis  (S.C.A.)  variance  component 

II 

U 

<  > 

the  reciprocal  variance  component 

<  > 
> 

II 

the  additive  genetic  variance  component 

II 

<  >Q 

the  dominance  genetic  variance  component 

II 

o 

<  > 

the  genetic  variance  component 

<  > 

II 

the  phenotypic  variance  component 

W<> 

O 

II 

the  block  x  genotype  variance  component 

/\ 

VE  = 

the  environmental  variance  component  (error  mean  square) 

Restrictions : 

The  subroutine  is  limited  to  an  input  diallel  table  matrix  of 
genotype  means  of  10  parents  x  10  parents  x  10  blocks. 

Program  Usage 

Calling  Sequence:  CALL  DIALL  (NP ,  NQ,  NR,  NS,  MO,  ME,  NA,  CH,  IE,  EMS , 
A,  IW,  IP) 

Parameters : 

Name  Purpose 


NP  - 

number  of  parents 

NQ  - 

rank  of  the  third  dimension  of  the  diallel  table 

NR  - 

number  of  blocks  in  the  original  design 

NS  - 

the  number  of  subsamples  per  cell  in  the  original  design 

MO  - 

diallel  model  number 

0  -  fixed  model  I  -  Gardner-Eberhart  constants 

1  -  fixed  model  I  -  Griffing  combining  ability  effects 

2  -  random  model  II  -  components  of  variation 

3  -  mixed  model  A  -  components  of  variation 

ME  - 

diallel  method  number 

1  =  parents,  F^Ts  anc*  recipi°cal  ' s 

' 
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2  =  parents  and  Fi’s 

3  =  Fi's  and  recipiocal  Fi's 

4  =  Fi’s 

0  =  parents  and  Fi’s  (Gardner-Eberhart  Analysis  III) 


NA  - 


a  vector  describing  the  experiment  (20A4) 


CH 


a  vector  describing  the  analysed  character  (10A4) 


IE 


error  degrees  of  freedom 


EMS- 


error  mean  square 


A 


matrix  of  block  sums  if  NP  x  NP  diallel  table  or 
subsample  means  if  NP  x  NP  x  NQ  diallel  table 


IW 


logical  I/O  unit  for  output  statements  WRITE  (IW,  FMT) . . . 


IP 


0  -  if  parents  are  inbred 
1  -  if  parents  are  not  inbred 


REFERENCES 

Gardner,  C.O.  1963.  Estimates  of  genetic  parameter  in  cross-fertilizing 
plants  and  their  implications  in  plant  breeding.  In  Statistical  genetics 
and  plant  breeding.  W.D.  Hanson  and  H.F.  Robinson  (eds.)  NAS-NRC 
Pub.  No.  982:225:252. 

Gardner,  C.O.  and  S.A.  Eberhart.  1966.  An  analysis  and  interpretation 
of  the  variety  cross  diallel  and  related  populations.  Biometrics 
22:439-452. 

Griffing,  B.  1956.  Concept  of  general  and  specific  combining  ability 
in  relation  to  diallel  crossing  systems.  Aust.  J.  Biol.  Sci.  9:463-493. 

Hayman ,  B.I.  1954.  The  analysis  of  variance  of  diallel  tables  Biomatrics 
10:235-244. 

Baker,  R.J.  1978.  Issues  in  Diallel  Analysis.  Crop  Science  18:533-536. 


# 

, 


onononnnonnnonnnnnoonononnnnoo^^ 


191 


SUBROUTINE  DIALL(NP ,NQ,NR,NS,MO,ME,NA,CH, IE,EMS,A, IW, IP) 
DIALLEL  ANALYSIS  SUBROUTINE 
VERSION  1.0,  JUNE  1979 


NP  =  NUMBER  OF  PARENTS  (MAX.  10) 

NQ  =  NUMBER  OF  DIMENSIONS  OF  DIALLEL  TABLE  (MAX.  10) 

(I.E.  THE  RANK  OF  THE  THIRD  DIMENSION  OF  ’A’) 

NR  =  NUMBER  OF  REPLICATIONS  (BLOCKS, LOCATIONS, YEARS) 

IN  THE  ORIGINAL  DESIGN 

NS  =  NUMBER  OF  SUBSAMPLES  PER  CELL  IN  THE  ORIGINAL  DESIGN 
MO  =  DIALLEL  MODEL  NUMBER 

0  =  FIXED  MODEL  I  -  GARDNER-EBERHART  CONSTANTS 

1  =  FIXED  MODEL  I  -  GRIFFING’S  COMBINING  ABILITY  EFFECTS 

2  =  RANDOM  MODEL  II  -  COMPONENTS  OF  VARIANCE 

3  =  MIXED  MODEL  A  -  COMPONENTS  OF  VARIANCE 
ME  =  DIALLEL  METHOD  NUMBER 

0  =  PARENTS  AND  FI’S  (GARDNER-EBERHART  ANALYSIS  III) 

1  =  PARENTS, FI’S, AND  RECIPROCAL  FI’S 

2  =  PARENTS  AND  FI’S 

3  =  FI’S  AND  RECIPROCAL  FI’S 

4  =  FI’S 

NA  =  VECTOR  DESCRIBING  THE  EXPERIMENT  (20A4) 

CH  =  VECTOR  DESCRIBING  THE  CHARACTER  BEING  ANALYSED  ( 10A4) 

IE  =  ERROR  DEGREES  OF  FREEDOM 
EMS  =  ERROR  MEAN  SQUARE 

A  =  MATRIX  OF  BLOCK  SUMS  IF  NP  X  NP  DIALLEL  TABLE  OR 
SUBSAMPLE  MEANS  IF  NP  X  NP  X  NQ  DIALLEL  TABLE 
IV  =  LOGICAL  I/O  UNIT  FOR  OUTPUT  STATEMENTS  :  VRITE( IV.FMT) . . . 
IP  =  0  IF  PARENTS  ARE  INBRED 

1  IF  PARENTS  ARE  NOT  INBRED 


PROGRAMMED  BY  BRIAN  R.  PINCHBECK 


DIMENSION  A (10, 10, 10) ,NA(20) ,CH(10) 

DIMENSION  S( 10, 10) ,R(10, 10) ,X(10) ,Y(1G) ,G(10) ,VS(10) ,VG(16) 

DOUBLE  PRECISION  RN,SN,A,S,R,X,Y,V,G,VS,VG,BT,BA,BD,BDD,BPD,BCD,BR 
& , BMR , BOR , SST , SSA , SSD , SSDD , SSPD , SSCD , SSR , SSMR , SSOR , TOT 
COMMON  W( 10, 10) 

IA=NP-1 

ID=NP* IA/2 

ICD=NP*(NP-3)/2 

I0R=(NP-l)*(NP-2)/2 

IF  (NS.LE.0)  NS= 1 

IF  (ME.EQ.0)  ID= ID- 1 

PP=2.0 

PU= 1 .0 

PPP= 1 .0 

IF  (IP.EQ. 1)  PP=4.0 

IF  (IP.EQ. 1)  PPP= 16.0 

IF  (IP.EQ. 1)  PU=4.0 

IB=NR-1 

IBA=IB*IA 

IBD=IB*ID 

IBCD=IB® I CD 

IBOR=IB*IOR 

AI  =  IA 

DI  =  ID 

CDI=ICD 

ORI= IOR 

BI  =  IB 

BAI= IBA 

BDI= IBD 

BCDI= IBCD 

BORI=IBOR 

SN=NS 

RN=NR 

PN=NP 

BT=0.0 

BA=O.0 

BD=0.0 

BDD=0.0 

BPD=0.0 

BCD=0.0 

BR=0.0 

BMR=0.0 

BOR=0.0 

BMT  =0 . 0 

BTI=0.0 

TOT=0.0 

ESS = EMS* IE 


. 


w>\  3,  <i  "  > 


LME=ME+ 1 

GOTO  ( 190, 10, 190,360,550) ,LME 
10  CONTINUE 

C  METHOD  1  ANALYSIS  OF  VARIANCE 
IT=NP**2-1 
TI  =  IT 

IF  (NQ.LE. 1)  GO  TO  40 

IBT= IT* IB 

BTI=IBT 

DO  30  K= 1 ,NQ 

DO  20  1=1, NP 

DO  20  J= 1 ,NP 

W( I , J) =A( I , J ,K)*DSQRT (SN) 

20  CONTINUE 

CALL  DAV 1 (NP , SST , SSA , SSD , SSDD , SSPD , SSCD , SSR , SSMR , SSOR) 
BT=BT+SST 
BA=BA+SSA 
BD=BD+SSD 
BDD=BDD+SSDD 
BPD=BPD+SSPD 
BCD=BCD+SSCD 
BR=BR+SSR 
BMR=BMR+SSMR 
BOR=BOR+SSOR 
30  CONTINUE 
40  CONTINUE 
DO  50  1=1, NP 
DO  50  J= 1 ,NP 
W(I, J)=0.0 
DO  50  K=1 ,NQ 

W ( I , J) =W ( I , J) +A( I , J , K) ®DSQRT (SN)/DSQRT (RN) 

50  CONTINUE 

CALL  DAV l ( NP , SST , SSA , SSD , SSDD , SSPD , SSCD , SSR , SSMR , SSOR) 

BT=BT-SST 

BA=BA-SSA 

BD=BD-SSD 

BDD=BDD-SSDD 

BPD=BPD-SSPD 

BCD=BCD-SSCD 

BR=BR-SSR 

BMR=BMR-SSMR 

BOR=BOR-SSOR 

SMT=SST/TI 

SMA=SSA/AI 

SMD=SSD/DI 

SMPD=SSPD/AI 

SMCD=SSCD/CD I 

SMR=SSR/DI 

SMMR=SSMR/AI 

SMOR=SSOR/OR I 

FAT=SMT/EMS 

FBT=SMT/BMT 

FAA=SMA/EMS 

FAD=SMD/EMS 

FADD=SSDD/EMS 

FAPD=SMPD/EMS 

FACD=SMCD/EMS 

FAR=SMR/EMS 

FAME =SMMR/EMS 

FAOR=SMOR/EMS 

PAT=FPROB( IT , IE, FAT) 

PAA=FPROB( I A, IE,FAA) 

PAD=FPROB( ID, IE, FAD) 

PADD=FPROB( 1 , IE, SSDD) 

PAPD=FPROB( IA, IE.FAPD) 

PACD=FPROB( ICD, IE,FACD) 

PAR=FPROB( ID, IE .FAR) 

PAMR=FPROB( I A, IE.FAMR) 

PAOR=FPROB( IOR, IE.FAOR) 

IF  (NQ.LE. 1)  GO  TO  60 

BMT=BT/BT I 

BMA=BA/BAI 

BMD=BD/BDI 

BMDD=BDD/BI 

BMPD=BPD/BAI 

BMCD= BCD/BCD I 

BRM=BR/BDI 

BMRM=BMR/BAI 

BORM=BOR/BOR I 

FABT=BMT/EMS 
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C 


C 


F  AB A = BM A/EMS 

FABD=BMD/EMS 

FABDD=BMDD/F.MS 

FABPD=BMPD/EMS 

FABCD=BMCD/EMS 

FABR=BRM/EMS 

FABMR=BMRM/EMS 

FABOR=BORM/EMS 

PABT  =FPROB ( I BT , IE,FABT) 

PABA=FPROB( IBA, IE,FABA) 

PABD=FPROB ( I BD , IE.FABD) 

PABDD=FPROB( IB, IE,FABDD) 

PABPD=FPROB( IBA, IE,FABPD) 

PABCD=FPROB ( I BCD , IE.FABCD) 

PABR=FPROB ( IBD , IE,FABR) 

PABMR=FPROB ( I BA , I E , FABMR ) 

PABOR=FPROB( IBCR, IE.FABOR) 

60  CONTINUE 

WRITE  ( IW,730)  ME,NA,CH 
( IW,720) 

( IW,740) 

( IW, 1020)  IT , SST , SMT , FAT , PAT 
( IW,750)  I A , SSA , SMA , FAA , PAA 
I D , SSD , SMD , FAD , PAD 
SSDD , SSDD , FADD , PADD 
I A , SSPD , SMPD , FAPD , PAPD 
I CD , SSCD , SMCD , FACD , PACD 
I D , SSR , SMR , FAR , PAR 

I A ,  SSMR , SMMR , FAMR , P AMR 
I  OR ,  SSOR ,  S'MOR ,  FAOR ,  P  AOR 

IF  (NQ.LE. 1)  GO  TO  70 
WRITE  ( IW, 830)  I BT , BT , BMT , FABT , PABT 

I BA , BA , BMA , FABA , PABA 
I BD , BD , BMD , FABD , PABD 

I B ,  BDD , BMDD , FABDD , PABDD 
IBA , BPD , BMPD , FABPD , PABPD 
I BCD , BCD , BMCD , FA BCD , PABCD 
I BD , BR , BRM , FABR , P ABR 
I BA , BMR , BMRM , FABMR , PABMR 
I BOR , BOR , BORM , FABOR , P ABOR 


WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 


WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 


( IW,760) 
( IW ,770) 
( IW,780) 
( IW,790) 
( IW , 800) 
( IW,810) 
( IW,820) 


<IW,840> 
( IW,850) 
( IW,860) 
( IW ,870) 
( IW,880) 
( IW,890) 
( IW,900) 
( IW,910) 


70  CONTINUE 

WRITE  ( IW,920) 
WRITE  ( IW,730) 
WRITE  ( IW, 1 160) 
DO  90  1=1, NP 
X ( I ) =0.0 
Y( I) =0.0 


I E , ESS , EMS 
ME,NA,CH 
( I , I = 1 , NP ) 


DO  80  J= 1 ,NP 

W ( I , J) =W( I , J) /DSQRT (SN) /DSQRT (RN) 

X(I)=X(I)+W(I,J) 

Y(I)=Y(I)+W(J,I) 

TOT=TOT+W ( I , J) 

80  CONTINUE 

WRITE  ( IW, 1080)  I, (W<I,J) , J=1,NP) 

90  CONTINUE 

XMU=T0T/PN**2 

SE 1 =SQRT ( EMS/PN**2) 

WRITE  ( IW,  1 130)  XMU , SE  1 
METHOD  1  COMPONENTS  OF  VARIANCE 
,  IF  (M0.LT.2)  GO  TO  100 
WRITE  ( I W , 730 )  ME,NA,CH 

CALL  COV(NP, NR, NS, MO, SMA, SMD, SMR, BMT, EMS, A I ,DI ,BTI, IE,PP ,PU,PPP) 
RETURN 
100  CONTINUE 

IF  (MO.NE.l)  GO  TO  170 

METHOD  1  GRIFFING’S  COMBINING  ABILITY  ANALYSIS 
EMS=EMS/RN/SN 
DO  120  1=1, NP 
R( I , I) =0.0 
VS( I) =0. 

G(I)=((X(I)+Y ( I) )/(2.*PN) )-XMU 

S?i!j?=((W(i!j)+W(J,I))/2.)-((X(I)+X(J)+Y(I)+Y(J))/(2.»PN))+XMU 
S(J, l>  =S< I , J) 

IF  (I-J)  110, 120,110 
110  CONTINUE 

VS(I)=VS(I)+S(I, J)**2 
R(I,J)=(W(I,J)-W(J,I))/2. 

R( J, I)=-R( I , J) 

120  CONTINUE 


c 


c 


SE2.SE3 
(1,1=1 ,NP) 

I,(S(I,J),J=1,I) 

SE4 , SE5 , SE6 , SE7 , SE8 

SE9,SE10 

(1,1=1  ,NP) 


DO  130  1=1, NP 

VG(I)=G(I)*G(I)- (PN- 1 . ) /PN/PN/2 . #EMS 
VS ( I ) = VS ( I ) /PN/2 . - (PN- 1 . ) /PN/2 . *EMS 
130  CONTINUE 

SE2=SQRT ( (EMS* (PN-1 . ) )/2./PN/PN) 

SE3=SQRT (EMS/PN) 

SE4=SQRT ( ( EMS* ( PN- 1 . ) • (PN- 1 . ) ) /PN/PN) 
SE5=SQRT( (EMS* (PN*PN-2.*PN+2. ) )/2. /PN/PN) 
SE6=SQRT ( EMS* 2 . * (PN-2 . ) /PN) 

SE7=SQRT (EMS* (PN- 1 . )/PN) 

SE8=SQRT (EMS* (PN-2. )/PN) 

SE9=SQRT (EMS* (3. *PN-2. )/2./PN) 

SE10=SQRT (EMS* 1.5* (PN-2. )/PN) 

SE 1 1 =SQRT (EMS/2 . ) 

SE12=SQRT (EMS) 

WRITE  (IW.730)  ME.NA.CH 
WRITE  ( IW, 1030) 

WRITE  ( IW, 1040) 

DO  140  1=1, NP 

WRITE  ( IW, 1050)  I ,G( I) , VG( I) , VS( I) 

140  CONTINUE 

WRITE  ( IW, 1060) 

WRITE  ( IW, 1070) 

DO  150  1=1, NP 
WRITE  ( IW, 1080) 

150  CONTINUE 

WRITE  ( I W , 1 090 ) 

WRITE  (IW, 1100) 

WRITE  (IW, 1110) 

DO  160  1=1, NP 
WRITE  ( IW, 1080)  I,(R(I,J),J=1,NP) 

160  CONTINUE 

WRITE  ( IW, i 120)  SE 1 1 , SE 1 2 
RETURN 
170  CONTINUE 

DO  180  1=1, NP 

DO  180  J  = I , NP 

W(I , J)=(W(I, J)+W(J, I>)/2.0 

W(J, I) =W( I , J) 

180  CONTINUE 

GARDNER-EBERHART  CONSTANTS 
WRITE  ( I W,730)  ME,NA,CH 
CALL  GARDCO(NP,ME, IW) 

RETURN 
190  CONTINUE 

METHOD  2  ANALYSIS  OF  VARIANCE 
IT=NP* (NP+ 1 )/2-l 
TI  =  IT 

IF  (NQ.LE.l)  GO  TO  220 

IBT= IT* IB 

BTI=IBT 

DO  210  K= 1 , NQ 

DO  200  1=1, NP 

DO  200  J  = I , NP 

W(I, J)=A(I,J,K)*DSQRT(SN) 

W(J, I)=W(I,J) 

200  CONTINUE 

CALL  DAV2 (NP , ME , SST , SSA , SSD , SSDD , SSPD , SSCD) 

BT=BT+SST 

BA=BA+SSA 


BD=BD+SSD 
^  BDD=BDD+SSDD 
BPD=BPD+SSPD 
BCD = BCD* SSCD 
210  CONTINUE 
220  CONTINUE 

DO  230  1=1, NP 
DO  230  J= I ,NP 
W( I , J) =0.0 
DO  230  K= 1 ,NQ 

W( I , J) =W( I , J) +A( I , J,K) *DSQRT (SN)/DSQRT (RN) 
W( J, I) =W( I , J) 

230  CONTINUE 

CALL  DA V2 ( NP , ME , SST , SSA , SSD , SSDD , SSPD , SSCD ) 

BT=BT-SST 

BD=BD-SSD 

BA=BA-SSA 

BDD=BDD-SSDD 

BPD=BPD-SSPD 

BCD=BCD-SSCD 


SMT=SST/TI 
SMA=SSA/AI 
SMD=SSD/DI 
SMPD=SSPD/A I 
SMCD=SSCD/CDI 
FAT=SMT/EMS 
FAA=SMA/EMS 
FAD=SMD/EMS 
FADD = SSDD/EMS 
FAPD=SMPD/EMS 
FACD=SMCD/EMS 
PAT=FPROB( IT, IE, FAT) 
PAA=FPROB ( I A , I E , FAA) 
PAD=FPROB( ID, IE, FAD) 
PADD=FFROB( 1 , IE, FADD) 
PAPD=FPROB( IA, IE,FAPD) 
PACD=FPROB( I CD, IE, FACD) 

IF  (NQ.LE. 1)  GO  TO  240 
BMT=BT/BTI 
BMA=BA/BAI 
BMD=BD/BDI 
BMDD=BDD/BI 
BMPD=BPD/BAI 
BMCD=BCD/BCD I 
FA3T=BMT/EMS 
FABA=BMA/EMS 
FABD=BMD/EHS 
FABDD=BMDD/EMS 
FABPD=BMPD/EMS 
FABCD=BMCD/EMS 
PABT=FPROB( IBT, IE,FABT) 
PABA=FPROB( I BA, IE.FABA) 
PABD=FPROB ( I BD , IE,FABD) 
PABDD=FPROB (IB, IE,FABDD) 
PABPD=FPROB( I BA, IE, FABPD) 
PABCD=FPROB( IBCD, IE,FABCD) 
240  CONT INUE 


WRITE  ( IW,730)  ME , NA , CH 
WRITE  ( IW,720) 

WRITE  ( IW,740) 

WRITE  ( IW, 1020)  IT , SST , SMT , FAT , PAT 
IF  (ME.EQ.2)  GO  TO  250 
WRITE  (IW,9S0)  I A , SSPD , SMPD , FAPD , P APD 
WRITE  <IW,940)  SSDD.SSDD, FADD, P ADD 
WRITE  ( IW,930)  I D , SSD , SMD , FAD , PAD 
WRITE  ( IW,960)  I A , SSA , SMA , FAA , PAA 
WRITE  ( IW,9S0)  I CD , SSCD , SMCD , FACD , P ACD 
IF  (NQ.LE. I)  GO  TO  260 
WRITE  ( IW,830)  I BT , B I , BMT , FABT , PAST 

(IW, 1010)  I BA , BPD , BMPD , FABPD , PABPD 
(IW, 1000)  I B , BDD , BMDD , FABDD , PABDD 
( IW,990)  I BD , BD , BMD . FABD , P ABD 
IBA,BA,BMA ,FABA,PABA 
I BCD , BCD , EMCD , FABCD , PABCD 


( IW,970) 
(IW,880) 


WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
GO  TO  260 
250  WRITE  ( IW,750) 
WRITE  ( IW ,760) 
WRITE 
WRITE 

WRITE  ( IW, 790) 


I A,  SSA, SKA, FAA, PAA 
ID, SSD, SMD, FAD, PAD 

( IW,770)  SSDD,SSDD, FADD, P ADD 
( IW,780)  I A, SSPD, SMPD, FAPD, P APD 
I CD , SSCD , SMCD , FACD , P ACD 
IF  (NQ.LE. 1)  GO  TO  260 
WRITE  ( IW,830)  I BT , BT , BMT , FABT , PABT 

I BA , BA , BMA , FABA , PABA 
I BD , BD , BMD , FABD , PABD 

IB ,  BDD , BMDD , FABDD , PABDD 
I BA , BPD , BMPD , FABPD , PABPD 
I BCD , BCD , BMCD , FABCD , PABCD 


( IW,840) 
( IW,850) 
( IW,860) 
( IW,870) 
( IW,S80) 


WRITE 
WRITE 
WRITE 
WRITE 
WRITE 

260  CONTINUE 

WRITE  ( IW,920) 
WRITE  ( IW,730) 
WRITE  (IW, 1 170) 
DO  290  1=1, NP 
X(I)=0.0 


I E , ESS , EMS 
ME,NA,CH 
(1,1=1 ,NP) 


DO  270  J  = 1 , NP 

W ( I , J ) = W ( I , J ) /DSQRT ( SN ) /DSQRT ( RN ) 
X(I)=X(I)+W(I, J) 

270  CONTINUE 

DO  280  J  = 1 , I 
TOT=TOT+W( I , J) 

280  CONTINUE 


196 


290 


I.<V(I,J),Jsl,I) 


300 


310 


320 


330 


WRITE  ( IW, 1080) 

CONTINUE 

XMU=T0T*2 . /PN/ (PN+ 1 . ) 

SE1«SQRT(EMS*2./PN/<PN+ 1 . ) ) 

WRITE  ( IW, 1 130)  XMU.SEl 
MET?gD  2  COMPOHENTS  OF  VARIANCE 
IF  (M0.LT.2)  GO  TO  300 
WRITE  <IV,730>  ME , NA , CH 

RETijRNOV<NP’NR’NS,MO'SMA’S>,D’0-0lBMT'EMS-A1-DI-BTI.IE,PP,PU.PPP) 
CONTINUE 
GARDNER-EBERHART  CONSTANTS 
IF  (MO.NE.0)  GO  TO  310 
WRITE  (IW, 730)  ME , NA , CH 
CALL  GARDCO(NP,ME, IW) 

RETURN 
CONTINUE 

METems=1ms/rn/snG  S  combining  ability  ANALYSIS 

DO  320  1=1  NP 
VS(I)=0.0 

nnI^<)5(J)xm(I,I))"(2-*TOT)/PN)/(I>N+2-> 
uu  SZy)  J  =  I ,  NP 

*?!)*(PN+2!)))~<<X<I)+X(J)+W<I’I)+W<J*  J))/(PN+2-))  +  ((2-*T0T)/((PN+ 
S(J,I)=S(I, J) 

IF  (I.EQ.J)  GO  TO  320 
VS(I)=VS(I)+S(I, J)**2 
CONTINUE 
DO  330  1=1, NP 

VG(I)=G(I) **2-( ( (PN- 1 . )/(PN* (PN-2. ) ) )*EMS) 

CONTINUE^' I)/<PN-2* } } ^ <PN-3. )/(PN-2. ) ) *EMS) 

SE2=SQRT ( ( (PN- 1 . )/(PN* (PN+2. ) ) )*EMS) 

SE3=SQRT ( (2 . *EMS) / (PN+2 . ) ) 

SE4=SQRT ( ( (PN* (PN- 1 . ) )/( (PN+ 1 . ) * (PN+2. ) ) ) *EMS) 

SE5=SQRT ( ( (PN**2+PN+2. )/( (PN+ 1 . )* (PN+2. ) ) ) *EMS) 

SE6=SQRT( ( (2. * (PN-2. ) )/(PN+2. ) )*EMS) 

SE7=SQRT ( ( (2. * (PN+ 1 . ) )/(PN+2. ) )*EMS) 

SE8=SQRT ( ( (2 . *PN) / (PN+2 . ) ) *EMS) 

WRITE  ( IW,730)  ME , NA , CH 
WRITE  ( IW, 1030) 

WRITE  ( IW, 1040) 

DO  340  1=1, NP 

WRITE  ( IW, 1050)  I,G(I) ,VG(I) ,VS(I) 

CONTINUE 
WRITE  ( IW, 1060) 

WRITE  ( IW, 1070) 

DO  350  1=1, NP 

WRITE  ( IW , 1080)  I , (S ( I , J) , J=1 , I) 

CONTINUE 

WRITE  ( IW , 1090)  SE4 , SE5 , SE6 , SE7 , SE8 
RETURN 
CONTINUE 

METHOD  3  ANALYSIS  OF  VARIANCE 
IT=NP**2-NP-1 
TI  =  IT 

IF  (NQ.LE. 1)  GO  TO  390 
IBT= IT* IB 
BT I = I BT 
DO  380  K= 1 ,NQ 
"  DO  370  I = 1 , NP 
DO  370  J= 1 , NP 
W(I, J)=A(I, J,K)*DSQRT(SN) 

CONTINUE 

CALL  DAV3 ( NP , SST , SSA , SSD , SSR , SSMR ,  SSOR ) 

BT=BT+SST 
BA=BA+SSA 
BD=BD+SSD 
BR=BR+SSR 
BMR=BMR+SSMR 


340 


350 


360 


SE2.SE3 
(1,1=1 ,NP) 


370 


BOR=BOR+SSOR 
380  CONTINUE 
390  CONTINUE 

DO  400  1=1, NP 
DO  400  J= 1 ,NP 
W(I, J)=0.0 
DO  400  K= 1 ,NQ 

W(I, J)=W( I , J)+A(I, J,K)*DSQRT(SN)/DSQRT(RN) 
400  CONTINUE 


I  (I)2»,«D0V.U)0.t  <««  m 


CALL  DAV3  (NP  »SST ,  SSA ,  SSD ,  SSR ,  SSMR ,  SSOR) 

BT=BT-SST 

BA=BA-SSA 

BD=BD-SSD 

BR=BR-SSR 

BMR=BMR-SSMR 

BOR=BOR-SSOR 

SMT=SST/TI 

SMA=SSA/A I 

SMD-SSD/CDI 

SMR=SSR/D I 

SMMR=SSMR/A  I 

SMOR=SSOR/OR I 

FAT=SMT/EMS 

FAA=SMA/EMS 

FAD=SMD/EMS 

FAR=SMR/EMS 

FAMR=SMMR/EMS 

FAOR=SMOR/EMS 

PAT=FPROB( ITt IE, FAT) 

PAA=FPROB( I A, IE,FAA) 

PAD=FPROB ( ID, IE, FAD) 

PAR=FPROB( ID, IE, FAR) 

PAMR=FPROB ( I A , IE.FAMR) 

PAOR=FPROB ( IOR , IE.FAOR) 

IF  (NQ.LE. 1)  GO  TO  410 

BMT=BT/BT I 

BMA= BA/BA  I 

BMD=BD/BCDI 

BRM=BR/BDI 

BMMR=BMR/BAI 

BMOR=BOR/BOR I 

FABT =BMT/EMS 

FABA=BMA/EMS 

FABD=BMD/EMS 

FABR=BRM/EMS 

FABMR=BMMR/EMS 

FABOR=BMOR/EMS 

PABT =FPROB( IBT, IE, FABT ) 

PABA=FPROB ( I BA , I E , FABA ) 

P ABD=FPROB ( I BD , IE,FABD) 

PABR=FPROB ( I BD , IE , FABR) 

PABMR=FPROB ( I BA , IE.FABMR) 

PABOR=FPROB( I BCR, IE.FABOR) 

410  CONTINUE 

WRITE  <IW,730)  ME,NA,CH 
WRITE  ( IW ,720) 

WRITE  (IW.74Q) 

WRITE  ( IW , 1020)  IT , SST , SMT , FAT . P AT 
WRITE  ( IW,750)  IA,SSA,SMA,FAA,PAA 
WRITE  ( IW,760)  I  CD,  SSD,  S>1D,  FAD,  PAD 
WRITE  <IW,800)  ID, SSR, SMR, FAR, PAR 
WRITE  ( IW,810)  I A , SSMR , SMMR , FAMR , PAMR 
WRITE  ( IW,820)  IOR,SSOR ,SMOR,FAOR,PAOR 
IF  (NQ.LE.  1)  GO  TO  420 
WRITE  ( IW ,830)  I BT , BT , BMT , FABT , PABT 
WRITE  (IW.840)  IBA , BA , BMA , FABA , PABA 
WRITE  ( IW,850)  I BCD , BD , BMD , FABD , PABD 
WPITE  (IW,890)  I BD , BR , BRM , FABR , PABR 
WRi  >E  (IW,900)  IBA , BMR , BMMR ,  FABMR ,  PAB>1R 
WRITE  ( IW ,910)  I BOR , BOR , BMOR , F ABOR , P ABOR 
420  CONTINUE 

WRITE  ( IW ,920)  I E , ESS , EMS 

EMS=EMS/RN/SN 

WRITE  ( IW ,730)  ME , NA , CH 

WRITE  ( IW, i 160)  ( I ,  I  =  1 , NP ) 

DO  440  1=1, NP 
W(I,  I)=G.0 
X  ( I ) =0 . 0 
Y(  I) =0-0 
DO  430  J  = 1 , NP 

W( I , J) =W( I , J) /DSQRT (SN ) /DSQRT ( RN> 

X( I) =X( I ) +W( I , J) 

Y(I)=Y(I)+W(J,1) 

TOT=TOT+W ( I , J) 

430  CONTINUE 

WRITE  < IW , 1080)  I, (W(I, J) ,J=1,\P) 

440  CONTINUE 

XMU=TOT/PN/ ( PN- 1 . ) 

SE 1 =SQRT ( EMS/PN/ ( PN- 1 . ) ) 


■ 
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WRITE  ( IW,  1 130)  XMU.SEI 
METHOD  3  COMPONENTS  OF  VARIANCE 
IF  (M0.LT.2)  GO  TO  450 
WRITE  (IV, 730)  ME,NA,CH 

CALL  COV (NP ,NR,NS,MO,SMA,SMD,SMR,BMT,EMS, AI ,DI ,BTI , IE,PP, 
RETURN 

450  CONTINUE 

IF  (MO.NE.l)  GO  TO  530 

METHOD  3  GRIFFING’S  COMBINING  ABILITY  ANALYSIS 
DO  460  1=1 ,NP 
R(  I , I) =0.0 
S  ( I , I ) =0 . 0 
VS(I)=0.0 
460  CONTINUE 
IG=NP- 1 

DO  470  1  =  1,  IG 
K=I+ 1 

DO  470  J=K,NP 

S(I ,  J)  =  (W(  I ,  J)+W(J,  I)  )/2.-(X(I)+X(J)+Y(I)+Y(J)  )/2./(PN-2. 
&-1 .  )/(PN-2. ) 

R(I,J)=(W(I,J)-W(J,I))/2. 

S(J,I)=S(I,J) 

R( J, I ) =-R( I , J) 

470  CONTINUE 

DO  490  1=1 ,NP 
DO  480  J=1 ,NP 
VS(I)=VS(I)+S(I, J)**2 
480  CONTINUE 

G(I)=(PN*(X(I)+Y(I) )-2.*T0T)/2./PN/(PN-2. ) 

VG( I) =G( I)** 2-EMS* ( PN- 1 . ) /PN/ ( PN-2 . ) /2 . 

VS ( I ) =' VS ( I ) / ( PN-2 . ) /2 . -EMS* ( PN-3 . ) / ( PN-2 . ) /2 . 

490  CONTINUE 

SE2=SQRT (EMS* (PN-1 . )/2./PN/(PN-2. ) ) 

SE3=SQRT (EMS/(PN-2. ) ) 

SE4 =SQRT ( EMS* ( PN-3 . ) /2 . / ( PN- 1 . ) ) 

SE5=SQRT (EMS* (PN-3 . ) / (PN-2 . ) ) 

SEL=SQRT (EMS* (PN-4. )/(PN-2. ) ) 

SE7=SQRT (EMS/2 . ) 

WRITE  ( IW,730)  ME , NA , CH 
WRITE  ( IW , 1030) 

WRITE  ( IW, 1040) 

DO  500  1  =  1 1  NP 
WRITE  ( I W , 1 050 ) 

500  CONTINUE 

WRITE  ( IW, 1060) 

WRITE  ( IW, 1070) 

DO  510  1=1, NP 
WRITE  ( IW,  1080) 

510  CONTINUE 

WRITE  ( IW, 1 140)  SE4 , SE5 , SE6 
WRITE  (IW, 1110) 

DO  520  1=1 ,NP 
WRITE  ( IW, 1080) 

520  CONTINUE 

WRITE  ( IW ,  1 150) 

RETURN 

530  CONTINUE 

DO  540  1=1, NP 
DO  540  J  = I , NP 
W(I, J)=(W(I, J)+W(J, I))/2.0 
W( J, I) =W( I , J) 

540  CONTINUE 

GARDNER-EBERHART  CONSTANTS 
WRITE  ( IW’,730)  ME,NA  ,CH 
CALL  GARDCO(NP ,ME, IW) 

RETURN 

METHOD  4  ANALYSIS  OF  VARIANCE 
550  CONTINUE 

IT=NP* (NP-1 )/2- 1 
TI  =  IT 

IF  (NQ.LE.l)  GO  TO  580 
IBT= IT* IB 
BT I = I BT 
DO  570  K= 1 ,NQ 
DO  560  I  =  1 ,  NP 
DO  560  J= I , NP 
W( I , J) =A( I , J ,K)*DSQRT (SN) 

W(J, I) =W( I , J) 

560  CONTINUE 

CALL  D AV4 ( NP , SST , SS A , SSD ) 


I ,G( I) , VG( I) ,VS(I) 

SE2 , SE3 
(I, 1=1, NP) 

I, (S(I, J) ,J=1,I> 


(I, 1=1, NP) 

I ,  (R( I , J) , J= 1 ,NP) 
SE7 


PU,PPP) 


) +TOT/(PN 


* 


(  s>2*.  ‘i  j  >  vaa 
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BT=BT+SST 
BA=BA+SSA 
BD=BD+SSD 
570  CONTINUE 
580  CONTINUE 

DO  590  1=1 ,NP 
DO  590  J  = I , NP 
W(I, J)=0.0 
DO  590  K=1,NQ 

I ,  J)=W(I ,  J)+A(I,  J,K)*DSQRT(SN) /DSQRT (RN) 

590  CONTINUE 

CALL  DAV4(NP,SST  ,SSA,SSD) 

BT=BT-SST 

BA=BA-SSA 

BD=BD-SSD 

SMT =SST/T I 

SMA=SSA/AI 

SMD=SSD/CDI 

FAT=SMT/EMS 

FAA=SMA/EMS 

FAD=SMD/EMS 

PAT=FPROB( IT, IE, FAT) 

PAA=FPROB ( IA , IE.FAA) 

PAD=FPROB( ID, IE, FAD) 

IF  (NQ.LE. 1)  GO  TO  600 
BMT=BT/BTI 
BMA=BA/BAI 
BMD=BD/BCDI 
F  ABT = BMT/EMS 
FABA=BMA/EMS 
FABD=BMD/EMS 
PABT=FPROB( IBT, IE.FABT) 

PABA=FPROB ( I BA , I E , FABA) 

PABD=FPROB( IBD, IE.FABD) 

600  CONTINUE 

WRITE  ( IW, 730)  ME,NA, CH 
WRITE  ( IW, 720) 

WRITE  ( IW ,740) 

WRITE  ( IW, 1020)  IT, SST, SMT, FAT, PAT 
WRITE  ( IW,750)  I A , SSA , SM A , FA A , PA A 
WRITE  ( IW,760)  I CD , SSD , SMD , FAD , PAD 
IF  (NQ.LE. 1)  GO  TO  610 
WRITE  ( IW,830)  I BT , BT , BMT , FABT , P ABT 
WRITE  (IW.S40)  I BA , BA , BMA , FABA , PABA 
WRITE  ( IW,850)  I BCD , BD , BMD , FABD , P ABD 
610  CONTINUE 

WRITE  (IW,920>  IE,ESS,EMS 
WRITE  ( IW ,730)  ME , NA , CH 
WRITE  (IW, 1170)  (1,1=1  ,NP) 

DO  640  1=1, NP 
W(I, I)=0.0 
X( I) =0.0 
DO  620  J= 1 ,NP 

W ( I ,  J)  =W ( I , J ) /DSQRT (SN) /DSQRT (RN) 

X(I)=X(I)+W(I, J) 

620  CONTINUE 

DO  630  J= 1 , NP 
TOT=TOT+W ( I , J) 

630  CONTINUE 

WRITE  ( IW, 1080)  1, (W(I, J) , J=1,I) 

640  CONTINUE 

SE 1 =SQRT ( EMS*2 . /PN/ ( PN- 1 . ) ) 

XMU=2. *TOT/PN/ ( PN- 1 . ) 

WRITE  ( IW , 1 130)  XMU.SEl 
C  METHOD  4  COMPONENTS  OF  VARIANCE 
IF  (M0.LT.2)  GO  TO  650 
WRITE  ( IW,730)  ME,NA,CH 

CALL  COV (NP ,  NR,NS,MO,SMA,S>!D,0.0,BMT,EMS,AI ,DI ,BTI , IE,PP,PU,PPP) 
RETURN 

650  CONTINUE 

C  GARDNER-EBERHART  CONSTANTS 
IF  (MO.NE.0)  GO  TO  660 
WRITE  ( IW,730)  ME,NA ,CH 
CALL  GARDCO ( NP , ME , I W ) 

RETURN 

660  CONTINUE 

C  METHOD  4  GRIFFING’S  COMBINING  ABILITY  ANALYSIS 
EMS=EMS/RN/SN 
DO  670  1=1, NP 


■ 
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S(  I ,  I ) =0.0 
DO  670  J  = I ,  NT 
IF  (I.EQ.J)  GO  TO  670 

S(I,  J)=W(I, J)-(X<I)+X(I))/(PN-2.)+2.»T0T/(PN-l.)/(PN-2.) 

S(J, I)=S(I, J) 

670  CONTINUE 

DO  690  1=1, NP 
DO  680  J  = I , NP 
VS(I)=VS(I)+S(I, J)**2 
680  CONTINUE 

G(I)=(PN*X(I)-2. *TOT ) /PN/(PN-2 . ) 

VG( I)  =G( I ) *  *2-EMS* ( PN- 1 . ) /PN/ (PN-2 . ) 

VS ( I ) =  VS ( I ) / ( PN-2 . ) -EMS* ( PN-3 . ) / ( PN-2 . ) 

690  CONTINUE 

SE2=SQRT (EMS* (PN- 1 . )/PN/(PN-2. ) ) 

SE3=SQRT (EMS* 2. /(PN-2. ) ) 

SE4=SQRT (EMS* (PN-3. )/(PN-l . ) ) 

SE5=SQRT (EMS»2 . * (PN-3. )/ (PN-2. ) ) 

SE6=SQRT (EMS *2 . * (PN-4 . ) / (PN-2 . )  ) 

WRITE  ( IW,730)  ME , NA , CH 
WRITE  ( IW, 1030) 

WRITE  ( I W , 1040) 

DO  700  1= I ,NP 

WRITE  ( IW, 1050)  I , G ( I ) ,VG(I) ,VS(I) 

700  CONTINUE 

WRITE  ( I W , 1060) 

WRITE  ( IW , 1070) 

DO  710  1=1, NP 
WRITE  ( IW , 1080) 

710  CONTINUE 

WRITE  ( IW, 1 140) 

RETURN 

720  FORMAT  (’  ’ , ’ANALYSIS  OF  VARIANCE’/) 

730  FORMAT  ( ’ 1 ’ , ’METHOD’ , 12, ’  DIALLEL  ANALYSIS ’//IX, 20 A 4// IX, 10A4//) 
740  FORMAT  (’0’,’SOURCE  OF  VARIATION’ , 12X, ’DF’ ,8X, ’SUM  OF  SQUARES’ ,8X, 
&’MEAN  SQUARES’,  18X, ’FM4X, ’P’> 

750  FORMAT  ( ’0’ ,2X, ’GENERAL  (GCA) ’ , 15X, I3.2F20.4, 10X.2F15.4) 

(’0’ ,2X, ’HETEROSIS  (SCA) ’ , 13X, I3,2F20.4, 10X.2F15.4) 

( ’0’ ,4X, ’AVERAGE  HETEROSIS’ , 1 IX, ’ 1 ’ ,2F20.4, 10X,2F15.4> 

( ’0’ ,4X , ’VARIETY  HETEROSIS’ ,9X, I3,2F20.4, 10X.2F15.4) 

( ’0’ ,4X, ’SPECIFIC  HETEROSIS’ ,8X, I3,2F20.4, 10X.2F1S.4) 

( ’0’ ,2X, ’RECIPROCAL’ , 18X, I3,2F20.4, 10X,2F15.4) 

( ’0’ ,4X, ’MATERNAL  RECIPROCAL’ ,7X, I3,2F20.4, 10X.2F15.4) 

( ’  0 ’ , 4X , ’ NON-MATERNAL  REC I PROCAL ’ , 3X , 1 3 , 2F20 . 4 , 1 0X , 2F 1 5 . 4 ) 

( ’0’ , ’GENOTYPE  INTERACTION’ , 10X, I3,2F20.4, 10X,2F15.4) 

( ’0’ ,2X , ’GENERAL  INTERACTION’ ,9X, I3,2F20.4, 10X,2F15.4) 

(’0’ ,2X, ’HETEROSIS  INTERACTION’ ,7X, I3,2F20.4, 10X,2F15.4) 

, 4X , ’AVERAGE  INTERACTION’ ,7 X, I3,2F20.4, 10X,2F15.4) 

,4X, ’VARIETY  INTERACTION’ ,7X, I3,2F20.4, 10X,2F15.4) 

, 4X , ’SPECIFIC  INTERACTION’ ,6X, I3.2F20.4, 10X,2F15.4) 

( ’0’ ,2X, ’RECIPROCAL  INTERACTION’ ,6X, I3,2F20.4, 10X.2F15.4) 

( ’0’ ,4X, ’MATERNAL  INTERACTION’ ,6X, I3,2F20.4, 10X,2F15.4) 

,4X , ’NON-MATERNAL  INTERACTION’ ,2X, I3,2F20.4, 10X,2F15.4) 
, ’ ERROR ’ , 24X , 14 , 2F20 . 4//) 

,2X, ’CROSSES’ , 2 1 X , I 3 , 2F20 . 4 , 10X.2F15.4) 

( ’0’ ,2X, ’VARIETIES  VS  CROSSES’ , 10X, ’ 1’ ,2F20. 4, 10X.2F15. 4) 

(  *0’ ,2X, ’VARIETIES’ , 19X, I3,2F20.4, 10X,2F15.4) 

(  ’0’ ,4X, ’GENERAL  (GCA) ’ , 13X, I3.2F20.4, 10X,2F15.4) 

(’0’ ,4X , ’GENERAL  INTERACTION’ ,7X, I3,2F20.4, 10X,2F15.4) 

(’0’ ,4X, ’SPECIFIC  (SCA) ’ , 12X, I3,2F20.4, 10X,2F15.4) 

( ’0’ ,2X , ’CROSSES  INTERACTION’ ,9X, I3.2F20.4, 10X.2F15.4) 

,2X, ’VARIETIES  VS  CROSSES  INT. ’ ,3X, I3,2F20.4, 10X,2F15.4 


SE2.SE3 
(I , 1=1, NP) 

I, (S(I, J) , J=l,  I) 

SE4 , SE5 , SE6 


760 

770 

780 

790 

800 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
810  FORMAT 
820  FORMAT 
830  FORMAT 
840  FORMAT 
850  FORMAT 
860  FORMAT 
870  FORMAT 
880  FORMAT 
890  FORMAT 
900  FORMAT 
910  FORMAT 
920  FORMAT 
930  FORMAT 
940  FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
&) 

1010  FORMAT 
1020  FORMAT 
1030  FORMAT 
1040  FORMAT 

&7X, ’S.C.A. 
1050  FORMAT  (’ 


( ’0’ 
(’0’ 
(  ’0’ 


( ’0’ 
( ’0’ 
(’0’ 


950 

960 

970 

980 

990 

1000 


( ’0’ 

(’0’ 
(’0’ 
( ’0’ 
( ’0’ 


,2X, ’VARIETIES  INTERACTION’ ,7X , I3,2F20.4, 10X,2F15.4) 

, ’GENOTYPES’ ,21X, I3,2F20.4, 10X.2F15.4) 

, ’ GR I FF I NG ” S  COMBINING  ABILITY  ANALYSIS  (MODEL  I)’/) 

, 10X, ’PARENT’ ,5X, ’G.C. A.  EFFECTS’ ,7X, ’G.C. A.  VARIANCE’ 
VARIANCE’/) 

12X , I2,9X,F10.4,2(10X,F12.4)) 


1060  FORMAT  ( ’0’ , 10X, ’S.E.  OF  THE  G.C. A.  EFFECTS  -’ ,20X,F10.4/1 IX, ’S.E. 

&  OF  THE  DIFFERENCES  BETWEEN  G.C. A.  EFFECTS  -\F10.4//) 

1070  FORMAT  ( ’0’ , 10X , ’S.C.A .  EFFECTS ’//6X, ’PARENTS’ ,4X, 15,91 10) 

1080  FORMAT  ( *0’ , 1 IX, 12, 2X, 10F1O.4) 

1090  FORMAT  ( ’0’ , 10X, ’S.E.  OF  THE  PARENTAL  S.C.A.  EFFECTS  -’ ,35X,F10.4/ 
&11X ,’S.E.  OF  THE  CROSS  S.C.A.  EFFECTS  -’ ,38X,F10.4/1 IX, ’S.E.  OF  TH 
&E  PARENTAL  DIFFERENCES  (S.C.A.  EFFECTS)  -’ ,21X,F10.4/1 IX, ’S.E.  OF 
&THE  COMMON  PARENT  CROSS  DIFFERENCES  (S.C.A.  EFFECTS)  -’ , 10X.F10.4/ 
&1 IX, ’S.E.  OF  THE  DIFFERENT  PARENT  CROSS  DIFFERENCES  (S.C.A.  EFFECT 
&S)  -’ ,7X,F10.4) 

1100  FORMAT  (’  ’ , 10X, ’S.E.  OF  THE  PARENT  MINUS  PROGENY  DIFFERENCES  (S.C 


. 
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&.A.  EFFECTS)  -’ ,9X,F10.4/1 IX, ’S.E.  OF  THE  PARENT  MINUS  NON-PROGENY 
&  DIFFERENCES  (S.C.A.  EFFECTS)  -\SX,F10.4) 

1110  FORMAT  ( ’0*  , 10X, ’RECIPROCAL  EFFECTS V/25X , *  MALE  PARENT’/’  FEMALE 
&PARENT’ , 18,91 10) 

1120  FORMAT  (’0M0X/S.E.  OF  THE  RECIPROCAL  EFFECTS  -*  ,40X,F  10.4/1  IX,  ’ 
&S.E.  OF  THE  DIFFERENCES  BETWEEN  RECIPROCAL  EFFECTS  ,20X,F10.4//) 
1130  FORMAT  (’0’ , ’POPULATION  MEAN  -  * ,F10.4,5X, ’S.E.  -  \F10.4//) 

1140  FORMAT  (’0’ , 10X, ’S.E.  OF  S.C.A.  EFFECTS  -  ’ ,48X,F10.4/1 IX, ’S.E.  OF 
&  THE  COMMON  PARENT  CROSS  DIFFERENCES  (S.C.A.  EFFECTS)  -  ’ , 10X, 
&F10.4/1 IX, ’S.E.  OF  THE  DIFFERENT  PARENT  CROSS  DIFFERENCES  (S.C.A. 
&EFFECTS)  -  ’  7X  F10.4) 

1 1S0  FORMAT  ( ’0* , ’S.E.  OF  THE  RECIPROCAL  EFFECTS  -  ’,40X,F10.4) 

1160  FORMAT  (’0’ , ’DIALLEL  MEANSV/25X, ’MALE  PARENT’/’  FEMALE  PARENT’, 18 
& , 91 10) 

1170  FORMAT  ( ’0’ , ’DIALLEL  MEANS ’//6X, ’PARENTS’ ,4X, 15,91 10) 

END 

SUBROUT  I NE  DAV 1 ( NP , SST , SSA , SSD , SSDD , SSPD , SSCD , SSR , SSMR , SSOR) 
DIMENSION  X( 10) , Y ( 10) 

DOUBLE  PREC I S I ON  PN , X , Y , W , SST , SSA , SSD , SSDD , SSPD , SSCD , SSR , SSMR , SSOR 
& , X 1 , X2 , X4 , X6 , X8 , H2 , H4 , H6 , TO , SP , SSP 
COMMON  W( 10, 10) 

PN=NP 

IA=NP-1 

X 1  =0 . 0 

X2=0.0 

X4=0.0 

X6=0.0 

X8=0.0 

H2=0.0 

H4=0.0 


H6=0.0 
TO=0.0 
SF  =0.0 
SSP =0.0 


DO  10  1=1, NP 

X(I)=0.0 

Y(I)=0.0 


DO  10  J=1 ,NP 
X1=X1+W(I, J)**2 
X( I) =X( I ) +W( J , I) 

Y( I ) =Y( I) +W( I , J) 

10  CONTINUE 
DO  20  1=1, NP 
X2=X2+ (X ( I ) +Y ( I ) ) e®2 
X4=X4+ (X ( I ) -Y ( I ) ) ®*2 
X6=X6+ (X( I ) +Y( I) -2. *W( I , I) ) *®2 
XS=X8+ (X( I ) +Y ( I ) -PN*W (1,1) )**2 


TO=TO+X( I ) 

SP=SP+W( 1,1) 

SSP=SSP+W(I, I)**2 
DO  20  J= 1 ,NP 

H2=H2+ (W(I,J)+W(J,I))**2 
20  CONTINUE 
DO  30  1=1, I A 
DO  30  K= I , I A 
J=K+ 1 

H4=H4+ (W( I , J) -W ( J , I ) ) **2 
H6=H6+ (W(I , J) +W( J, I) ) **2 


30  CONTINUE 

SST=X1-T0«*2/PN/PN 

SSA= ( X2/2 . /PN) - ( 2 . *T0**2/PN/PN) 

SSD= (H2/4 . ) - (X2/2 ./PN) + (T0**2/PN/PN) 

SSDD= ( TO- ( PN*SP ) ) *  *2/ ( PN- 1 . ) /PN**2 

SSPD= (X8- ( ( (2. «TO-PN*SP ) *  *2) /PN ) ) /PN/ (PN-2 . ) 

SSCD= ( H6/2 . ) - ( X6/2 . / ( PN-2 . ) )  +  ( (TO-SP ) ^  *2/ (PN-2 . ) / (PN- 1 


)) 


SSR=H4/2. 

SSMR=X4/2./PN 

SSOR=SSR-SSMR 


RETURN 

FND 

SUBROUT  I NE  DAV2 ( NP , ME , SST , SSA , SSD , SSDD , SSPD , SSCD ) 

DOUBlI^RECIS^ON  PN,X,W, SST, SSA, SSD, SSDD, SSPD, SSCD, SSAA.SP, TO, X2, 
&X4 , X6 , X8 , H2 , SSP 
COMMON  W( 10, 10) 

PN=NP 


SP=0.0 

TO=0.0 

X2=0.0 


X4=0.0 


. 

- 


. 

- 


. 
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X6=0.0 
XS=0.0 
H2=0.0 
SSP=0.0 
DO  20  1=1, NP 
X( I) =0.0 
SP=SP+W( I , I) 

SSP=SSP+W( I , I)##2 
DO  10  J= 1 ,KP 
X(I)=X(I)+W(I, J) 

10  CONTINUE 

X2=X2+  (X(I)+W(I ,  I ) )**2 
X4=X4+ (X(I)-W(I,!)) *92 

X6=X6+ (X(I)+W(I,I)-( (PN+2. )/2. )*W( I , I) )**2 
DO  20  J=I ,NP 
H2=H2+W( I , J) *#2 
TO=TO+W( I , J) 

20  CONTINUE 

SST=H2-2. *T0**2/PN/(PN+ 1 . ) 

SSA= (X2-4 . *TO*  *  2/PN ) / ( PN+2 . ) 

SSAA=X4/ (PN-2 . ) -4.* (TO-SP) ** 2/PN/ (PN-2. ) 

SSD=H2-X2/ (PN+2. ) +2. *T0**2/(PN+ 1 . )/(PN+2. ) 

SSDD= ( TO- (PN+ 1 . )/2. ®SP)**2*4. /PN/ ( PN*  «2- 1 . ) 

SSPD=4 . * (X6- ( (2 . *TO- (PN+2 . ) /2 . *SP) «*2/PN) ) / (PN**2-4 . ) 

SSCD= (H2-SSP) -X4/CPN-2. ) +2. * (T0-SP)**2/(PN- 1 . )/(PN-2. ) 

IF  (ME.NE.0)  RETURN 

SSPD=SSPD+SSA-SSAA 

SSA=SSAA 

SSD=SSA+SSCD 

RETURN 

END 

SUBROUT  I NE  DA V3 ( NP , SST , SSA , SSD , SSR , SSMR , SSOR ) 

DIMENSION  X( 10) , Y ( 10) 

DOUBLE  PRECISION  PN,X, Y,W,TO,SST,SSA,SSD,SSR,SSMR,SSOR,Xl ,X2,X4,H2 
&,H4 

COMMON  W ( 10 , 10) 

PN=NP 

TO=0 . 0 

X1=0.0 

X2=0.0 

X4=0.0 

B2=0.0 

H4=0.0 

DO  20  I = 1 , NP 
W(I,  I)=0.0 
X( I ) =0.0 
Y( I ) =0.0 
DO  10  J= 1 ,NP 
X1=X1+W(I, J)**2 
X(I) =X( I) +W( J , I) 

Y( I)=Y( I) +W( I , J) 

10  CONTINUE 
TO=TO+X( I ) 

X2=X2+ (X(I)+Y(I))*®2 
X4=X4+ (X(I)-Y(I))**2 
DO  20  J  = I , NP 

H2=H2+ (W(I,J)+W(JII))#*2 
H4=H4+ (W( I , J)-W( J, I) )*®2 
20  CONTINUE 

SST =X 1 -TO* *2/PN/ ( PN- 1 . ) 

SSA=X2/2 . / ( PN-2 . ) -2 . *TO* • 2/PN/ ( PN-2 . ) 

SSD=H2/2 . -X2/2 ./ (PN-2 . ) +T0**2/ (PN-1 . )/(PN-2. ) 

SSR=H4/2 . 

SSMR=X4/2./PN 

SSOR=SSR-SSMR 

RETURN 

END 

SUBROUTINE  DAV4(NP,SST,SSA,SSD) 

DIMENSION  X(10) 

DOUBLE  PRECISION  PN,X,W,SST,SSA,SSD,X2,H2,T0 
COMMON  W( 10, 10) 

PN=NP 

X2=0.0 

H2=0.0 

TO=0.0 

DO  20  I = 1 , NP 

w  ( r ,  n=0.0 
x ( I ) =0 . 0 
DO  10  J  =  1 ,  NP 
X(I)=X(I)+V(I,J> 


. 
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10  CONTINUE 

X2=X2+X( I)**2 
DO  20  J  = I , NP 
H2=H2+W ( I ,  J ) *#2 
T0=T0+W( I ,  J) 

20  CONTINUE 

SST=H2-2. *T0**2/PN/(PN-1 . ) 

SSA=X2/ ( PN-2 . ) -4 . *TO*  *2/PN/ ( PN-2 . ) 

SSD=H2-X2/(PN-2. )+2.*T0**2/(PN-l . )/(PN-2. ) 

RETURN 

END 

SUBROUT I NE  GARDCO ( NP , ME , I W ) 

GARDNER-EBERHART  DIALLEL  CONSTANTS  SUBROUTINE 
NP  =  NUMBER  OF  PARENTS 
ME  =  DIALLEL  METHOD  NUMBER 

0  =  PARENTS  AND  FI’S  (GARDNER-EBERHART  ANALYSIS  III) 

1  =  PARENTS, FI’S, AND  RECIPROCAL  FI’S 

2  =  PARENTS  AND  FI’S 

3  =  FI’S  AND  RECIPROCAL  FI’S 

4  =  FI’S 

IW  =  LOGICAL  I/O  UNIT  FOR  OUTPUT  STATEMENTS  :  WRITE( IW,FMT) . . . 
W  =  (INPUT)  DIALLEL  TABLE  (NP  X  NP) 


NA 


10 


20 

30 


40 

50 


60 


70 


80 

90 

100 


DIMENSION  X( 10) , Y( 10) 

DOUBLE  PRECISION  VM,CM, AH,X, Y,W 
COMMON  W  ( 10 , 10) 

NA=0 

IF  (ME.EQ.0)  NA= 1 
=  TYPE  OF  ANALYSIS 
0  =  ANALYSIS  II 
1  =  ANALYSIS  III 
WRITE  ( IW,90) 

IF  (NA.EQ.0)  WRITE  (IW,160) 

IF  (NA.NE.0)  WRITE  (IW,170) 

VM=0.0 

CM=0.0 

DO  10  1=1, NP 

VM=VM+W( I , I)/FLOAT(NP) 

DO  10  J= 1 ,NP 
IF  (I.EQ.J)  GO  TO  10 
CM=CM+W ( I , J) /FLOAT (NP* (NP- 1 ) ) 
CONTINUE 


AH=CM-VM 

IF  (ME.GT.2)  GO  TO  30 
WRITE  ( IW, 100)  VM , CM , AH 
DO  20  1=1, NP 
X(I)=W(I,  D-VM 
WRITE  ( IW, 120)  I,X(I) 
CONTINUE 
GO  TO  50 
CONTINUE 


VM=0.0 


AH=0.0 

WRITE  ( IW, 1 10)  CM 

DO  40  1=1, NP 

X( I) =0.0 

CONTINUE 

DO  60  1=1, NP 

Y( I ) =0.0 

DO  60  J= 1 ,NP 

IF  (I.EQ.J)  GO  TO  60 

Y(I)=Y(I)+(W(I , J) -CM- (X(I)+X(J))/2. ) /FLOAT (NP-1) 

CONTINUE 

WRITE  ( IW, 130) 

DO  70  1=1, NP 

Y( I ) =Y( I ) +Y( I) /FLOAT (NP-2) 


VH=Y ( I) 

IF  (NA.NE.0)  VH=Y ( I ) +X( I )/2. 
WRITE  ( IW, 120)  I , VH 
CONTINUE 


WRITE  ( IW, 140) 

DO  80  1=1, NP 
DO  80  J  = I , NP 
IF  (I.EQ.J)  GO  TO  80 

W(I, J)=W(I, J)-CM-(X(I)+X(J))/2.-Y(I)-Y(J) 
WRITE  (IW, 150)  I , J,W( I , J) 

CONTINUE 

RETURN 

FORMAT  ( ’ 0GARDNER-EBERHART  CONSTANTS ’ > 
FORMAT  ( ’ 0VAR I ET Y  MEAN ’ , F 1 5 . 4/ ’  CROSS  MEAN 


.F17.4/’ 


AVERAGE  HETERO 


. 


■ 


c 


c 


&SIS’ ,F10.4//’  VARIETY  CONSTANTS ’/’  PARENT* ,5X, ’CONSTANT* ) 

110  FORMAT  (’0CROSS  MEAN’ ,F17.4) 

120  FORMAT  (’  * ,15, F 14. 4) 

130  FORMAT  ( ’©VARIETY  HETEROSIS  CONSTANTS  (II)  OR  GENERAL  (III)*/’  PAR 
SENT’ ,5X, ’CONSTANT’) 

140  FORMAT  ( 1 ©SPECIFIC  HETEROSIS  CONSTANTS’/*  CROSS ’, 5X CONSTANT ’ ) 

150  FORMAT  (*  ’,12,’  X  \I2,F13.4) 

160  FORMAT  (’0 ANALYSIS  II’/) 

170  FORMAT  (’ ©ANALYSIS  III’/) 

END 

SUBROUTINE  COV(NP,NR,NS,MO,SMA,SMD,SMR,BMT,EMS,AI ,DI ,BTI , IE,PP,PU, 
&PPP) 

PN=NP 

RN=NR 

SN=NS 

MODEL  A  DESIGN  COMPONENTS  OF  VARIATION 
CVA=0 .0 
CVD=0.0 
CVR=0 .0 
VCA=0.0 
VCD=0.0 
VCR=0.0 
HF=0.0 

IF  (SMA.GT.EMS)  CVA= (SMA-EMS)/2./PN/RN/SN 

IF  (SMD.GT.EMS)  CVD= (SMD-EMS)*PN®*2/2./(PN**2-PN+ 1 . )/RN/SN 

IF  (SMR.GT.EMS)  CVR= (SMR-EMS)/RN/SN/2. 

VE=SQRT (2 . *EMS**2/IE) 

IF  (SMA.GT.EMS)  VCA=SQRT(2.*CVA**2/AI) 

IF  (SMD.GT.EMS)  VCD=SQRT(2.*CVD*»2/DI) 

IF  (SMR.GT.EMS)  VCR=SQRT (2. *CVR**2/DI ) 

AV=PP*CVA 

DV=PU*CVD+CVR 

GV=DV+AV 

PV=GV+EMS 

VAV=PP*PP*VCA 

VDV=PPP*VCD+VCR 

VGV=VDV+VAV 

VPV=VGV+VE 

IHN=AV/PV*100 

IHB=GV/PV* 100 

IF  ( DV . GT . 0 . 0 . AND . AV . GT . 0 . 0 )  HF=SQRT (DV/AV) 

IF  (M0.EQ.2)  GO  TO  10 
WRITE  ( IW ,20) 

WRITE  (IW.30)  CVA , VCA , CVD , VCD 
WRITE  (IW  40)  CVR  VCR 

write  ( iw' 50)  ems!ve,av,vav,dv,vdv,gv,vgv,pv,vpv,ems,ve, IHN,IHB,HF 
RETURN 
10  CONTINUE 

MODEL  II  DESIGN  COMPONENTS  OF  VARIATION 
CVA=0.0 
VCA=0 .0 
CVI=0.0 
VCVI=0.0 


HF=0. 0 

IF  (SMA.GT.SMD)  CVA= (SMA-SMD1/2./PN/RN/SN 
IF  (SMA.GT.SMD)  VCA=SQRT(2.*CVA*«2/AI) 

IF  (BMT . GT . EMS . AND . NR . GT . 1 . AND . NS . GT . 1 )  CVI= (BMT-EMS)/SN 

IF  ( BMT . GT . EMS . AND . NR . GT . 1 . AND . NS . GT . 1 )  VCVI=SQRT(2.»CVI*®2/BTI) 

AV=PP*CVA 


GV=DV+AV 
PV=GV+EMS+CVI 
VAV=PP*PP*VCA 
VGV=VDV+VAV 
VPV=VGV+ VE+VCV I 
IHN=AV/PV® 100 
IHB=GV/PV® 100 

IF  ( DV . GT . 0 . 0 . AND . AV . GT . 0 . 0)  HF=SQRT (DV/AV) 


WRITE  ( IW,60) 

WRITE  ( IW ,30)  CVA, VCA, CVD, VCD 
WRITE  ( IW , 40)  CVR, VCR 

WRITE  ( I W , 50)  EMS , VE , AV , VAV , DV , VDV , GV , VGV , PV , VPV , EMS , VE, IHN , IHB , HF 
RETURN 


20  FORMAT  ( ’0’ , ’COMPONENTS  OF  GENETIC  VARIANCE  (MODEL  A)’/) 

30  FORMAT  (’0\9X, ’GENERAL’, F20. 4,*  .F15.4//10X, ’HETEROSIS’ ,F18. 4, 

+- ’ ,F15 . 4/) 

40  FORMAT  (’  ’ ,9X, ’RECIPROCAL’ ,F17.4, ’  ,F15.4/) 

50  FORMAT  (’  ’ ,9X, ’ERROR’ ,F22.4, ’  ,F15.4////10X, ’ADDITIVE’ ,F18. 4, ’ 

&  +-’ ,F15.4//10X, ’DOMINANCE’ ,F17.4, ’  ,F15.4//10X, ’GENETIC’ ,F19.4 

, F 1 5 . 4// 1 0X , ’PHENOTYPIC’ ,F16. 4, ’  +-\F15.4//10X, ’ENVIRONMENT 
ML’  .F13.4,  ’  .F15.4////10X,  ’HERITABILITYV/15X,  ’NARROW’  ,  1 10,  ’%’/ 


■ 


A/1SX,  *  BROAD  *  ,  1 1 1,  ’%’///10X,  ’HETEROSIS  FACTOR  ( SQRT  ( DOM  I  NANCE/ADD  I T 
AIVE)) \F10.4/> 

60  FORMAT  <’ 0  \  ’ COMPONENTS  OF  GENETIC  VARIANCE  (MODEL  II)’/) 

END 

FUNCTION  FPROB(M,N,X) 

COMPUTES  PROBABILITY  OF  F-STATISTIC  X  WITH  M  AND  N  DEGREES  OF 
FREEDOM  IN  NUMERATOR  AND  DENOMINATOR  RESPECTIVELY. 

INTEGER  A,B 
FPROB= 1 .0 

IF  (X.LE.0. .OR.M.LE.0.OR.N.LE.0)  RETURN 

A=2* (M/2)“M+2 

B=2* (N/2) -N+2 

W=X*M/N 

Z=l./(1.+W) 

IF  (A.NE. 1)  GO  TO  20 
IF  (B.NE. 1)  GO  TO  10 
P=SQRT (W) 

Y= 0.3 183099 

D=Y*Z/P 

P=2. *Y*ATAN (P) 

GO  TO  40 
10  P=SQRT(W*Z) 

D=0.5«P*Z/W 
GO  TO  40 

20  IF  (B.NE. 1)  GO  TO  30 
P=SQRT(Z) 

D=0.5*Z*P 
P=l.-P 
GO  TO  40 
30  D=Z*Z 
P=W#Z 

40  Y=2.*W/Z 
J=B+2 

IF  ( J.GT.N)  GO  TO  60 
DO  50  K=J,N,2 
D=(1.+A/(K-2.))*D«Z 
IF  (A.EQ. 1)  P=P+D*Y/(K-1) 

IF  (A.NE. 1)  P=(P+W)*Z 
50  CONTINUE 
60  Y=W*Z 
Z=2./Z 
B=N-2 
I=A+2 

IF  (I.GT.M)  GO  TO  80 
DO  70  K=I ,M,2 
J=K+B 

D=Y*D*J/(K-2) 

P=P-Z«D/J 
70  CONTINUE 
80  FPROB=l.-P 
RETURN 
END 


. 
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Appendix  C*  Hayman's  analysis  of  population 
subroutine* 


1:  HAYMAN  subroutine  listing* 


2;  HAYMAN  sample  output* 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


10 

c 


20 

30 

c 


40 


50 

C 


C 

C 

C 


60 

70 


80 


90 


subroutine:  hapmuw,x,v,pms,  ie,ems,name,ch,nc,nv> 

HAYMAN  S  ANALYSIS  OF  POPULATION  MEANS 
VERSION  1.0,  JUNE  1979 


IW 
X 
V 

PMS 

IE 

EMS 

NAME 

CH 

NC 

NV 


I/O  UNIT  FOR  WRITE  CONTROL  (I.E.  WRITE( IW, . . . ) ) 

VECTOR  OF  MEANS  (ORDER  =  PI ,P2,F1 ,F2,B1 ,B2) 

VECTOR  OF  ST.  DEVIATIONS  (ORDER  =  P 1 ,P2,F1 ,F2,B1 ,B2) 
POPULATIONS  MEAN  SQUARE 
ERROR  DEGREES  OF  FREEDOM 
ERROR  MEAN  SQUARE 

NAME  OF  THE  EXPERIMENT  (ALPHANUMERIC  20A4) 

CHARACTER  IDENTIFICATION  (ALPHANUMERIC  20A4) 

MAXIMUM  NUMBER  OF  ’WORDS’  IN  ’NAME’  AND  ’CH’ 

0  =  UNWEIGHTED,  1  =  VARIANCE  WEIGHTED 

REQUIRES  SUBROGRAM  -  FPROB 

DIMENSION  X( 1) ,V( 1) , NAME( 1) ,CH( 1) 

DIMENSION  A(6,6) ,Q(6) ,P(6> ,R(6) ,Y(6) 

REAL  B (6, 6) /l. 0,1. 0,0. 0,1. 0,0. 0,0. 0,1. 0,-1. 0,0. 0,1. 0,0. 0,0. 0,1.0, 
&0. 0,1. 0,0. 0,0. 0,1. 0,1. 0,0. 0,0. 5, 0.0, 0.0, 0.25, 1.0, 0.5, 0.5, 0.25, 0.25 
&,0.25, 1 .0,-0.5,0.5,0.25,-0.25,0.25/ 

REAL  P0P(6)/’P1’ , ’ P2 ’ , ’ F 1 ’ , ’ F2 ’ , ’ BC 1 ’ , ’BC2 V 
PSS=PMS*5. 

ESS=EMS*IE 
DO  10  1=1,6 
B(  1 , 1) =X( I) 

WEIGHTING  BY  STANDARD  ERROR 
IF  (NV.EQ.0)  GO  TO  30 
DO  20  J= 1 ,6 
DO  20  1=1,6 
B(I,J)=B(I,J)/V(J) 

CONTINUE 

CONTINUE 

COMPUTE  CORRELATIONS 
DO  50  1=1,6 
DO  50  J= I , 6 
TX=0.0 
TY=0.0 
SX=0.0 
SY=0.0 


SXY=0.0 
A( I , J) =0.0 
DO  40  K= 1 ,6 
TX=TX+B( I ,K) 

TY=TY+B( J,K) 

SX=SX+B(I,K)**2 
SY=SY+B(J,K)**2 
SXY=SXY+B( I ,K)*B( J,K) 

CONTINUE 

SMX=SX-TX**2/6. 

SMY=SY-TY#*2/6. 

SMXY =SXY-TX*TY/6 . 

A( I , J) =SMXY/SQRT (SMX*SMY) 

A( J , I ) =A( I , J) 

CONTINUE 

COMPUTE  MULTIPLE-R  SQUARED  AND  SUMS  OF  SQUARES 
AUP=  1 .0 
R(  1)  =  1 .0 
DO  110  L=2,6 
K=L 

IF  (L.LT.4)  GO  TO  70 

FINDS  AND  ENTERS  THE  INTERACTION  COMPONENT  WITH  THE 
LARGEST  PARTIAL  CORRELATION  WHEREAS  THE  MAIN  EFFECTS 
ARE  ENTERED  FIRST  DIRECTLY 
AM=0.0 
DO  60  M=4,6 

VC=(A( 1 ,M)*A(M, 1) )/A(M,M) 

IF  (VC.LE.AM)  GO  TO  60 

K=M 

AM=VC 

CONTINUE 

PE=A(K,K) 

DO  80  1=1,6 
Q(  I) =A( I ,K) 

P ( I) =A(K, I) 

CONTINUE 
DO  90  1=1,6 
DO  90  J= 1 ,6 

A(I, J)=A(I, J)-(Q(I)#P(J))/PE 

CONTINUE 

DO  100  1=1,6 

A(K, I ) =P< I )/PE 


100  CONTINUE 

R(K)=AUP-A( 1,1) 

Y (K) =PSS*R(K) 

AUP=A( 1,1) 

110  CONTINUE 
C  COMPUTE  ANOVA  COMPONENTS 
EPSS=Y (4) +Y (5) +Y (6) 

ERSQ=R(4) +R(5) +R(6) 

EPMS=EPSS/3. 

PF=PMS/EMS 
PP=FPR0B(5, IE,PF) 

AF=Y (2)/EMS 
AP=FPROB( 1 , IE,AF) 

DF=Y (3)/EMS 
DP=FPROB( 1 , IE,DF) 

EPF=EPMS/EMS 
EPP=FPR0B(3, IE, EPF) 

AAF=Y (4)/EMS 
AAP=FPROB( 1 , IE,AAF) 

ADF=Y (5)/EMS 
ADP=FPROB( 1 , IE,ADF) 

DDF=Y (6)/EMS 
DDP=FPROB ( 1 , IE,DDF) 

WRITE  ( IW, 130) 

WRITE  ( IW, 140)  (NAME(I) , 1=1 ,NC) 
WRITE  ( IW, 140)  <CH(I) , 1=1, NC) 

IF  (NV.EQ.l)  WRITE  (IW,150) 

C  PRINT  MEANS 

WRITE  ( IW, 160) 

DO  120  1=1,6 

WRITE  ( IW, 170)  POP(I),X(I) 

IF  (NV.EQ.0)  GO  TO  120 
WRITE  (IW, 180)  V(I) 

120  CONTINUE 
C  PRINT  ANOVA  TABLE 
WRITE  (IW, 190) 

WRITE  ( IW,200) 

WRITE  (IW.210) 

WRITE  ( IW,220) 

WRITE  ( IW,230) 

WRITE  ( IW,240) 

WRITE  ( IW,250) 

WRITE  ( IW,260) 

WRITE  ( IW,270) 


PSS , PMS , PF , PP , R ( 1 ) 

Y (2) , Y (2) , AF, AP ,R(2) 

Y (3) , Y (3) , DF,DP,R(3) 
EPSS , EPMS , EPF , EPP , ERSQ 
Y (4) , Y (4) , AAF, AAP ,R(4) 
Y (5)  , Y (5)  ,  ADF ,  ADP , R  (5) 
Y (6) , Y(6) ,DDF,DDP,R(6) 
I E , ESS , EMS 


130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 


C 

C 


RETURN 

FORMAT  ( ’ 1HAYMANS  ANALYSIS  OF  POPULATION  MEANS’/1 
&EAN  EFFECTS)’) 

FORMAT  (’0\20A4) 

FORMAT  (’©WEIGHTED  (STANDARD  ERROR) 

FORMAT  (’-POPULATION  MEANS’/) 

FORMAT  (’  ’ ,3X,A4,F10.4) 

FORMAT  ( ’+’ ,20X, ’+-’ ,F10.4) 

FORMAT  (’-ANALYSIS  OF  VARIANCE 


ANALYSIS’) 


&MS’ , 19X, ’F’ , 19X, ’P’ , 18X , ’ RSQ ’ ) 


’//’  SOURCE’ ,8X, ’DF1 


(’0 

(’0 

(’0 

(’0 

(’0 

(’0 


ADDITIVE’ ,5X, 
DOMINANCE’ ,4X, 
EPISTASIS’ , 4X, ’3 
ADD-ADD’ ,4X, 
ADD-DOM’ ,4X, 
DOM-DOM’ ,4X, 


’ 1’ ,5F20.4) 

’ 1 ’ , 5F20.4) 


,5F20.4) 
’ 1 ’ , 5F20.4) 
’ 1 ’ , 5F20.4) 
’ 1 ’ , 5F20. 4) 


FORMAT  ( ’ 0POPULAT IONS’ ,4X, ’5’ , 5F20 . 4) 

FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 

FORMAT  ( ’ 0ERROR ’ , 6X , 15 , 2F20 . 4) 

END 

FUNCTION  FPROB(M,N,X) 

COMPUTES  PROBABILITY  OF  F-STATISTIC  X  WITH  M  AND  N 
FREEDOM  IN  NUMERATOR  AND  DENOMINATOR  RESPECTIVELY. 
INTEGER  A,B 
FPROB=  1 . 0 

IF  (X.LE.0. .OR.M.LE.0.OR.N.LE.0)  RETURN 

A=2* (M/2) -M+2 

B=2* (N/2)-N+2 

W=X*M/N 

Z=l./(1.+W) 

IF  (A.NE. 1)  GO  TO  20 
IF  (B.NE.l)  GO  TO  10 
P=SQRT (W) 

Y=0. 3 183099 

D=Y*Z/P 

P=2.*Y*ATAN(P) 

GO  TO  40 


(ADJUSTED  FOR  M 


, 14X, ’SS’ , 18X, ’ 


DEGREES  OF 


P=SQRT(W*Z) 

D=0.5*P*Z/W 
GO  TO  40 

IF  (B.NE. 1)  GO  TO  30 
P=SQRT (Z) 

D=0.5*Z*P 

P=l.-P 

GO  TO  40 

D=Z*Z 

P=W*Z 

Y=2.*W/Z 

J=B+2 

IF  ( J.GT.N)  GO  TO  60 

DO  50  K= J ,N,2 

D= ( 1 . +A/(K-2. ) ) *D*Z 

IF  (A.EQ.l)  P=P+D®Y/(K-I) 

IF  (A.NE.l)  P=(P+W)*Z 

CONTINUE 

Y=W*Z 

Z=2./Z 

B=N-2 

I  =A+2 

IF  (I.GT.M)  GO  TO  80 
DO  70  K=I ,M,2 
J=K+B 

D=Y*D® J/<K-2) 

P=P-Z*D/J 
CONTINUE 
FPROB= 1 . -P 
RETURN 
END 
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Appendix  D.  Publications. 
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STRAIN-SPECIFIC  INEFFECTIVE 
NODULATION  OF  THE  TROPICAL  LEGUME 
DESMODIUM  BY  RHIZOBIUM 1 

B.  R.  Pinchbeck2 
ABSTRACT 

The  Rhizobiunt  strain  CBI789  originally  isolated  from 
Desmodium  sandwicense  E.  Mey  produces  a  strain- 
specific  ineffective  nodulation  response  with  several  lines 
of  D.  aparines  (Link)  D.  C.  Study  of  Fs  and  backcross 
populations  from  D.  sandwicense  x  D.  aparines  hybrids 
provided  evidence  for  segregation  of  a  dominant  gene 
at  a  single  locus  conditioning  ineffective  nodulation  in 
D.  aparines.  The  D.  aparines  lines  expressing  ineffective 
nodulation  are  used  in  commerce. 

Additional  index  words:  Genetics,  Inoculation,  N  fixa¬ 
tion. 


INEFFECTIVE  nodulation  of  soybeans  [Glycine 
max  (L.)  Merr.]  by  Rhizobium  japonicum  (Kirch- 
ner)  Buchanan  has  been  shown  to  be  controlled  by  the 
three  dominant  genes  Rj2,  Rj3,  and  Rj4  (2,  7,  8). 
Major  genes  conditioning  strain-specific  nodulation 


•This  investigation  was  supported  by  a  Ph.D.  research  grant 
from  the  International  Development  Research  Ctr.,  Ottawa, 
Canada.  Received  3  June  1977. 

•Statistical  Analyst,  Computing  Services,  Univ.  of  Alberta, 
Edmonton,  Canada  (formerly  graduate  student,  Plant  Breeding 
Institute,  Univ.  of  Sydney,  Sydney,  Australia). 

•Commonwealth  of  Australia  identification  number. 

Table  1.  Qualitative  interactions'  of  Desmodium  lines  with  two 
Rhizobium  strains. 


Effectiveness! 


Species  t 

CB627 

CB1789 

D.  sandwicense  (CPI18156) 

E 

E 

D.  sandwicense  (CPI18163) 

E 

E 

D.  sandwicense  (N993) 

E 

E 

D.  sandwicense  (CSL) 

E 

E 

D.  sandwicense  (CWF) 

E 

E 

D.  sandwicense  (CLL) 

E 

E 

D.  aparines  (CPI23189) 

E 

I 

D.  aparines  (Q8382) 

E 

I 

D.  aparines  (CPI18010) 

E 

I 

D.  aparines  (Desica-1) 

E 

E 

D.  aparines  (CPI33814) 

E 

E 

D.  pringlei  (CPI38722) 

E 

E 

D.  uncinatum  (CPI8990) 

E 

E 

D  sp.  (Ecuador) 

E 

E§ 

t  Accession  or  line  identification. 

%  E— effective 

nodulation 

(green),  I— ineffective  nodulation  (yellow-chlorotic), 
duced  chlorotic  reaction  that  appears  to  be  pH  sensitive. 

§  Strain  in- 

have  been  reported  in  Trifoliutn  pratense  L.  (6)  and 
T.  subterraneum  L.  (3).  In  preliminary  studies  with 
the  tropical  legume,  Desmodium,  an  ineffective  nodu¬ 
lation  interaction  was  found  for  several  lines  of  D. 
aparines  (Link)  D.C.  [D.  intortum  (Mill.)Urb.j  when 
inoculated  with  Rhizobium  strain  CB17893  (Table  1). 
Strain  CB627,  referred  to  in  Table  1,  is  a  commercial 
inoculum  for  Desmodium  and  effectively  induced  nod¬ 
ulation  of  all  lines  studied.  Strain  CB1789  was  isolated 
from  D.  sandwicense  E.  Mey  and  strain  CB627  was 
isolated  from  D.  aparines.  Strains  CB627  and  CB1789, 
members  of  the  cowpea  group  of  Rhizobium,  are  slow- 
growing;  however,  strain  CB627  is  relatively  faster 
growing  than  strain  CB1789  (J.M.  Vincent,  personal 
communication). 

This  study  was  undertaken  to  determine  the  mode  of 
inheritance  of  the  ineffective  nodulation  interaction 
between  D.  aparines  and  Strain  CB1789. 

MATERIALS  AND  METHODS 

Fa  populations  were  produced  from  a  D.  sandwicense  male-ster¬ 
ile  line  derived  from  CPI18156  (5)  which  was  crossed  to  the  D. 
aparines  lines  CPI23189  and  Q8382.  Backcrosses  were  obtained 
by  crossing  the  Fa  hybrid  onto  the  male-sterile  line.  Crosses  be¬ 
tween  D.  sandwicense  and  D.  aparines  produce  fully  fertile 
progeny  and  are  considered  to  be  closely  allied  species  (4). 

The  F2  and  backcross  plants  were  grown  successively  in  “modi¬ 
fied”  Leonard  jar  assemblies  (9)  under  N-free  growth  conditions 
and  inoculated  with  5  ml  of  a  suspension  of  strain  CB1789.  Un¬ 
inoculated  controls  and  parental  genotypes  were  placed  adjacent 
to  the  experiment.  Eight  weeks  following  planting,  the  plants 
were  scored  for  green  vs.  yellow-chlorotic  color  induced  by  N 
deficiency.  Seventy  randomly  chosen  F2  plants  were  removed  from 
the  Leonard  jars,  inspected  for  nodulation  response  and  replanted 
in  pots  to  check  their  response  to  N  fertilizer. 

RESULTS  AND  DISCUSSION 

The  nodulation  scores  of  the  F2  and  backcross  plants 
are  presented  in  Table  2.  The  ineffectively  nodulated 
plants  were  extremely  stunted  and  chlorotic  at  8  weeks 
and  were  similar  in  appearance  to  the  uninoculated 
and  ineffectively  nodulated  controls.  The  roots  of  the 
ineffectively  nodulated  plants  exhibited  various  de¬ 
grees  of  nodule  formation.  The  majority  of  nodules 
were  small,  whereas  a  few  were  of  normal  size.  All 
nodules  sampled  lacked  pink  pigmentation.  When 
these  ineffectively  nodulated  plants  were  transferred 
to  pots,  all  responded  to  N  fertilization. 

Chi-square  analysis  of  the  segregation  data  support 
the  hypothesis  that  the  strain-specific  ineffective  nodu¬ 
lation  interaction  between  lines  of  D.  aparines  and 
strain  CB1789  is  produced  by  a  dominant  gene  at  a 


Table  2.  Observed  nodulation  scores  and  Chi-square  probability  of  F,  and  backcross  plants  inoculated  with  strain  CB1789. 


F, 

Backcross 

Nodulation 

Chi-square 

probability 

3:1 

Nodulation 

Chi-square 

probability 

Cross 

Ineffective 

Effective 

Ineffective 

Effective 

1:1 

D.  sandwicense  (CPI  18156)  x 

D.  aparines  (CPI  23189) 

88 

36 

0.80-0.70 

35 

36 

0.95-0.90 

D.  sandwicense  (CPI  18156)  x 

D  aparines  (Q  8382) 

Deviation  x1 

Heterogeneity  x’ 

91 

179 

24 

60 

0.80-0.70 

0.99-0.95 

0.30-0.20 

41 

76 

32 

68 

0.50-0.40 

0.70-0.60 

0.60-0.50 

. 
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single  locus.  These  results  are  similar  to  those  report¬ 
ed  for  soybeans  (2,  7,  8). 

The  symbols  Ni/ni  are  suggested  for  the  designation 
of  this  ineffective  nodulation  gene.  The  Desmodium 
lines  which  are  effectively  nodulated  have  the  geno¬ 
type  ninj.  The  ineffectively  nodulated  lines  would  be 
expected  to  have  the  genotype  NjNi;  however,  further 
experimentation  would  be  required  to  test  for  func¬ 
tional  allelism  between  the  two  D.  aparines  lines  stu¬ 
died.  The  ineffectively  nodulated  D.  aparines  lines 
have  small  seeds  while  effectively  nodulated  D.  apa¬ 
rines  lines  have  larger  seeds.  The  linkage  of  the  Ni 
allele  with  small  seed  size  has  not  been  studied. 

One  ineffectively  nodulated  line  of  D.  aparines  has 
attracted  interest  as  the  commercial  cultivar  —  ‘Green- 
leaf’  (1).  Ineffective  nodulation  alleles  carried  by  this 
cultivar  may  be  at  a  disadvantage  if  strain  CB1789 
were  to  attain  a  dominant  status  in  the  field. 
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